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A.  Introduction 

It  is  well  known  that  the  maintenance  and  adaptation  of  bone  integrity  is  dependent  on  a 
complex  interaction  of  metabolic  and  environmental  factors  (mechanical  stresses, 
nutritional  status).  Unfortunately,  the  specific  relationship  between  these  factors  and  the 
biomechanical  competence  of  bone  tissue  remains  incompletely  quantified.  As  a  result, 
strategies  for  preventing  or  effectively  treating  bone  fragility  or  enhancing  general  bone 
health  are  far  from  being  optimized.  The  specific  goals  of  this  research  program  were  to 
contribute  to  two  major  objectives  in  support  of  reducing  the  incidence  of  fracture: 

a.  The  development  and  application  of  micro-imaging  and  testing  techniques 
in  animal  models  to  study  bone  structure  function  properties. 

b.  Exploring  the  influence  of  calcium  and  vitamin  D  metabolism  and 
physical  forces  on  bone  integrity. 


B.  Body 

The  final  progress  report  of  this  research  program  is  described  below.  As  will  be  gleaned 
from  the  report,  essentially  all  of  the  work  and  tasks  that  were  approved  by  the 
USAMRMC  were  accomplished,  although  some  of  the  specific  methodologic  details 
varied  from  the  original  proposal.  These  few  variations  (as  noted  within  the  text)  were 
implemented  as  a  result  of  the  ongoing  findings  of  the  study.  It  is  also  important  to  note 
that  we  undertook  and  have  nearly  completed  some  additional  studies  (funded  from  our 
own  discretionary  sources)  that  were  inspired  by  the  data  from  our  study  and  others  in  the 
field.  These  will  also  be  noted  in  the  report. 

A  summary  of  the  important  research  program  results  is  provided  below.  The  majority  of 
the  results  and  accomplishments  are  presented  in  the  form  of  figures  and  tables. 
Methodologic  details  are  freely  available  to  any  interested  parties  by  contacting  us 
directly.  While  the  body  of  the  report  has  a  summary  of  selected  research  results,  all  of 
the  data  is  presented  in  the  Appendix. 

1.  Calcuim  Intake  Studies 

These  studies  were  designed  to  test  three  primary  hypotheses  and  2  subhypotheses  and 
used  two  specific  mouse  strains  (C57BL/6J  and  C3H/HeJ): 

HI:  Reductions  in  dietary  calcium  results  in  changes  in  mechanical  properties. 

H1A:  The  effects  are  dependent  on  genetic  background 
H1B:  The  effects  are  dependent  on  dose  of  calcium 

H2:  The  effects  of  altering  calcium  is  gender  dependent. 

H3:  The  effects  of  altered  calcium  intake  is  influenced  by  intense  exercise. 
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There  were  24  experimental  groups  with  variations  in  calcium  dose  (low  intake  (0.1%), 
normal  (0.4%  to  0.6%)  and  high  (1.0%)),  gender,  and  exercise  (cage  exercise  and  intense 
treadmill  running). 

The  exercise  involved  running  on  a  treadmill  system  developed  as  part  of  the  research 
program.  Since  a  primary  aim  of  the  program  was  to  explore  the  use  of  microimaging 
and  micromechanical  testing  methods  for  evaluating  skeletal  alteration  in  mouse  models, 
the  studies  utilized  numerous  outcome  measures  including;  MicroCT,  whole  bone 
mechanical  tests,  microspecimen  preparation  and  mechanical  tests,  Raman  spectroscopy 
and  histologic  analysis.  By  way  of  summary,  the  following  numbers  of  animals  were 
completed  in  the  program: 


Femoral  Analyses 


Studies 

MicroCT  scan 

MicroCT 

analysis 

Mechanical 

testing 

Raman 

spectroscopy 

Completed 

392 

392 

368 

130 

Remaining  to 
complete 

0 

0 

0 

0 

Vertebral  Body  analyses 


Studies 

MicroCT  scan 

MicroCT 

analysis 

Mechanical 

testing 

Completed 

382 

381 

394 

Remaining  to 
complete 

0 

0 

0 

Microspecimen  production 

Approximately  100  microspecimens  were  produced  for  the  program  from  these  animals 
and  the  testing  procedures  were  calibrated  and  validated.  While  the  aim  of  exploring 
techniques  for  measuring  ECM/tissue  properties  was  completed,  no  changes  in 
fundamental  tissue  properties  were  found.  As  a  result,  no  further  testing  was  performed. 

Raman  spectroscopy 

More  than  130  specimens  were  prepared  and  analyzed  using  Raman  spectroscopy.  The 
results  of  these  tests  will  be  presented  graphically  below.  This  very  time  consuming 
procedure  demonstrated  very  valuable  methodology  that  we  have  already  employed  in 
other  studies  of  mouse  skeleton  phenotyping. 


6 


Histologic  analysis 

The  histologic  studies  were  performed  as  a  method  of  examining  the  character  of  the 
tissue  morphology  to  determine  if  any  specific  pathologic  conditions  were  induced  by  the 
exercise  or  diet  alterations.  More  than  125  bones  were  examined  histologically; 
representative  animals  from  every  experimental  group.  As  will  be  shown  below,  no 
pathologic  features  were  noted. 

Additional  studies 


As  will  be  seen  in  the  data  graphs  provided  below,  the  effect  of  exercise  was 
unexpectedly  low  in  the  femurs  of  the  animals.  This  surprising  finding  lead  us  to  explore 
other  work  designed  to  evaluate  the  effect  of  loading  on  bone  by  others  in  the  scientific 
community.  As  a  result  of  this  search,  we  found  at  least  one  study  that  found  exercise 
induced  morphologic  alterations  in  the  tibia  of  mice.  Based  on  a  hypothesis  that  there 
might  be  a  site-specific  (anatomic  site)  effect,  we  extracted  the  tibias  from  60  mice  and 
did  MicroCT  analysis  on  them.  As  will  be  described,  we  found  only  a  small  effect  of 
exercise  on  the  bone  morphology  of  the  tibia  as  well. 

Selected  Results 


A  series  of  graphical  results  are  presented  below  that  summarize  the  results  of  the  studies 
on  nutrition,  gender,  and  exercise  on  bone  properties.  A  brief  statement  is  provided  with 
each  to  summarize  the  important  attribute  of  the  findings. 
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Figure  1:  Example  of  mechanical  femoral  data.  The  results  demonstrate  no  effect  of  exercise,  although 
there  were  significant  differences  in  properties  as  a  function  of  strain  and  calcium  intake.  Similar 
effects  were  found  for  ultimate  and  yield  properties,  as  well  as  all  post-yield  data. 


Stiffness  by  Ex 


nd  Diet  (Female  C57  Only) 


;  by  Exercise  and  Diet  (Male  C57  Only) 


Figure  2:  Example  of  bone  mineral  density  data  (from  MicroCT).  Once  again,  there  was  no  effect  of 
exercise,  but  there  were  significant  differences  across  the  mouse  strains. 


BMD  by  Die 


!  (Male  C57  only) 


i  by  Diet  &  Exercise  (Female  C57  only) 
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Figure  3:  Raman  spectroscopy  data  is  illustrated.  There  is  both  a  dose  (calcium  nutrition)  and  strain 
effect  demonstrated.  Interesting,  the  direction  of  the  dose  response  is  opposite  in  the  C57  verses  the 
C3He  animals.  This  might  suggest  a  fundamental  genetic  variation  in  response  as  a  function  of  mouse 
genotype.  It  also  demonstrates  the  power  of  the  Raman  technology  to  evaluate  bone  response. 


Mineral  to  Matrix  Ratio  by  Exercise  and  Diet  (Female  C57  Only) 


Figure  4:  The  effect  of  diet  and  exercise  on  mid-diaphyseal  bending  moment  of  inertia.  There  was  a 
significant  effect  by  level  of  calcium  and  mouse  strain.  There  was  no  significant  effect  of  exercise. 


Bending  Moment  ot  Inertia  by  Diet  a 


Bending  Moment  of  Inertia  by  Diet  and 


I  Exercise  (Female  C57  Only) 


I  Exercise  (Male  C57  Only) 


Figure  5:  The  results  for  the  vertebral  bodies  paralleled  those  of  the  femur.  As  demonstrated  here, 
there  was  a  significant  effect  of  mouse  strain  and  diet,  but  no  effect  of  exercise.  As  noted  earlier,  there 
are  also  baseline  differences  between  males  and  females,  but  their  responses  are  similar.  These  effects 
were  observed  for  all  morphologic  and  mechanical  data  from  the  vertebra. 


Bone  Volume  Fraction  by  Diet  and  Excercise  (Female  C57  Only) 


Bone  Volume  Fraction  by  Diet  and  Excercise  (Male  C57  Only) 


Bone  Volume  Fraction  by  Diet  and  Excercise  (Female  C3  Only)  Bone  Volume  Fraction  bV  Die<  and  Excercise  (Male  C3  Only) 


MM 


Figure  6:  As  noted  earlier,  the  lack  of  exercise  induced  response  in  the  femurs,  stimulated  us  to 
perform  additional  studies  in  the  tibias.  We  sampled  60  animals  and  examined  their  diaphyseal  and 
metaphyseal  bone.  As  noted  in  these  figures,  the  mouse  strain  and  diet  effects  were  replicated  in  the 
tibias.  In  contrast,  there  was  a  potential  effect  of  exercise  seen  in  the  tibias,  but  the  results  were  not 
consistent.  In  the  graph  below,  it  appears  that  the  males  demonstrated  an  increase  in  bone  volume 
fraction  with  exercise,  in  contrast  to  the  females.  These  results,  however,  were  not  seen  in  many  other 
variables  and  as  a  result,  it  is  difficult  to  support  a  powerful  effect  of  exercise. 


Tibia  Metapbyieal  Bona  Volume  Fraction  by  Sei  and  Ererclat  Tibia  Metaphyw  al  Trabecular  Number  by  Sei  and  Eaerclie 

(Noraal  diet.  C57  Strain)  (Normal  0*t  CST  Strain) 
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Figure  7:  The  histologic  analysis  demonstrated  that  the  morphology,  cellularity,  and  ECM  organization 
appeared  normal  in  all  groups.  An  example  of  the  histologic  findings  is  illustrated  below. 


2.  Vitamin  D  studies 

The  purpose  of  these  studies  was  to  explore  the  effect  of  reduced  vitamin  D  on  the 
mechanical  and  morphologic  properties  of  the  bone.  The  specific  hypotheses  to  be  tested 
were: 

HI:  Reduction  in  vitamin  D  metabolism  results  in  reduction  in  mechanical 

properties. 

H2:  The  effects  of  altered  vitamin  D  metabolism  and  gender. 

H3:  The  effects  of  altered  vitamin  D  metabolism  is  independent  of  intense  exercise. 

The  study  was  designed  to  utilize  three  specific  groups  of  animals: 

1.  C57B1/6J  Mouse  (background  control) 

2.  DBP-/-  Mouse  with  normal  diet 

3.  DBP-/-  Mouse  with  vitamin  D  deficient  diet  (Purina  diet  No.  5826C-5) 

The  DBP  (vitamin  D  binding  protein  knockout)  mouse  was  described  by  our  consulting 
investigator,  Dr.  Nancy  Cooke  and  has  already  been  shown  to  demonstrate  significant 
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bone  alterations  when  the  vitamin  D  binding  protein  deficiency  is  combined  with  a 
vitamin  D  deficient  diet. 


Progress/Results 

As  noted  in  several  of  our  yearly  progress  reports,  the  breeding  of  these  mice  was 
difficult  and  we  had  a  viral  infection  outbreak  during  the  study.  This  was  managed,  first 
by  delaying  the  vitamin  D  studies  until  later  in  the  program  as  well  as  we  acquired 
several  founders  from  different  sources  to  try  and  raise  the  colony.  Despite  our  enormous 
efforts  on  breeding  and  creating  the  colony  (approximately  an  18  month  period  of  time 
with  numerous  colony  starts  and  numerous  alterations  in  breeding  conditions)  we  were 
unable  to  establish  a  colony  of  these  animals  sufficient  to  conduct  the  proposed  studies. 
As  a  result  of  this  result,  which  occurred  during  the  last  year  of  funding  (year  3.5  through 
year  4),  we  decided  to  take  a  new  strategy  in  an  effort  to  collect  data  in  the  face  of 
vitamin  deficiency.  These  studies  (as  will  be  described  below)  were  conducted  with  the 
Pi's  own  discretionary  funds,  in  order  to  have  a  more  complete  set  of  data  from  this 
program.  These  studies  are  on-going  and  will  likely  be  completed  within  the  next  6 
months.  A  brief  summary  is  provided  below. 

Vitamin  D  deficiency  study  design 

Using  the  C57BL/6J  mice  on  the  ‘Normal’  Ca2+  from  part  one  of  the  study  as  a  historical 
control  group  (60  mice),  this  study  was  designed  to  take  60  C57BL/6J  mice  and  subject 
them  to  a  modified  diet.  The  diet  was  similar  to  the  control  (Normal)  calcium  diet  from 
the  groups  described  above,  but  was  then  altered  to  have  reduced  vitamin  D.  The  mice 
were  then  exercised  or  not  exercised  according  to  the  group  assigned  at  time  of  delivery. 
The  animals  were  randomly  assigned  to  one  of  the  groups  illustrated  below.  The  exercise 
protocol  was  the  same  as  in  our  previous  studies.  A  total  of  60  animals  were  entered  into 
the  study. 
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The  experimental  phase  of  the  study  is  completed,  and  29  mice  successfully  completed 
the  exercise  regime  along  with  the  30  non-exercised  animals.  All  animals  have  been 
euthanized  and  their  skeletal  structures  extracted. 

Currently  the  left  tibia  from  each  animal  has  been  fixed  in  EtOH  and  being  stored  until 
histological  analysis  can  be  performed..  45  of  the  left  femurs  have  been  scanned  using  the 
GE  MicroCT  system  to  date  and  14  have  already  been  analyzed. 

Within  the  next  6  months  all  the  MicroCT  studies  and  analyses  will  be  complete  as  well 
as  the  mechanical  testing  of  the  femurs  (4  pt.  bending)  and  tibias  (torsion).  We  will 
expect  to  add  these  data  to  our  program  database  as  well  as  publish  the  results. 

C.  Key  Research  Findings 

•  394  mice  that  were  entered  into  the  study  to  evaluate  effects  of  gender,  strain, 
nutrition  (vitamin  C)  and  exercise  on  morphology  and  mechanical  properties  of 
bone  (femurs,  vertebra,  tibia)  have  been  studied. 

•  An  additional  60  animals  have  been  entered  into  a  new  vitamin  D  deficiency  study 
and  have  been  subjected  to  exercise  while  being  supported  by  a  vitamin  D  deficient 
diet.  The  complete  results  will  be  available  in  the  next  few  months.  This  additional 
study  was  conducted  after  the  failure  of  breeding  the  DBP  mice. 
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•  Mouse  strain  had  a  significant  effect  on  geometric,  morphologic  and  mechanical 
properties  in  the  femurs,  vertebra  and  tibias.  C57BL/6J  mice  had  a  significantly 
smaller  cross-sectional  area  and  cortical  thickness  as  compared  to  the  C3H/HeJ 
mice 

•  Gender  had  a  significant  effect  on  all  parameters  (mechanical  and  morphologic) 
tested.  In  general  the  females  were  smaller  and  had  lower  mechanical  properties. 
There  were  no  consistent  differences  in  the  response  to  diet  or  exercise  across 
genders. 

•  Exercise  did  not  have  a  significant  effect  on  any  of  the  femoral  or  vertebral 
mechanical  parameters. 

•  There  were  a  few  significant  effects  of  exercise  on  tibial  morphologic  properties. 
This  observation  was  not  consistent,  suggesting  the  effect  was  not  robust. 

•  In  vertebra,  the  normal  dietary  calcium  group  had  a  significantly  higher  bone 
volume  fraction  as  compared  to  the  low  (p=0.003)  and  high  (p=0.003)  groups.  The 
bone  surface  to  volume  ratio  of  the  normal  group  was  shown  to  be  significantly 
lower  than  the  low  dietary  calcium  group  (p=0.030),  however  no  significant 
differences  were  seen  between  any  of  the  other  groupings. 

•  In  vertebra,  the  normal  dietary  calcium  group  had  a  significantly  higher 
trabecular  thickness  as  compared  to  the  low  group  (p<0.000)  and  marginally 
significant  as  compared  to  the  high  (p=0.051)  group. 

•  Raman  spectroscopy  demonstrated  significant  variations  in  organic  and  inorganic 
constituents  of  the  bone  ECM  across  mouse  strain  and  in  response  to  calcium  intake 
alterations. 

•  There  were  no  histologic  differences  across  any  of  the  groups  or  treatments. 

•  The  application  of  microimaging  and  micromechanical  testing  protocols  was 
found  to  be  outstanding  with  respect  to  characterizing  the  phenotypical  properties 
of  mouse  bone. 


D.  Reportable  Outcomes  (Publications) 

There  were  a  number  of  publications  (abstracts  and  full  length  papers)  that  have  already 
been  completed  as  a  result  of  this  research  program.  As  can  be  gleaned  from  the 
numerous  graphs  of  data  in  the  appendix,  we  are  continuing  to  work  on  the  data  and  will 
have  more  publications  directly  from  the  data  within  the  next  8  months.  Three  categories 
of  outcomes  are  summarized  below: 

a.  Abstracts  and  papers  that  are  a  direct  result  of  the  data  collected 
in  this  study 

b.  Papers  that  were  supported  by  the  program  that  present  the 
methodologies  that  are  new  to  the  field,  but  used  specimens  form 
a  variety  of  studies. 

c.  Papers  that  utilized  the  technologies  developed  in  this  program, 
but  were  not  supported  by  the  program.  Note:  since  the  major 
aim  of  this  study  was  to  develop  techniques  for  evaluating  the 
micro-properties  of  bon  beyond  traditional  density  measures,  the 
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success  of  this  program  has  already  lead  to  other  studies  using 
the  developed  techniques. 


Abstracts  and  papers  directly  supported  and  using  data  generated  from  the  program 

1.  Influence  of  nutrition  and  physical  forces  on  bone  structure/function  properties. 
*Kriegl,  JM,  *Oyserman,  S,  *Roller,  SA,  *Blumenfeld  J,  *Volkman  SK,  *Nashi 
S,  *Hall  JM,  *McCreadie,  BR  and  *Goldstein,  SA.  50th  Annual  Orthopaedic 
Research  Society,  march  2004,  San  Francisco  (abstract) 

2.  Mechanical  loading  effects  structural  morphology  of  male  C57BL/6J  tibia. 
Submitted  for  the  52  Orthopaedic  Research  Society  Meeting,  to  be  held  in 
Chicago  Ill.,  March  2006.  (abstract) 

3.  Morris,  MD,  Crane,  NJ,  Gomez,  LE.  Ignelzi,  MA.  Compatibility  of  staining 
protocols  for  bone  tissue  with  Raman  spectroscopy.  Calcified  Tissue  Int.  74:86- 
94,  2004 

4.  Additional  papers  in  preparation  to  describe  combined  effects  of  exercise,  gender 
and  strain  as  described  in  the  abstracts. 

Papers  supported  by  the  program  utilizing  methods  development 

1.  Carden,  A.,  Morris,  MD.  Rajachar,  RM,  Kohn,  DH.  Ultrastuctural  changes 
accompanying  the  mechanical  deformations:  a  Raman  imaging  study  of  bone 
tissue.  Calcified  Tissue  Int.  72:166-175,  2003. 

Papers  using  the  successfully  developed  techniques,  but  not  supported  by  the  study  funds. 

1.  Chen  T,  Kozloff  KM,  Goldstein  SA,  Morris  MD:  Bone  tissue  ultrastructural 
defects  in  a  mouse  model  for  osteogenesis  imperfecta:  A  Raman  spectroscopy 
study,  Proc.  SPIE,  5321,  85-92,  2004 

2.  Volkman  SK,  Galecki  AJ,  Burke  DT,  Miller  RA,  Goldstein  SA:  Quantitative  trait 
loci  that  modulate  femoral  mechanical  properties  in  a  genetically  heterogeneous 
mouse  population.  J  Bone  Miner  Res,  19(9):1497-1505,  2004 

3.  Kozloff  K,  Carden  A,  Bergwitz  C,  Forlino  A,  Uveges  T,  Morris  M,  Marini  JC, 
Goldstein  SA:  Brittle  IV  mouse  model  for  osteogenesis  imperfecta  IV 
demonstrates  postpubertal  adaptations  to  improve  whole  bone  strength.  J  Bone 
Miner  Res,  19(4):614-622,  2004. 


E.  Conclusions 

The  research  program  has  been  very  successful.  The  primary  aim  was  to  investigate  the 
potential  of  microimaging,  micromechanical  and  Raman  spectroscopic  techniques  in 
characterizing  the  characteristics  of  bone  in  animal  models.  These  studies  demonstrated 
that  the  use  of 

1.  MicroCT 

2.  micromechanical  whole  bone  tests  (femoral  4  pt.  bending,  vertebral 
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compression,  and  tibial  torsion) 

3.  microspecimen  (100  micron  beams)  and  nanoindentation  studies  of  ECM 
tissues  moduli 

4.  Raman  spectroscopy 

5.  histologic  morphology 

can  provide  outstanding  measures  of  bone  phenotypic  properties  in  mouse  models.  The 
studies  also  demonstrated  several  other  key  findings  with  respect  to  the  influence  of 
gender,  nutrition  and  exercise. 

1.  Mouse  strain  is  a  significant  factor  influencing  mechanical  and  morphologic 
properties  of  bone.  The  response  to  nutritional  challenges  is  also  dependent  on 
mouse  strain. 

2.  While  gender  effects  are  substantial  in  absolute  morphologic  and  mechanical 
properties,  the  response  to  nutritional  or  other  influences  does  not  seem  to  be 
gender  dependent.  The  one  possible  variation  in  these  findings  may  be  the 
response  to  exercise.  At  the  protocol  implemented,  there  was  some  suggestion  of 
a  male  response  to  exercise  that  was  not  seen  in  females  (only  in  tibia). 

3.  In  general,  the  exercise  protocol  utilized  in  this  program  (treadmill,  8  weeks)  was 
not  sufficient  to  induce  significant  bone  morphologic  or  mechanical  property 
alterations.  There  was  some  evidence  that  the  tibia  of  males  responded.  This  may 
suggest  that  the  loading  thresh  hold  was  too  load  to  elicit  a  consistent  response, 
and  that  the  tibia  are  more  sensitive  to  the  treadmill  stimulation.  This  is  likely  a 
local  stress  state  phenomena  that  might  be  overcome  by  a  variation  in  the 
boundary  loading  conditions  (magnitude,  rate,  loading  trajectory). 
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All  data  generated  in  study  represented  in  graphical  form 
Histologic  images  from  all  experimental  groups 


INFLUENCES  OF  NUTRITION  AND  PHYSICAL  FORCES  ON  BONE  STRUCTURE/FUNCTION  PROPERTIES 


*Kriegl,  J  M;  *Oyserman,  S  ;  *Roller,  S  A;  *Blumenfeld  J  ;  *Volkman  S  K;  *Nashi  S  ;  *Hall  J  M;  *McCreadie,  B  R;  +*Goldstein,  S  A 
+*Orthopaedic  Research  Laboratories,  University  of  Michigan,  Ann  Arbor,  MI 


INTRODUCTION 

Maintenance  of  bone  integrity  is  dependent  on  a  complex  interaction 
of  metabolic  (hormones,  cytokines,  growth  factors)  and  environmental 
factors  (mechanical  forces,  nutrition  availability)  .  Unfortunately,  the 
specific  relationships  between  these  factors  and  the  biomechanical 
properties  of  bone  tissue  remain  incompletely  quantified.  Understanding 
the  influence  of  nutritional  status  and  mechanical  usage  on  the  properties 
of  bone  may  enhance  our  ability  to  prevent  stress  fractures  associated 
with  intense  training  or  fragility  fractures  accompanying  aging  bone 
loss. 

The  purpose  of  this  study  was  to  investigate  the  interaction  between 
calcium  metabolism  and  exercise  mediated  mechanical  load  on  the 
biomechanical  properties  of  bone. 

METHODS 

C57BL/6J  (n=84  male,  n=90  female)  and  C3H/HeJ  (n=69  male,  n= 

67  female)  mice  were  purchased  from  Jackson  Laboratories  (Bar 
Harbor,  ME)  and  separated  into  groups  based  on  strain,  gender,  diet,  and 
exercise  regime.  A  synthetic  diet  consisting  of  low  (0.02%  Ca),  normal 
(0.95%  Ca),  and  high  (2.0%  Ca)  basal  feed  with  10%  lactose  and  0.67% 
phosphorus  was  purchased  from  Purina  Test  Diet  (Richmond,  IN)2.  An 
exercise  protocol  was  established  and  implemented,  which  consisted  of 
running  the  mice  on  a  custom  built  treadmill  with  a  0°  incline.  The  eight 
week  regime  consisted  of  ramping  the  speed  from  10  to  17  m/min  and 
increasing  the  duration  from  15  to  30  min  over  the  first  4  wks  and 
holding  the  speed  constant  for  an  additional  4  wks3.  The  exercise  regime 
and  the  controlled  calcium  diets  (low,  normal,  high)  were  begun  when 
the  animals  were  12  weeks  old.  The  mice  were  humanely  euthanized  at 
20  wks  of  age,  the  left  femur  and  eighth  caudal  vertebrae  were  dissected 
free  of  soft  tissue  and  frozen  in  LRS.  This  work  was  approved  by  The 
Animal  Care  and  Use  Committee. 

Left  femora  and  Cd8  vertebrae  were  scanned  in  distilled  water  on  a 
GE.EVS  Micro-CT  system.  Each  scan  was  reconstructed  at  a  mesh  size 
of  18pm  x  18pm  x  18pm  and  a  3D  digitized  image  was  generated  for 
each  specimen.  The  femora  and  vertebrae  images  were  rotated  into  a 
standard  orientation  and  thresholded4.  Geometric  analyses  were 
performed  on  a  3mm  mid-diaphysis  segment  in  order  to  obtain  the  cross- 
sectional  area,  cortical  thickness,  and  moment  of  inertia  (Iyy)  for  each 
femur.  Two  standardized  volumes  of  trabecular  bone  were  segmented 
from  the  proximal  and  distal  ends  of  the  vertebral  bodies  for  analysis. 
The  bone  volume  fraction,  bone  surface  to  volume  ratio,  trabecular 
thickness,  trabecular  number,  trabecular  spacing,  and  the  degree  of 
anisotropy  were  acquired  from  the  vertebral  morphologic  analyses. 

The  data  was  analyzed  using  SPSS  statistical  software  (Chicago,  IL). 
Multivariate  general  linear  models  were  used  to  test  for  main  effects  of 
four  factors  (gender,  strain,  diet  and  exercise)  and  for  interaction  effects 
among  the  four  factors.  Post  hoc  tests  were  used  to  compare  the  three 
levels  of  dietary  calcium.  Correlations  were  considered  significant  with 
a  p-value  less  than  0.05. 

RESULTS 

Femora:  Geometric  data  was  collected  from  323  femora.  Mouse 
strain  had  a  significant  effect  on  geometric  parameters.  C57BL/6J  mice 
had  a  significantly  smaller  cross-sectional  area  and  cortical  thickness  as 
compared  to  the  C3H/HeJ  mice.  Exercise  did  not  have  a  significant 
effect  on  any  of  the  geometric  parameters,  as  shown  in  Table  1.  Table  2 
shows  that  differences  in  micro-CT  measures  by  dietary  calcium  level 
are  not  significant. 


Table  1 :  Exercise  Effects  on  Cortical  Parameters 


Exercise 

N 

Cross  Sectional 
Area  (mm2) 

Cortical  Thickness 
(mm) 

Bending  Iyy 
(mm4) 

No 

168 

0.892  (0.20) 

0.256  (0.08) 

0.116(0.03) 

Yes 

155 

0.877  (0.20) 

0.250  (0.08) 

0.114(0.02) 

Data  are  given  as  mean  (STD) 


Table  2:  Dietary  Effects  on  Cortical  Parameters 


Dietary 

Calcium 

N 

Cross  Sectional 
Area  (mm2) 

Cortical  Thickness 
(mm) 

Bending  Iw 
(mm4) 

Low 

105 

0.891  (0.19) 

0.252  (0.08) 

0.119(0.03) 

Normal 

112 

0.884  (0.20) 

0.254  (0.08) 

0.114(0.02) 

High 

106 

0.878  (0.21) 

0.254  (0.08) 

0.112(0.02) 

Data  are  given  as  mean  (STD) 


Vertebrae:  Morphologic  data  was  collected  from  310  vertebrae. 
There  was  no  significant  effect  for  exercise  on  the  trabecular  parameters, 
as  shown  in  Table  3.  The  differences  in  proximal  trabecular  parameters 
for  the  three  dietary  calcium  groups  are  shown  in  Table  4.  The  normal 
dietary  calcium  group  had  a  significantly  higher  bone  volume  fraction  as 
compared  to  the  low  (p=0.003)  and  high  (p=0.003)  groups.  The  bone 
surface  to  volume  ratio  of  the  normal  group  was  shown  to  be 
significantly  lower  than  the  low  dietary  calcium  group  (p=0.030), 
however  no  significant  differences  were  seen  between  any  of  the  other 
groupings.  The  normal  dietary  calcium  group  had  a  significantly  higher 
trabecular  thickness  as  compared  to  the  low  group  (p<0.000)  and 
marginally  significant  as  compared  to  the  high  (p=0.051)  group.  Similar 
results  were  observed  in  the  distal  trabecular  parameters. 


Table  3:  Exercise  Effects  on  Proximal  Trabecular  Parameters 


Bone  Volume 
Fraction  (%) 

Bone  Surface  to 

Trabecular 

Exercise 

N 

Volume  Ratio 
(mm2/mm3) 

Thickness 

(mm) 

No 

161 

40.9  (0.07) 

30.08  (4.57) 

0.074  (0.01) 

Yes 

149 

41.9(0.06) 

29.30  (4.44) 

0.077  (0.01) 

Data  are  given  as  mean  (STD) 

Table  4:  Dietary  Effects  on  Proximal  Trabecular  Parameters 

Dietary 

N 

Bone  Volume 

Bone  Surface  to 
Volume  Ratio 
(mm2/mm3) 

Trabecular 

Thickness 

(mm) 

Calcium 

Fraction  (%) 

Low 

101 

40.6  (0.07)* 

30.29  (4.01)* 

0.073  (0.01)* 

Normal 

109 

42.8  (0.07)* 

29.18(5.26)* 

0.078  (0.02)* 

High 

100 

40.6  (0.06)* 

29.69  (4.07)* 

0.075  (0.01)* 

Data  are  given  as  mean  (STD)  *  Significant,  p  <  0.05 


DISCUSSION 

As  expected,  the  results  indicate  that  geometric  properties  are 
dependent  on  genetic  background.  C57BL/6J  mice  had  a  significantly 
smaller  cross-sectional  area  and  cortical  thickness  as  compared  to  the 
C3H/HeJ  mice,  as  noted  in  previous  studies5. 

Differences  due  to  nutritional  status  were  seen  only  in  trabecular 
bone  of  the  vertebrae,  where  modeling  phenomena  are  expected  to  occur 
more  rapidly  than  in  femoral  cortical  bone  due  to  the  greater  biodynamic 
activity  of  trabecular  bone.  The  loading  regimen  in  this  study  appears  to 
not  have  an  effect  on  the  femoral  geometric  or  vertebral  morphologic 
properties.  The  level  of  exercise  did  not  compensate  for  the  decrease  in 
the  trabeacular  parameters  during  dietary  calcium  alterations.  This  may 
be  caused  by  an  inadequate  load  stimulus.  Preliminary  biomechanical 
testing  data  suggests  a  difference  in  post-yield  behavior  in  exercise 
groups.  This  may  possibly  suggest  an  interaction  at  the  tissue  level, 
which  we  will  continue  to  explore. 
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MECHANICAL  LOADING  EFFECTS  STRUCTURAL  MORPHOLOGY  OF  MALE  C57BL/6J  TIBIAE 


INTRODUCTION 

Although  it  is  known  that  bone  integrity  is  dependent  on  a  complex 
interaction  of  both  metabolic  and  environmental  factors,  the  specific 
relationship  between  these  factors  and  the  biomechanical  properties  of 
bone  tissue  remains  incompletely  quantified.  Previous  work1  has 
illustrated  a  nutritional  dependent  response  of  vertebral  trabecular  and 
femoral  cortical  bone  biomechanical  properties  of  C57BL/6J  and 
C3H/HeJ  mice.  The  purpose  of  this  study  was  to  investigate  the 
influence  of  forced  exercise  on  the  morphologic  properties  of  cortical 
and  trabecular  bone  in  the  tibiae  of  both  male  and  female  C57BL/6J 
mice. 

METHODS 

C57BL/6J  (n=3 1  male,  n=29  female)  mice  were  purchased  from 
Jackson  Laboratories  (Bar  Harbor,  ME).  All  protocols  were  approved 
by  the  Animal  Care  and  Use  Committee.  Starting  at  12  weeks  of  age, 
C57BL/6J  mice  were  fed  a  standard  synthetic  diet  composed  of  0.95% 
Ca  basal  feed  with  10%  lactose  and  0.67%  phosphorus  purchased  from 
Purina  Test  Diet  (Richmond,  IN).  Mice  were  allocated  randomly  to 
either  exercise  or  control  groups.  Exercise  involved  running  eight  weeks 
on  a  custom  built  treadmill  with  a  0°  incline.  Over  the  8-wk  period, 
speed  was  ramped  up  from  10  to  17  m/min  and  duration  was  increased 
from  15  to  30  min  in  the  first  4  wks  and  held  constant  for  the  second  4 
wks.  Exercise  treatment  began  at  12wks  of  age.  At  the  termination  of 
the  exercise  regime,  mice  were  humanely  euthanized  and  the  right  tibiae 
were  removed  and  dissected  free  of  soft  tissue  and  frozen  in  LRS. 

The  right  tibiae  were  scanned  in  distilled  water  on  a  3D  Micro- 
Computed  Tomography  (micro-CT)  system  (GE  Healthcare 
BioSciences).  Scans  were  reconstructed  at  a  mesh  size  of  25  pm  x  25  pm 
x  25pm.  The  tibiae  were  rotated  into  a  standard  orientation  and 
thresholded  based  on  calibrated  histograms  to  separate  bone  and  no-bone 
voxels.  In  each  tibia  a  1.5mm  long  diaphyseal  region  of  interest  was 
defined  just  proximal  to  the  distal  tibiofibular  joint  and  the  bone  was 
evaluated  for  cross-sectional  area,  cortical  thickness,  and  volumetric 
bone  mineral  density  (vBMD).  Bone  volume  fraction,  bone  surface  to 
volume  ratio,  trabecular  thickness,  trabecular  separation,  and  trabecular 
number  were  calculated  from  morphometric  analysis  of  a  2mm  thick 
trabecular  region  just  distal  to  the  proximal  growth  plate.  A  3D  spline 
function  was  utilized  to  define  the  largest  possible  volume  of  trabecular 
bone  in  this  region. 

The  data  was  analyzed  using  SPSS  statistical  software  (Chicago,  IL). 
Multivariate  general  linear  models  were  used  to  test  for  main  effects  of 
two  factors  (gender  and  exercise)  and  for  their  interactions. 

Comparisons  were  considered  significant  with  p-values  of  less  than  0.05. 

RESULTS 

Cortical :  Exercise  induced  a  significant  increase  in  cross-sectional 
area  and  cortical  thickness  (Table  1).  The  results  demonstrated  a 
significant  gender  effect  on  cross-sectional  area  and  cortical  thickness 
(Table  2). 


Table  1 :  Exercise  Effects  on  Cortical  Parameters 


Exercise 

N 

Cross- 

Sectional 

Area  (mm2) 

Cortical 

Thickness 

(mm) 

Volumetric  Bone 
Mineral  Density 
(mg/cc) 

No 

30 

0.554(0.08)* 

0.200(0.02)* 

1441.328(64.8) 

Yes 

30 

0.574(0.07)* 

0.205(0.01)* 

1417.172(38.4) 

Data  are  given  as  mean  (SD)  *  Significant,  p  <  0.05 

Table  2:  Gender  Effects  on  Cortical  Parameters 

Cross- 

Cortical 

Volumetric  Bone 

Gender 

N 

Sectional 

Thickness 

Mineral  Density 

Area  (mm2) 

(mm) 

(mg/cc) 

Male 

31 

0.622(0.05)* 

0.212(0.01)* 

1431.834(57.2) 

Female 

29 

0.502(0.03)* 

0.193(0.01)* 

1426.488(51.7) 

Data  are  given  as  mean  (SD)  *  Significant,  p  <  0.05 


Trabecular  Morphology :  Exercise  was  found  to  have  a  significant 
effect  on  bone  volume  fraction  and  trabecular  number  (Table  3).  A 
gender  effect  was  also  demonstrated  in  the  trabecular  bone  including 
bone  volume  fraction,  trabecular  number  and  vBMD  (Table  4).  Most 
importantly,  evaluation  of  the  interactions  showed  that  only  the  male 


C57BL/6J  mice  showed  a  significant  increase  in  bone  volume  fraction, 
and  trabecular  number  with  exercise  (Table  5).  No  significant  effect 
was  seen  on  trabecular  parameters  in  female  C57BL/6J  mice  (data  not 
shown). 


Table  3:  Exercise  Effects  on  Trabecular  Parameters 


Bone 

Trabecular 

Trabecular 

Volumetric 
Bone  Mineral 
Density 
(mg/cc) 

Exercise 

N 

Volume 
Fraction  (%) 

Thickness 

(mm) 

Number 

(1/mm) 

No 

30 

0.346(0.07)* 

0.608(0.01) 

5.66(0.53)* 

544.496(20.8) 

Yes 

30 

0.375(0.06)* 

0.063(0.01) 

5.90(0.69)* 

547.698(24.1) 

Data  are  given  as  mean  (SD)  *  Significant,  p  <  0.05 


Table  4:  Gender  Effects  on  Trabecular  Parameters 


Bone 

Gende  Volume 

r  Fraction 

(%) 

Trabecula 

r 

Thickness 

(mm) 

Trabecular 

Number 

(1/mm) 

Volumetric 

Bone 

Mineral 

Density 

(mg/cc) 

Male  3  °-400(0-07 

0.063(0.0 

6.276(0.37 

533.519(19.3 

]  \* 

1) 

)* 

)* 

Femal  2  0.319(0.04 

0.060(0.0 

5.258(0.34 

559.541(17.0 

e  9  )* 

1) 

)* 

)* 

Data  are  given  as  mean  (SD)  *  Significant,  p  <  0.05 

Table  5:  Exercise  Effects  of  Trabecular  Parameters 

-  Males 

Bone 

Exercis  Volume 

e  Fraction 

(%) 

Trabecula 

r 

Thickness 

(mm) 

Trabecular 

Number 

(1/mm) 

Volumetric 

Bone 

Mineral 

Density 

(mg/cc) 

3  0.376(0.08 

0.062(0.0 

6.044(0.32 

535.750(21. 

NO  0  )* 

1) 

)* 

1) 

v  3  0.452(0.03 

0.065(0.0 

6.523(0.23 

531.139(17. 

Yes  0  )* 

1) 

)* 

6) 

Data  are  given  as  mean  (SD)  *  Significant,  p  <  0.05 


DISCUSSION 

Following  earlier  work1,2’3,  we  anticipated  that  exercise  would  have  a 
positive  influence  on  the  morphometric  properties  of  the  tibiae.  In 
trabecular  bone,  the  increases  in  both  volume  fraction  and  trabecular 
number  suggest  that  the  effect  of  exercise  is  dominantly  expressed  as  an 
increase  in  the  number  of  trabecular  struts  and  not  an  increase  in  the 
thickness  of  individual  struts.  Importantly,  gender  related  differences 
were  documented.  The  sensitivity  of  males  to  exercise  induced 
increases  in  bone  mass  and  morphology  (in  contrast  to  females)  may 
represent  important  differences  in  hormonal  influences  on  mechano- 
regulation.  These  issues  are  important  areas  for  future  investigation  as 
we  continue  to  search  for  methods  to  enhance  bone  properties  in  the  of 
face  age  or  disease. 
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Abstract.  We  report  the  use  of  Raman  microscopy  to 
image  mouse  calvaria  stained  with  hematoxylin,  eosin 
and  toluidine  blue.  Raman  imaging  of  stained  speci¬ 
mens  allows  for  direct  correlation  of  histological  and 
spectral  information.  A  line-focus  785  nm  laser  imag¬ 
ing  system  with  specialized  near-infrared  (NIR)  micro¬ 
scope  objectives  and  CCD  detector  were  used  to 
collect  approximately  100  x  450  pm  Raman  images. 
Principal  components  analysis,  a  multivariate  analysis 
technique,  was  used  to  determine  whether  the  histolog¬ 
ical  stains  cause  spectral  interference  (band  shifts  or 
intensity  changes)  or  result  in  thermal  damage  to  the 
examined  tissue.  Image  analysis  revealed  factors  for 
tissue  components  and  the  embedding  medium,  glycol 
methacrylate,  only.  Thus,  Raman  imaging  proved  to  be 
compatible  with  histological  stains  such  as  hematoxylin, 
eosin  and  toluidine  blue. 

Key  words:  Raman  spectroscopy  —  Histological  stain¬ 
ing  —  Mineralized  tissue  —  Bone  —  H&E  —  Toluidine 
blue 


Raman  microspectroscopy  and  Raman  imaging  are  in¬ 
creasingly  used  to  provide  chemical  structure  informa¬ 
tion  on  tissue  specimens.  While  spectroscopy  provides 
valuable  information,  biomedical  scientists  usually  wish 
to  link  the  information  to  histological  analysis  [1-4]. 
Depending  upon  the  staining  protocols  employed,  his¬ 
tology  may  provide  morphological  information,  identi¬ 
fication  of  broad  classes  of  molecules  such  as  proteins 
and  nucleic  acids,  or  may  provide  highly  specific  iden¬ 
tification  by  means  of  labeled  antibodies  or  other  mol¬ 
ecules  that  bind  uniquely  to  one  protein.  Histological 
analysis  can  also  identity  regions  of  a  specimen  that 
merit  further  investigation  by  spectroscopy.  Correla¬ 
tive  histological/Raman  imaging  may  provide  infor¬ 
mation  that  is  beyond  what  is  available  by  either 
methodology.  For  such  reasons  it  is  important  to 
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understand  the  compatibility  of  histological  and  spect¬ 
roscopic  protocols  and  to  identify  any  sources  of  in¬ 
compatibility. 

In  this  study  we  investigated  the  compatibility  of 
hematoxylin  and  eosin  (H&E)  staining  and  toluidine 
blue  staining  with  near-infrared  (NIR)  Raman  imaging 
of  fetal  day  18.5  and  postnatal  day  21  mice  calvaria,  or 
flat  bones  that  comprise  the  skull.  Raman  microspectr¬ 
oscopy  and  imaging  provide  molecular  structure  infor¬ 
mation  with  light  microscopy  spatial  resolution  [5]. 
Raman  spectra  are  obtained  by  focusing  a  laser  on  the 
specimen.  While  low  laser  power  is  usually  used,  the 
possibility  of  thermal  damage  is  always  present.  For 
Raman  microscopy  it  is  advantageous  to  use  line-fo¬ 
cused  lasers  that  result  in  low  power  at  any  point  along 
the  line.  This  configuration  simultaneously  reduces  the 
likelihood  of  laser  damage  and  provides  a  quick  survey 
of  the  tissue  along  the  line. 

We  evaluated  H&E  staining  because  it  is  among  the 
most  widely  used  staining  techniques.  Hematoxylin  is  an 
acidophilic  stain  derived  from  the  dye  hematein.  Stain¬ 
ing  with  hematoxylin  results  in  nuclei  and  basophilic 
cytoplasm  colored  blue.  Eosin,  a  basophilic  dye,  is  the 
common  counterstain  for  hematoxylin.  Acidophilic  cy¬ 
toplasm  and  connective  tissue  turn  varying  shades  of 
pink  to  red  with  eosin.  Together,  H&E  clearly  contrast 
the  different  tissue  components.  We  have  also  in¬ 
vestigated  the  effects  of  toluidine  blue,  which  is  often 
used  alone,  but  can  be  substituted  for  hematoxylin. 
Toluidine  blue  also  stains  nuclei  and  basophilic  cyto¬ 
plasm  blue  [6]. 

Raman  spectroscopy  of  tissue  specimens  is  generally 
excited  with  NIR  lasers.  NIR  lasers  are  preferred  to  the 
more  common  visible  wavelength  lasers  because  NIR 
light  does  not  excite  native  fluorescence  from  most 
fluorophores  in  tissues,  although  fluorescence  is  excited 
over  most  of  the  range  of  visible  wavelengths  (400-700 
nm).  Red  stains,  such  as  eosin,  absorb  at  blue-green  to 
ultraviolet  wavelengths  (250-500  nm)  but  not  at  red  or 
near-infrared  wavelengths.  However,  blue  stains,  such  as 
hematoxylin,  absorb  at  red  wavelengths  (650-700  nm) 
and  their  absorbance  spectra  may  extend  into  the 
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near-infrared  (>700  nm).  Thus,  blue  stains  cannot 
be  considered  innocuous  for  Raman  spectroscopy 
until  their  absorption  spectra  have  been  examined  for 
absorbance  at  the  laser  excitation  wavelength  being 
used. 

We  have  previously  shown  that  Raman  image  con¬ 
trast  of  bone  tissue  shows  both  morphology  and  chem¬ 
istry  [7-10].  Using  multivariate  data  analysis  (principal 
component  analysis,  PCA)  even  minor  components  of 
the  specimen  can  be  identified  [1]  and  in  some  cases  can 
be  imaged.  Serial  sections  can  be  used  to  compare 
spectroscopic  results  with  histological  results.  One  sec¬ 
tion  is  examined  by  spectroscopy  and  an  adjacent  sec¬ 
tion  is  subjected  to  histological  analysis.  It  is  accepted 
that  there  are  small  histological  variations  between  serial 
sections,  which  makes  it  somewhat  inaccurate  to  com¬ 
pare  sections.  A  more  objective  approach  of  correlative 
studies  between  spectroscopy  and  histopathology  would 
be  to  examine  spectroscopically  the  same  sections  that 
are  examined  histologically.  In  many  previous  studies, 
spectroscopic  sample  preparation  consists  of  snap 
freezing  tissue  and  sectioning  with  a  cryostat  [11-14]. 
The  examined  tissue  is  unfixed  and  unstained  [14-16]. 
The  sections  subjected  to  histological  analysis  are  fixed, 
typically  in  formalin  which  retains  tissue  rigidity  and 
prevents  bacterial  growth,  then  embedded  and  stained 
with  H&E. 

Correlative  studies  depend  on  the  staining  being 
compatible  with  the  spectroscopic  technique.  There  are 
two  main  concerns  with  staining  and  Raman  micros¬ 
copy.  The  first  is  whether  or  not  the  stain  will  interfere 
spectroscopically.  If  there  is  Raman  signal  from  the 
stain,  it  must  be  distinguishable  from  other  tissue  com¬ 
ponents  by  PCA.  The  second  concern  is  whether  the 
stain  will  absorb  at  the  laser  wavelength  and  generate 
heat  which  in  turn  will  destroy  the  specimen.  Previous 
studies  show  that  basic  fuchsin  does  not  interfere  spec¬ 
troscopically  nor  does  it  cause  thermal  damage  [9].  Here 
we  determine  the  compatibility  of  Raman  microscopy 
and  tissues  stained  with  H&E  and  toluidine  blue. 


Materials  and  Methods 

Experimental  Subjects 

B6CBA  FfiJ  wild-type  mice  (The  Jackson  Laboratory,  Bar 
Harbor,  ME)  were  harvested  at  fetal  day  18.5  and  postnatal 
day  21.  We  used  two  ages  of  mice  to  determine  compatability 
of  staining  and  Raman  imaging  with  delicate  tissue  (fetal  day 
18.5)  and  more  mature  tissue  with  more  mineral  (postnatal  day 
21).  The  calvaria,  flat  bones  from  the  topmost  portion  of  the 
skull,  were  excised  and  rinsed  with  a  phosphate-buffered  saline 
solution  at  pH  7.2.  Specimens  were  then  fixed  in  either  for¬ 
malin  or  ethanol  and  stored  at  4°C.  Prior  to  embedding  with 
glycol  methacrylate  (Technovit  1800  GMA:  Energy  Beam 
Sciences,  Agawam,  MA),  calvaria  were  dissected  and  placed  in 
a  50/50  mixture  of  infiltrating  solution  and  absolute  ethanol 
for  2  hours  and  then  in  infiltrating  solution  overnight.  Infil¬ 


trated  specimens  were  oriented  in  the  polyethylene  mold  tray, 
submerged  with  embedding  solution  and  topped  with  a  plastic 
holder  for  sectioning.  Specimens  were  left  to  cure  overnight 
at  4°C. 

Embedded  specimens  were  sectioned  at  20  pm  with  glass 
knives  in  a  JB-4  microtome  (Energy  Beam  Sciences)  and  then 
mounted  onto  a  glass  microscope  slide  (Probe-On  Plus,  Fisher 
Scientific,  Pittsburgh,  PA)  and  allowed  to  air  dry.  H&E 
staining  procedures  for  GMA-embedded  tissue  were  followed 
for  a  portion  of  the  sections  [17].  Other  sections  were  used 
to  study  modified  staining  protocols  in  which  submersion 
time  in  hematoxylin  solution  was  varied  (30  sec,  1  min,  2  min, 
5  min  and  7  min).  Other  sections  were  stained  with  toluidine 
blue  and  eosin,  again  following  the  GMA-embedded  tissue 
protocol. 


Raman  Analysis 

The  Raman  imaging  system  has  been  described  previously  [7], 
Briefly,  about  160  mW  of  line-focused  laser  light  from  a  785 
run  diode  laser  (Invictus,  Kaiser  Optical  Systems,  Inc.,  Ann 
Arbor,  MI)  is  delivered  to  the  sample  via  epi-illumination  with 
a  microscope  (BH-2;  Olympus,  Inc.,  Melville,  NY)  equipped 
with  a  NIR-optimized  20  x  0.75  NA  objective  (Fluar  Series, 
Carl  Zeiss,  Inc.,  Thornwood,  NY).  Raman  scatter  collected 
through  the  objective  is  focused  into  an  axial-transmissive 
spectrograph  (HoloSpec  f/1.8i,  Kaiser  Optical  Systems,  Inc., 
Ann  Arbor,  MI)  operated  at  3-4  cm-1  resolution  (25  pm  en¬ 
trance  slit).  The  dispersed  signal  is  collected  by  a  1024  x  128 
pixel  charge  coupled  device  camera  (DU401-BR-DD;  Andor 
Technology,  South  Windsor,  CT)  operated  at  -75°C.  To  ob¬ 
tain  an  image,  a  stepper-motor  controlled  x-y  translation  mi¬ 
croscope  stage  (NEAT,  Salem,  NH)  moves  the  sample  in  1.5 
pm  steps  for  the  desired  number  of  frames,  providing  3  pm 
spatial  resolution,  with  3  min  acquisitions  per  frame.  In  most 
cases  150  frames  were  collected,  to  yield  a  126  x  150  pixel 
image  containing  information  from  all  18,900  Raman  spectra. 
Brightfield  images  were  collected  by  a  color  video  camera 
(TMC-7;  Pulnix,  Sunnyvale,  CA)  mounted  on  the  microscope 
frame.  All  absorbance  spectra  were  collected  with  a  Shimadzu 
2101  Scanning  UV-Visible  Spectrophotometer.  All  scans  were 
performed  from  550-800  nm,  with  0.1  nm  slit  width,  0.5  nm 
sampling  intervals  and  a  fast  scan  speed. 


Factor  Analysis 

Factor  analysis  [19],  which  has  been  discussed  previously  [10, 
18],  was  performed  with  Matlab  5.3  (MathWorks,  Inc.,  Na¬ 
tick,  MA),  using  both  manufacturer-supplied  routines  as  well 
as  in-house  developed  routines.  The  entire  set  of  18,900  Ra¬ 
man  spectra  that  comprise  an  image  is  treated  as  a  matrix  of 
18,900  rows  (spectra)  and  1024  columns  (Raman  shifts).  The 
principal  components,  which  are  the  eigenvectors  of  this  ma¬ 
trix,  are  calculated  and  then  sorted  by  their  magnitudes,  called 
eigenvalues.  The  standard  procedure  for  calculating  eigenvec¬ 
tors  is  called  "singular  value  decomposition”  and  is  included  in 
Matlab.  The  eigenvectors  are  a  mathematically  orthogonal 
coordinate  system  (i.e.,  extensions  of  the  familiar  3-dimen- 
sional  Cartesian  coordinate  system)  from  which  the  original 
data  set  can  be  calculated.  The  number  of  eigenvectors  found 
is  equal  to  the  number  of  Raman  spectra.  However,  almost  all 
of  the  non-random  variation  in  the  image  data  set  can  be  de¬ 
scribed  by  the  eigenvectors  with  the  largest  magnitudes.  Usu¬ 
ally  between  four  to  eight  eigenvectors  are  sufficient.  The 
remaining  eigenvectors  describe  the  random  noise  in  the  data 
set  and  are  not  used  in  the  data  analysis.  The  criterion  for 
eigenvector  retention  is  that  addition  of  the  next  eigenvector 
describes  no  more  than  0.1%  of  the  variation  in  the  data  set. 
The  intensities  of  these  eigenvectors  are  called  scores. 

Principal  components  are  mathematically  orthogonal  and  so 
may  contain  both  positive  and  negative  bands.  They  are  sums 
and  differences  of  the  underlying  Raman  spectra.  To  convert 
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Fig.  1.  Absorbance  spectra  of  0.01%  hematoxylin  solution  at  pH  5.2  and  8,  0.01%  eosin  solution  and  0.01%  toluidine  blue 
solution.  The  arrow  points  to  absorbances  at  the  wavelength  of  excitation,  approximately  785  nm. 


the  retained  principal  components  to  Raman  spectra,  weighted 
sums  and  differences  are  calculated,  with  the  weighting  subject 
to  two  constraints.  The  first  constraint  is  that  at  each  Raman 
shift  the  calculated  weighted  sums  and  differences  be  non¬ 
negative,  as  are  the  Raman  spectra.  A  second  constraint  is 
that  at  each  point  the  reconstructed  Raman  spectra  must  add 
up  to  the  original  intensity.  These  constraints  are  sufficient  to 
guarantee  that  the  reconstructed  spectra,  called  factors,  will  be 
directly  proportional  to  the  real  Raman  spectra.  The  rela¬ 
tive  intensity  (score)  images  will  be  positive  and  will  be  maps  of 
the  spatial  distribution  and  relative  amounts  of  indicating 
where  and  how  much  of  the  components  are  present  in  the 
specimen. 

Results 

Absorption  of  Histological  Stains 

In  order  to  determine  if  eosin,  hematoxylin  and  tolui¬ 
dine  blue  would  cause  sufficient  heating  upon  laser  ex¬ 
posure  to  cause  tissue  damage,  the  absorption  at  the 
laser  wavelength  was  determined.  Figure  1  shows  the 
absorption  spectra  of  eosin,  hematoxylin  and  toluidine 
blue  from  550  nm  to  800  nm.  Although  the  spectra  were 
obtained  in  solutions  of  the  stains  for  experimental 
convenience,  the  spectra  of  the  solid  phase  stains  are 
similar.  Eosin  shows  no  detectable  absorbance  at  785 
nm,  as  would  be  expected  for  a  red  stain.  Toluidine  blue 
absorbs  weakly  at  785  nm  (0.04  from  a  0.01%  solution). 
Hematoxylin  absorbs  more  strongly  at  785  nm  (0.4  at 
pH  5.2  and  0.3  at  pH  8,  from  a  0.01%  solution). 


Light  Microscopy  of  Stained  Specimens 

Brightfield  images  of  mouse  calvaria  embedded  in  GMA 
and  stained  according  to  the  protocols  are  shown  (Fig. 
2).  Raman  images  of  regions  of  these  specimens  are 
presented  and  discussed  below.  The  image  of  a  control 
section  of  GMA  (no  tissue),  part  of  which  is  stained  with 
eosin,  is  shown  in  Figure  2a.  The  image  of  a  second 
control,  a  section  of  GMA-embedded  and  unstained 
calvarial  tissue  from  a  21 -day-old  wild-type  mouse  is 
shown  in  Figure  2b.  Figure  2c  shows  the  same  section 
after  staining  with  eosin.  Note,  the  holes  in  2b  and  2c  are 
a  result  of  the  embedding  process,  not  tissue  damage  by 
the  laser.  Consistent  with  the  data  presented  in  Figure  1, 
no  thermal  damage  was  detectable  after  Raman  imaging 
was  performed  on  these  sections.  Figure  2d  shows  a 
section  of  fetal  day  18.5  calvarial  tissue  stained  with 
eosin  and  toluidine  blue.  The  presence  of  toluidine  blue 
did  not  result  in  any  detectable  thermal  damage,  even 
though  it  is  weakly  absorbing  at  785  nm. 

Several  H&E  protocols  resulted  in  enough  bound 
hematoxylin  to  cause  thermal  damage  under  785  nm  la¬ 
ser  illumination.  Figure  2e  shows  a  section  of  calvarial 
tissue  from  a  21 -day-old  wild-type  mouse  stained  by  an 
H&E  protocol  that  consisted  of  requiring  15  min  im¬ 
mersion  in  hematoxylin.  Note  the  laser  damage  to  the 
upper  right  hand  corner  of  the  specimen.  Visible  damage, 
burning  of  the  tissue,  often  accompanied  by  noticeable 
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Fig.  2.  Brightfield  images  of  sections  studied  with  Raman 
imaging.  Sections  were  of  the  following:  (a)  GMA  partially 
stained  with  eosin,  (b)  postnatal  day  21  mouse  calvaria  em¬ 
bedded  in  GMA  and  unstained,  (c)  postnatal  day  21  mouse 
calvaria  embedded  in  GMA  and  stained  with  eosin,  (d)  fetal 

odor,  occurred  as  soon  as  the  laser  was  incident  on 
the  specimen.  A  video  image  was  obtained  less  than  1 
second  later.  Not  surprisingly  then,  protocols  requiring 
longer  immersion  in  hematoxylin  also  resulted  in 
damage. 

Hematoxylin  immersion  was  reduced  in  a  series 
of  steps  to  find  conditions  that  resulted  in  good  contrast 
in  reflectance  imaging  with  no  visible  or  spectroscopically 
detectable  thermal  damage.  Fetal  tissue  was  used  for  test¬ 
ing  these  protocols  because  it  is  more  easily  damaged  than 
mature  tissue.  Immersion  times  of  7,  5,  2,  1  and  0.5  min 
were  examined.  Sections  immersed  in  hematoxylin  for  2, 
5,  and  7  min  burned  immediately.  Specimens  immersed  in 
hematoxylin  for  1  min  sustained  1.5  min  of  laser  exposure 
before  visible  damage  was  observed.  Immersion  in 
hematoxylin  solution  for  0.5  min  resulted  in  no  thermal 
damage,  even  after  more  than  5  min  of  exposure  to  the 
laser.  There  was  little  difference  in  contrast  obtained  after 
the  0.5-minute  immersion  and  the  longer  H&E  staining 
protocols,  as  shown  by  comparison  (Fig.  2e,f). 


day  18.5  mouse  calvaria  embedded  in  GMA  and  stained  with 
toluidine  blue  and  eosin,  (e)  fetal  day  18.5  mouse  calvaria 
embedded  in  GMA  and  stained  by  standard  H&E  protocol,  (f) 
fetal  day  18.5  mouse  calvaria  embedded  in  GMA  and  stained 
by  modified  H&E  protocol. 

Raman  Imaging  of  GMA  Stained  with  Alcoholic  Eosin 

Figure  2a  shows  a  reflectance  brightfield  image  of  a 
GMA  section  (no  tissue)  that  was  used  to  test  for  spec¬ 
troscopic  interference  front  eosin  (only  part  of  the  section 
was  stained  with  eosin).  Figure  3a  shows  the  only  Raman 
factor  recovered  front  the  Raman  imaging  experiment 
and  3b  shows  the  corresponding  score  image. 

In  Figure  3a  only  the  spectral  bands  between  420 
cm-1  and  1045  cm-1  are  shown  in  order  to  bring  out 
detail  that  emphasizes  the  similarities  between  corre¬ 
sponding  factors  obtained  here  and  in  the  staining 
protocol  experiments  described  below.  Spectra  were 
obtained  over  the  420-1800  cm"1  range  and  contain 
additional  Raman  bands  that  are  not  shown.  For  con¬ 
sistency,  the  Raman  images  throughout  the  paper  have 
similar  dimensions  to  the  image  of  Figure  3b,  approxi¬ 
mately  100  pm  x  450  pm. 

If  the  eosin  stain  produces  a  Raman  spectrum  that  is 
detectable,  there  will  be  eosin  bands  in  the  Raman 
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Fig.  3.  Raman  factor  for  a  section  of  GMA  partially  stained  with  eosin  (a)  and  the  corresponding  score  image  (b).  Note  only 
GMA  is  detectable.  The  section  examined  was  20  pm  thick,  mounted  on  a  glass  microscope  slide.  Laser  power  incident  to  the 
sample  was  185  mW. 


spectra  from  this  part  of  the  section  but  not  from  the 
unstained  part.  The  extra  bands  should  be  resolvable  by 
factor  analysis  even  if  they  are  not  readily  visible  upon 
inspection.  However,  as  indicated  in  Figure  3,  the  eosin 
is  not  visible  spectroscopically.  Only  one  factor  con¬ 
taining  Raman  bands  is  obtained.  This  factor  is  repre¬ 
sentative  of  the  typical  GMA  spectrum,  with  signature 
bands  at  600  cm-1  and  968  cm-1.  The  score  (intensity)  is 
uniform  across  the  section,  with  no  difference  between 
the  stained  and  unstained  regions.  Only  striations  from 
the  microtome  glass  knife  and  one  surface  irregularity 
mar  the  otherwise  flat  surface.  We  conclude  that  if  eosin 
spectra  are  present,  they  must  contribute  much  less  than 
1%  to  the  total  observed  Raman  intensity. 

Raman  Imaging  of  Calvaria  Embedded  in  GMA  and  Unstained 

To  determine  whether  spectral  interference  would  result 
from  the  application  of  histological  stains,  we  first  re¬ 
solved  the  factors  of  the  GMA  embedded  tissue  without 
staining.  Figure  2b  is  the  reflectance  brightfield  image  of 
the  region  of  a  section  of  calvaria  from  a  21 -day  post¬ 
natal  mouse.  The  specimen  here,  calvaria  from  a  21 -day 
postnatal  mouse,  is  embedded  in  GMA  and  unstained. 
It  serves  as  a  control  for  comparison  with  images  from 
stained  specimens.  Three  non-background  Raman 
spectrum  factors  are  observed.  The  factors  and  the 
corresponding  score  images  are  shown  in  Figure  4. 

The  GMA  factor  (Fig.  4a)  is  identical  to  the  GMA 
factor  shown  in  Figure  3a.  Bone  matrix  (Fig.  4c)  and 
bone  mineral  (Fig.  4e)  factors  are  also  obtained.  These 
are  identical  to  factors  that  we  have  previously  reported 
for  similar  mouse  calvaria  specimens  [7].  The  corre¬ 
sponding  score  images  are  Figures  4b  (GMA),  4d  (bone 
matrix)  and  4f  (bone  mineral).  There  is  some  overlap 
between  the  GMA  image  and  the  bone  mineral  and 
matrix  images,  but  it  is  not  great.  Though  there  may  be 
some  infiltration  into  the  tissue,  the  GMA  shows 
greatest  intensity  where  bone  is  not  present  in  the  sec¬ 


tion.  It  is  clear  that  factor  analysis  adequately  separates 
the  components  of  the  specimen,  which  is  composed  of 
GMA  (Fig.  4a),  bone  matrix  (Fig.  4c)  and  bone  mineral 
(Fig.  4e). 

Raman  Imaging  of  Calvaria  Embedded  in  GMA  and  Stained 

In  Figure  5  Raman  score  images  from  a  selected  region 
of  eosin-stained  tissue  specimen  of  Figure  2c  are  shown. 
Three  non-background  factors  that  contain  Raman 
spectral  information  for  the  specimen  were  found.  The 
remaining  factors  contained  only  background  lumines¬ 
cence  and  were  not  included.  The  three  Raman  fac¬ 
tors  correspond  to  GMA  (Fig.  5a),  bone  mineral  (Fig. 
5c)  and  bone  matrix  (Fig.  5e).  As  found  for  the  un¬ 
stained  GMA-embedded  calvaria,  the  GMA  is  found 
(Fig.  5b)  mostly  where  the  bone  mineral  (Fig.  5d)  is 
absent,  although  it  does  fill  small  depressions  in 
the  tissue.  In  comparing  the  bone  mineral  and  matrix 
(Fig.  5f)  score  images  there  is  a  relatively  uniform 
distribution  of  mineral  on  the  matrix,  as  would  be 
expected  of  a  small  region  of  mature  tissue  such  as  this 
21-day-old  mouse.  There  was  no  factor  found  for  the 
eosin.  The  stain  is  probably  present  in  a  sufficiently 
enough  thin  layer  that  it  cannot  be  detected  spectrosc¬ 
opically. 

A  fetal  day  18.5  mouse  calvarial  section  stained  with 
toluidine  blue  and  eosin  is  shown  in  Figure  6.  As  dis¬ 
cussed  above,  the  tissue  was  chosen  because  it  is  more 
easily  damaged  than  mature  tissue.  There  are  two  non¬ 
background  factors,  GMA  (Fig.  6a)  and  bone  (Fig.  6c). 
The  corresponding  score  images  are  shown  as  Fig.  6b 
and  6d.  Although  mineralization  is  expected  and  present 
in  this  tissue,  it  is  not  as  extensive  as  in  the  postnatal 
day  21  tissue.  The  bone  factor  appears  at  lower  signal/ 
noise  ratio  than  in  the  more  mature  tissue  specimens 
used  in  the  Figs.  4  and  5.  No  Raman  factors  are  found 
for  toluidine  blue  because  the  stain  is  present  in  a  thin 
layer. 


M.  D.  Morris  et  al.:  Compatibility  of  Staining  Protocols  for  Bone  Tissue  with  Raman  Imaging 


91 


High 


LOW 


Fig.  4.  Raman  factors  and  score  images  for  postnatal  day  21 
mouse  calvaria  embedded  in  GMA  and  unstained.  Shown  is 
the  Raman  factor  for  GMA  (a),  the  Raman  factor  for  bone 
matrix  (c)  and  the  Raman  factor  for  bone  mineral  (e).  The- 
corresponding  score  images  are  shown  in  b,  d  and  f.  Note,  in 


the  score  image,  the  GMA  is  least  intense  where  the  bone 
mineral  is  the  most  intense.  The  section  examined  was  20  pm 
thick  and  mounted  on  a  glass  microscope  slide.  Laser  power 
incident  to  the  sample  was  153  mW. 


Fig.  5.  Raman  factors  and  score  images  for  postnatal  day  21 
mouse  calvaria  embedded  in  GMA  and  stained  with  eosin. 
Shown  is  the  Raman  factor  for  GMA  (a),  the  Raman  factor 
for  bone  mineral  (c)  and  the  Raman  factor  for  bone  matrix  (e). 


The  corresponding  score  images  are  shown  in  b,  d  and  f.  The 
section  examined  was  20  pm  thick  and  mounted  on  a  glass 
microscope  slide.  Laser  power  incident  to  the  sample  was  157 
mW. 
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Fig.  6.  Raman  factors  and  score  images  for  fetal  day  18.5  score  images  are  shown  in  (b)  and  (d).  The  section  examined 

mouse  calvaria  embedded  in  GMA  and  stained  with  toluidine  was  20  pm  thick  and  mounted  on  a  glass  microscope  slide, 

blue  and  eosin.  Shown  is  the  Raman  factor  for  GMA  (a)  and  Laser  power  incident  to  the  sample  was  185  mW. 

the  Raman  factor  for  bone  mineral  (c).  The-corresponding 


Raman  Shift  (cnr') 

Fig.  7.  Raman  factors  and  score  images  for  fetal  day  18.5  mouse  calvaria  embedded  in  GMA  and  stained  with  H&E.  The  Raman 
factor  for  bone  matrix  (a)  and  the  Raman  factor  for  GMA  (c).  The  corresponding  score  images  are  shown  in  (b)  and  (d).  The 
section  examined  was  20  pm  thick  and  mounted  on  a  glass  microscope  slide.  Laser  power  incident  to  the  sample  was  185  mW. 

Figure  7  shows  data  from  the  suture  of  fetal  day  18.5  damaged  areas.  However,  spectral  images  were  easily 

tissue  stained  with  H&E.  This  is  the  most  delicate  and  obtainable  from  tissue  stained  by  immersion  for  30 

spectroscopically  challenging  tissue  used  in  this  study  seconds  in  hematoxylin.  No  hematoxylin  factor  was 

because  collagen  is  a  significantly  weaker  Raman  scat-  found  and  no  hematoxylin  bands  are  present  in  the 

terer  than  the  mineral  in  bone.  The  factor  shown  in  factors. 

Figure  7a  is  the  spectrum  of  bone  matrix.  The  score 

image  is  shown  in  7b.  The  other  non-background  factor 

is  GMA  (Fig.  7c);  its  score  image  is  seen  in  Fig.  7d.  Discussion 

Calvarial  sections  stained  according  to  the  1,  2,  5,  7  and 

15  min  H&E  protocols  suffered  immediate  thermal  Because  light  absorption  can  result  in  sufficient  local 

damage.  No  Raman  spectra  were  obtainable  from  heating  to  damage  tissues,  staining  protocols  were  tested 
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for  thermal  compatibility  with  Raman  imaging.  From 
the  absorption  spectra,  one  would  expect  that  eosin 
staining  would  not  result  in  thermal  damage,  even  to 
delicate  tissue,  but  that  damage  could  result  from  ab¬ 
sorption  by  either  of  the  blue  stains,  hematoxylin  or 
toluidine  blue.  Stained  sections  should  always  be  in¬ 
spected  for  thermal  damage  after  exposure  to  the  laser 
light.  Additionally,  even  if  there  is  no  visual  evidence  of 
thermal  damage,  Raman  spectra  from  control  unstained 
sections  and  adjacent  stained  sections  should  be  com¬ 
pared  for  subtle  differences.  In  particular,  Raman  shift 
frequencies  and  band  intensities  change  when  a  protein 
is  denatured  thermally  or  chemically.  Any  changes  must 
be  presumed  to  be  evidence  of  thermal  damage  to  the 
tissue,  unless  further  study  demonstrates  otherwise. 
Similarly,  chemical  alterations  as  a  result  of  staining 
could  result  in  Raman  factors  different  from  the  Raman 
factors  of  unaffected  tissue.  However,  the  stain  is  usually 
a  superficial  layer  that  affects  a  small  fraction  of  the 
specimen.  All  band  frequency  and  intensity  changes  in 
this  study  are  attributed  to  sample-to-sample  variations. 

General  conclusions  can  be  drawn  about  the  com¬ 
patibility  of  Raman  microspectroscopy  or  imaging  and 
histological  staining.  There  are  no  compatibility  prob¬ 
lems  if  a  stain  does  not  absorb  laser  radiation.  If  near- 
infrared  lasers  are  used,  one  would  not  expect  thermal 
damage  resulting  from  red  stains  such  as  eosin  or  basic 
fuchsin.  Blue  stains  cannot  be  assumed  to  be  innocuous. 
With  the  785  nrn  laser  used  in  these  experiments,  tolui¬ 
dine  blue  does  not  cause  any  thermal  damage  because  it 
is  only  weakly  absorbing  at  this  wavelength.  Hematox¬ 
ylin  absorbs  785  nm  more  strongly  but  by  reducing  the 
hematoxylin  immersion  time  to  0.5  min,  thermal  dam¬ 
age  was  eliminated  with  only  slight  loss  in  contrast.  It 
may  be  necessary  to  modify  standard  staining  protocols 
when  using  a  laser  line  that  lies  within  the  absorption 
spectrum  of  a  given  stain. 

A  second  concern  when  using  Raman  imaging  is 
whether  or  not  the  stain(s)  will  spectroscopically  inter¬ 
fere  with  the  Raman  signal  of  tissue  components.  It  is 
clear  that  stains,  i.e.,  eosin,  hematoxylin  and  toluidine 
blue,  do  not  interfere  with  Raman  imaging  as  evidenced 
by  lack  of  Raman  factors  in  Figures  3,  4,  5,  6,  7.  Im¬ 
mersion  in  these  stains  does  not  interfere  with  resolution 
of  tissue  components  by  factor  analysis.  Comparisons 
with  unstained  specimens  revealed  that  the  factors  rep¬ 
resenting  bone  mineral  and  bone  matrix  are  unchanged. 
Close  examination  and  comparison  of  Figures  4  (un¬ 
stained  mouse  calvaria  embedded  in  GMA)  and  5 
(mouse  calvaria  embedded  in  GMA  and  stained  with 
eosin)  reveal  no  spectral  differences. 

NIR  Raman  spectroscopy  is  favored  for  tissue  spec¬ 
imens  because  tissue  fluorescence  is  only  weakly  excited. 
Line  focusing  reduces  local  power  density  more  than 
lOOx  compared  to  point  focusing  of  the  same  laser  beam 
[20,  21].  This  advantage  is  obtained  even  if  the  goal  of  an 


experiment  is  not  acquisition  of  a  complete  two-di¬ 
mensional  spectroscopic  image  of  a  specimen.  While  line 
focusing  reduces  the  likelihood  of  damage  to  delicate 
specimens,  it  cannot  eliminate  the  problem  completely. 
Therefore,  it  is  prudent  to  test  for  thermal  damage  be¬ 
fore  attempting  Raman  microspectroscopy  or  imaging 
of  tissue  that  has  been  subjected  to  blue  staining. 
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Abstract.  Raman  spectroscopy  and  imaging  are  known 
to  be  valuable  tools  for  the  analysis  of  bone,  the  deter¬ 
mination  of  protein  secondary  structure,  and  the  study 
of  the  composition  of  crystalline  materials.  We  have 
utilized  all  of  these  attributes  to  examine  how  mechan¬ 
ical  loading  and  the  resulting  deformation  alfects  bone 
ultrastructure,  addressing  the  hypothesis  that  bone 
spectra  are  altered,  in  both  the  organic  and  inorganic 
regions,  in  response  to  mechanical  loading/deformation. 
Using  a  cylindrical  indenter,  we  have  permanently  de¬ 
formed  bovine  cortical  bone  specimens  and  investigated 
the  ultrastructure  in  and  around  the  deformed  areas 
using  hyperspectral  Raman  imaging  coupled  with  mul¬ 
tivariate  analysis  techniques.  Indent  morphology  was 
further  examined  using  scanning  electron  microscopy. 
Raman  images  taken  at  the  edge  of  the  indents  show 
increases  in  the  low-frequency  component  of  the  amide 
III  band  and  high-frequency  component  of  the  amide  I 
band.  These  changes  are  indicative  of  the  rupture  of 
collagen  crosslinks  due  to  shear  forces  exerted  by  the 
indenter  passing  through  the  bone.  However,  within  the 
indent  itself  no  evidence  was  seen  of  crosslink  rupture, 
indicating  that  only  compression  of  the  organic  matrix 
takes  place  in  this  region.  We  also  present  evidence  of 
what  is  possibly  a  pressure-induced  structural  transfor¬ 
mation  occurring  in  the  bone  mineral  within  the  indents, 
as  indicated  by  the  appearance  of  additional  mineral 
factors  in  Raman  image  data  from  indented  areas.  These 
results  give  new  insight  into  the  mechanisms  and  causes 
of  bone  failure  at  the  ultrastructural  level. 

Key  words:  Raman  —  Bone  —  Biomechanics  —  Spec¬ 
troscopic  imaging  —  Fracture 


Like  most  tissues,  bone  is  hierarchical  in  structure  [1], 
Any  study  of  bone  properties  must  investigate  the  tissue 
at  several  levels  of  organization  in  order  to  gain  a 
complete  understanding  of  the  influence  of  structure  and 
composition  on  these  properties.  This  is  especially  true 
when  studying  mechanical  properties.  Many  studies 
have  looked  at  the  effect  of  applied  load  on  bone,  but  on 
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a  macroscopic  level  [1],  A  few  have  looked  at  a  lower, 
microstructural  level;  for  example,  studies  of  micro¬ 
cracking  [2^1],  tensile  properties  of  single  osteons  [5], 
and  nanoindentation  studies  of  individual  lamellae  [6, 
7],  However,  studies  at  an  ultrastructural  or  chemical 
level  have  been  impossible  using  the  most  common 
techniques.  As  a  result  of  these  limitations,  it  is  not 
known  definitively  how  individual  collagen  fibrils  or 
mineral  crystallites  react  to  mechanical  load  in  either  the 
elastic  or  plastic  deformation  regimes.  There  are  many 
hypotheses — rupture  of  crosslinks  in  the  collagen  fibrils, 
mineral  phase  transformations  or  rearrangements,  sep¬ 
aration  of  fibrils,  and/or  relative  motion  (“gliding”)  of 
fibrils  [8,  9].  However,  new  technical  approaches  are 
needed  to  address  these  hypotheses  and  to  study  the 
ultrastructural  mechanisms  guiding  the  response  of  bone 
to  mechanical  loading. 

In  recent  years,  vibrational  spectroscopy  has  enjoyed 
increasing  success  as  a  technique  for  studying  bone 
properties  [10].  Infrared  and  Raman  spectroscopies  have 
been  used  to  study  normal  and  osteoporotic  tissue,  as 
well  as  to  investigate  compositional  changes  associated 
with  various  bone  diseases  [10-13].  Raman  spectroscopy 
has  proven  particularly  useful  because  of  its  ability  to  be 
used  on  intact  specimens  as  well  as  its  lack  of  strong 
interferences  from  water.  Raman  spectra  can  be  col¬ 
lected  with  micron-scale  spatial  resolution  and  give  si¬ 
multaneous  information  about  the  mineral  and  organic 
matrix  components  of  bone.  Raman  spectra  can  also  be 
collected  over  a  two-dimensional  area  to  give  spatially 
relevant  image  information.  With  the  use  of  multivariate 
statistics,  small  local  variations  in  bone  composition  can 
be  resolved  [14-16].  Information  from  the  mineral 
component  and  the  organic  matrix  is  obtained  simulta¬ 
neously,  creating  a  complete  picture  of  bone  composi¬ 
tion  in  the  area  surveyed.  The  phosphate  V|  band  at 
approximately  960  cnf 1  and  the  bands  associated  with 
collagen  (amide  III  at  ~  1250  cnf1,  CH2  wag  at  ~  1450 
cm”1,  and  amide  I  at  ~1665  cm”1)  are  of  particular  in¬ 
terest  for  bone  compositional  studies  [10]. 
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Fig.  1.  Raman  imaging  apparatus. 


Raman  spectroscopy  is  a  well-established  technique 
for  looking  at  changes  brought  on  by  mechanical  loading 
in  such  materials  as  polymers  [17,  18]  and  protein  films 

[19] ,  as  well  as  changes  in  protein  structure  in  solution 

[20] .  Generally,  structural  changes  induced  by  an  applied 
load  cause  shifts  in  bands  in  the  Raman  spectrum  [17- 
19].  The  amide  I  band,  at  about  1665  cm-1,  and  the  amide 
III  band,  at  about  1250  cm-1,  are  known  to  be  especially 
good  indicators  of  protein  conformation  because  of  the 
amide  moiety’s  role  in  crosslinking  and  bonding  [21]. 
Spectra  of  crystalline  materials  such  as  silicon  or  calcium 
hydroxyapatite  also  show  shifts  in  spectral  bands  with 
applied  load,  presumably  as  a  result  of  pressure-induced 
rearrangement  of  the  crystalline  lattice  [22-25]. 

Recent  experiments  in  this  lab  using  Raman  spec¬ 
troscopy  have  started  to  provide  insight  into  damage 
mechanisms  at  the  ultrastructural  level.  Increases  in  the 
number  of  mineral  factors  (phosphate  Vi  spectral  com¬ 
ponents)  have  been  observed  at  the  leading  edge  of  fa¬ 
tigue-induced  microcracks  in  cortical  bone  [14].  It  was 
not  clear  from  these  studies  whether  the  spectroscopic 
changes  were  a  cause  or  an  effect  of  the  damage.  Sub¬ 
sequent  studies,  which  analyzed  spectral  changes  in  real¬ 
time  during  compressive  and  tensile  loading  of  bone, 
provided  support  for  the  hypothesis  that  spectral  shifts 
result  from  mechanical  damage,  and  demonstrated  that 
changes  in  the  organic  matrix  also  result  from  me¬ 
chanical  loading  [26].  In  the  present  study,  we  address 
the  hypothesis  that  bone  spectra  are  altered,  in  both  the 
organic  and  inorganic  regions,  in  response  to  mechani¬ 
cal  loading/deformation. 

We  report  here  the  results  of  experiments  on  bovine 
tissue  that  was  mechanically  loaded  and  locally  deformed 


via  cylindrical  indentation.  Our  results  show  that  most  of 
the  damage  to  the  organic  matrix  occurs  at  the  edges  of 
indents  rather  than  in  the  center.  This  damage  is  in  the 
form  of  ruptured  cross-links,  as  indicated  by  spectro¬ 
scopic  shifts  in  amide  I  and  amide  III  bands.  Addition¬ 
ally,  there  is  evidence  that  mineral  structural  transitions 
are  occurring  within  the  indented  areas,  as  evidenced  by 
the  increased  numbers  of  mineral  factors  obtained  from 
the  analysis  of  Raman  image  data  from  these  regions 
relative  to  control  areas.  These  results  support  the  hy¬ 
pothesis  that  additional  mineral  factors  observed  in  areas 
of  damage  are  an  effect  rather  than  a  cause  of  the  damage. 

Materials  and  Methods 

Specimens 

Bovine  bone  from  skeletally  mature  animals  was  obtained 
from  a  local  abattoir.  A  section  from  the  lateral  anterior 
portion  of  the  central  annulus  of  a  femur  was  milled  into  a  flat 
rectangular  beam  approximately  3  cm  x  1  cm  x  15  cm  and 
stored  frozen  at  -64°C.  For  these  experiments,  the  beam  was 
cut  into  4  pieces,  1  of  which  was  reserved  for  a  separate  ex¬ 
periment  and  was  not  used  further  in  these  studies.  The  re¬ 
maining  specimens  were  polished  on  one  side  using  a  graded 
series  of  grits,  followed  by  a  0.3  pm  alumina  slurry.  The 
specimens  were  sonicated  briefly  to  remove  any  particles  and 
were  then  placed  in  a  calcium-buffered  antibacterial  saline 
solution.  Storage  in  such  solutions  at  concentrations  previ¬ 
ously  determined  to  maintain  equilibrium  has  been  found  to 
maintain  bone  mechanical  properties  by  minimizing  the  dif¬ 
fusion  of  calcium  from  the  tissue  [27],  The  specimens  were 
stored  at  4°C  until  needed. 

Indentations 

Cylindrical  indentations  were  made  with  a  custom-built  cy¬ 
lindrical  indenter  tip  attached  to  an  fnstron  servo-hydraulic 
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Fig.  2.  Advantages  of  factor  analysis:  signal-to-noise  im¬ 
provement  and  background  removal.  (A)  Raw  Raman  spec¬ 
trum  of  bone  from  one  of  the  image  data  sets  discussed  in  this 
work.  (B)  Factor,  representing  the  Raman  spectrum  of  bone, 
derived  from  factor  analysis  of  the  entire  image  data  set  (ap¬ 
proximately  3100  spectra  total).  For  clarity,  additional  back¬ 
ground  factors  resulting  from  luminescence  are  not  shown. 

mechanical  testing  machine  (model  8521).  Using  a  cylindrical 
indenter  is  useful  for  spectroscopic  imaging  purposes  because 
the  floor  of  the  indent  is  flat,  allowing  spectroscopic  investi¬ 
gation  within  the  indent  as  well  as  at  its  edges.  Indentations 
resulting  from  two  different  loads  were  made  upon  each 
specimen.  The  shallower  (low  load)  indents  (approximate 
depth  10  pm)  were  made  with  loads  ranging  from  70  to  110  N, 
and  the  loads  for  the  deeper  (high  load)  indents  (approximate 
depth  100-150  pm)  were  in  the  range  of  230-240  N. 

Raman  Imaging 

Our  Raman  imaging  system  (Fig.  1)  has  been  previously  de¬ 
scribed  [14,  16,  28,  29].  Briefly,  it  consists  of  a  785  nm  external 
cavity  diode  laser  (Kaiser  Optical  Systems,  Inc.)  with  a  high- 
aspect-ratio  rectangular  beam  profile.  The  beam  size  is  reduced 
using  a  telescope  and  the  beam  is  coupled  into  a  modified  epi- 
fluorescence  microscope  (Olympus,  Inc.).  A  1 0x/0.5  NA  ob¬ 
jective  (Zeiss)  focuses  the  line-shaped  laser  beam  (length  ~200 
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pm)  onto  the  specimen.  A  complete  Raman  spectrum  is  col¬ 
lected  simultaneously  front  every  point  along  the  line  by  the 
objective,  focused  by  a  collection  lens,  and  directed  into  an 
axial-transmissive  spectrograph  (Kaiser  Optical  Systems,  Inc.) 
which  disperses  the  Raman  scatter.  The  signal  is  then  detected 
by  a  thermo-electrically  cooled  CCD  camera  (Andor  Tech¬ 
nologies,  Inc.)  and  stored  on  a  personal  computer.  A  two- 
dimensional  image  is  built  up  by  moving  the  specimen  in 
approximately  1  pm  steps  under  the  laser  beam  using  a 
translation  stage  (New  England  Affiliated  Technologies,  Inc.) 
and  acquiring  a  new  line  of  spectra  at  each  new  stage  position. 
An  image  data  set  thus  consists  of  a  matrix  with  every  point 
containing  a  complete  Raman  spectrum. 

For  these  studies,  Raman  images  were  taken  in  three  lo¬ 
cations  for  each  indent:  within  the  indent  on  the  indent  floor, 
at  the  very  edge  of  the  indent  on  the  bone  surface,  and  at  least 
500  pm  outside  the  indent  area.  A  total  of  18  images  (3  images/ 
indent  x  2  indents/specimen  x  3  specimens)  were  obtained. 
Each  image  contained  25  lines  of  spectra;  the  approximate  size 
of  each  area  imaged  was  200  x  35  pm.  The  total  laser  power  at 
the  specimen  was  approximately  230  mW,  spread  into  a  rect¬ 
angle  approximately  200  pm  long  and  6  pm  wide.  It  should  be 
noted  that  because  of  the  shape  of  the  laser  beam,  power 
density  is  relatively  low,  and  heat  dissipation  is  highly  efficient 

[30] .  Heat  damage  to  mineralized  tissue  specimens  of  the  size 
used  in  this  study  has  not  been  observed  with  this  technique 

[31] . 


Data  Analysis 

Both  univariate  and  multivariate  analysis  techniques  were  used 
to  examine  the  data.  A  spectroscopic  parameter  proportional 
to  the  mineral-to-organic-matrix  ratio  was  used.  Similar 
parameters  have  been  frequently  employed  by  infrared  spec- 
troscopists  [13,  32],  This  spectroscopic  mineral-to-organic- 
matrix  ratio  was  calculated  for  all  of  the  areas  imaged  by 
dividing  the  spectra  from  each  image  into  10  random  groups, 
calculating  the  mean  spectrum  from  each  group,  fitting  this 
mean  spectrum  to  a  third-order  polynomial  to  remove  the 
effects  of  luminescent  background,  then  dividing  the  maximum 
height  of  the  phosphate  V!  band  at  ~959  cm-1  by  the  maxi 
mum  height  of  the  CH2  wag  band  at  ~1450  cm-1.  The  CH2 
wag  band  was  chosen  since  it  is  less  likely  than  the  amide 
bands  to  change  shape  or  shift  due  to  changes  in  protein 
structure  caused  by  applied  load  [33],  Ten  mineral-to-organic- 
matrix  ratio  values  were  generated  for  each  image,  and  a 
standard  deviation  was  calculated  from  these  10  values. 

Factor  analysis  was  also  performed  on  each  image  data  set. 
Factor  analysis  is  a  well-established  multivariate  image  anal¬ 
ysis  technique  [34]  used  in  such  fields  as  magnetic  resonance 
imaging  [35—37]  and  fluorescence  imaging  [38],  We  have  fre¬ 
quently  used  this  data  analysis  approach  for  the  analysis  of 
Raman  spectroscopy  and  image  data  [39-41].  In  this  tech¬ 
nique,  the  data  are  considered  as  matrices  with  one  spatial  and 
one  spectral  dimension.  The  covariance  matrix  of  this  matrix  is 
calculated,  and  a  singular  value  decomposition  is  performed  to 
obtain  the  eigenvalues  and  eigenvectors  of  the  covariance 
matrix.  The  eigenvectors  can  then  be  combined  in  linear 
combinations  (or  “rotated”)  using  non-negativity  and  band- 
shape  constraints  to  produce  Raman  spectra.  This  is  not  a 
random  process.  Spectra  cannot  be  created  by  rotation;  how¬ 
ever,  they  can  be  separated  from  backgrounds  very  effectively. 
In  addition,  because  factor  analysis  is  an  averaging  technique, 
the  signal-to-noise  ratio  in  factors  is  generally  significantly 
better  than  in  raw  spectra  (Fig.  2).  Factor  analysis  is  useful  for 
separating  species  that  may  have  overlapping  spectra,  but 
differ  in  their  spatial  distributions.  Thus,  if  two  species  have 
different  spectra  but  identical  spatial  distributions,  their  spec¬ 
tra  will  not  be  separated. 

Once  the  eigenvectors  (now  called  factors)  have  been  ef¬ 
fectively  rotated,  each  individual  factor’s  distribution  in  the 
area  surveyed  can  be  calculated  by  regressing  the  factors  back 
onto  the  data.  The  resulting  image  is  called  a  score  image  and 
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Fig.  3.  SEM  images  of  cylindrical  indents.  (A)  Low-load  in¬ 
dent  on  specimen  #3.  Note  the  shear  bands  visible  in  the  high- 
magnification  image  on  the  right.  (B)  High-load  indent  on  the 
same  specimen.  Again,  shear  bands  are  visible  in  the  higher- 

is  a  plot  of  the  spatial  distribution  of  factor  intensity  within  the 
area  imaged. 

For  these  experiments  the  image  data  were  divided  into  two 
spectral  subregions:  a  “mineral  subregion’’  (approximately 
500-1200  cm-f)  containing  the  most  intense  and  characteristic 
bands  of  the  mineral  components,  and  an  “organic  matrix 
subregion”  (approximately  1200-1800  cm-1)  containing  the 
most  intense  bands  attributable  to  the  collagenous  matrix. 
Each  spectral  subregion  was  analyzed  separately.  Factor 
analysis  was  performed  in  MATLAB  software  (Mathworks, 
Inc.);  both  vendor-supplied  and  locally  written  scripts  were 
used. 


SEM  Imaging 

Indentation  specimens  were  prepared  for  SEM  analysis  with  a 
thin  coating  of  colloidal  carbon  for  electron  conductivity.  A 
Hitachi  S3200N  Scanning  Electron  Microscope  was  used  for 
all  imaging  under  conventional  high-vacuum  mode  and  sec¬ 
ondary  electron  scintillator  detection  mode.  The  electron  beam 
energy  was  kept  at  20  kV  and  the  operating  distance  was  ad¬ 
justed  to  15  mm.  Digital  images  were  captured  and  stored  with 
Quartz  PCI  acquisition  software. 


magnification  image.  Shear-induced  radial  cracking  is  also 
visible  in  both  lower-magnification  images.  (C)  High-magnifi¬ 
cation  view  of  shear  bands  along  the  indent  wall. 


Results 

SEM  Images 

Images  of  both  low-  and  high-load  indents  clearly 
showed  shear  bands  around  the  circumference  of  the 
indent,  indicating  the  presence  of  shear  forces  caused  by 
the  indenting  process  (Fig.  3). 

Miner al-t o-Organic-Matrix  Ratio 

The  spectroscopic  mineral-to-organic-matrix  ratio  for 
image  areas  was  calculated  in  order  to  confirm  that  the 
areas  imaged  were  comparable  in  terms  of  composition, 
and  that  any  changes  we  observed  were  a  result  of  me¬ 
chanical  loading  and  deformation,  not  compositional 
variations.  Because  Raman  band  intensity  varies  linearly 
with  the  number  of  scattering  centers,  the  phosphate 
Vi/CH2  wag  band  height  ratio  is  proportional  to  the 
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Fig.  4.  Bar  graph  of  spectroscopic  mineral-to-organic-matrix  ratio  (phosphate  V!  maximum  intensity  divided  by  CH2  wag  band 
maximum  intensity)  calculated  for  all  areas  imaged.  Error  bars  are  ±  1  SD. 


chemical  mineral-to-organic-matrix  ratio.  Figure  4 
shows  a  bar  graph  of  the  calculated  mineral-to-organic- 
matrix  values  for  all  the  areas  imaged.  The  values  are 
relatively  consistent.  Two-way  ANOVA  revealed  that 
while  there  were  no  statistically  significant  differences 
between  any  of  the  areas  on  specimens  2  and  3 
(_P<0.05),  specimen  1  was  significantly  different  from 
the  other  two  specimens  in  the  edge  and  indent  areas.  It 
should  be  noted  that  there  were  no  statistically  sig¬ 
nificant  differences  in  the  spectroscopic  mineral-to- 
organic-matrix  ratio  among  the  widely  spaced  control 
areas  on  any  of  the  specimens  [P<  0.207).  Statistical 
variations  in  the  indent  and  edge  regions  therefore  could 
be  the  result  of  damage  resulting  from  the  indent  pro¬ 
cess.  Based  on  these  results,  we  concluded  that  the  dif¬ 
ferences  in  spectroscopic  mineral-to-organic-matrix 
ratio  among  the  three  specimens  were  not  significant, 
and  that  comparisons  among  the  different  specimens 
were  valid. 

Factor  Analysis:  Organic  Matrix  Subregion 

In  all  images  of  control  areas  and  areas  within  inden¬ 
tations,  only  one  organic  matrix  factor  was  observed 
(Fig.  5).  This  factor  contained  the  bands  commonly 
observed  in  Raman  spectra  of  bone  organic  matrix 
[10]:  the  amide  III  band  at  ~1250  cm-1,  the  CH2  wag  at 
~1447  cm"1,  and  the  amide  I  band  at  ~1667  cm-1.  Score 
images  associated  with  this  single  factor  showed  a 


relatively  uniform  distribution  throughout  the  imaged 
area. 

Different  results  were  obtained  from  the  images 
taken  at  the  edge  of  the  indents.  Here,  two  organic 
matrix  factors  were  observed  (Fig.  6).  One  represented 
the  commonly  observed  Raman  spectrum  of  bone  or¬ 
ganic  matrix  described  above.  The  second  factor  dis¬ 
played  similar  bands,  but  the  amide  III  band  showed  a 
distinct  increase  in  its  lower-frequency  component, 
whereas  the  amide  I  band  was  shifted  to  higher  wave- 
numbers  (~1678  cm-1).  The  score  images  of  these  two 
factors  were  complementary,  with  the  first  factor  having 
higher  intensity  closer  to  the  indent  and  the  second 
having  higher  intensity  further  away  from  the  indent. 

Factor  Analysis:  Mineral  Subregion 

Multiple  mineral  factors  containing  the  phosphate  Vj 
peak  in  different  positions  were  observed  in  most  of  the 
areas  imaged  (Fig.  7).  It  is  not  uncommon  to  observe 
more  than  one  mineral  factor  in  Raman  data  from  a 
given  bone  specimen  owing  to  tissue  heterogeneity  re¬ 
sulting  from  the  aging  and  remodeling  processes  [42]. 
Slight  local  variations  in  the  chemical  composition  of  the 
mineral  result  in  Raman  spectra  with  varied  phosphate  Vj 
band  positions,  which  are  easily  separated  by  factor 
analysis.  However,  an  interesting  trend  was  observed  in 
the  entire  set  of  images.  In  almost  every  image  group, 
more  mineral  factors  were  observed  within  the  indented 
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Fig.  5.  Images  of  control  regions  and  regions  with  indenta¬ 
tions  for  the  high-load  indent  on  specimen  #3.  (A)  Brightfield 
image  showing  the  control  area  imaged  for  the  230  N  indent  on 
specimen  #3  (the  indent,  which  is  approximately  500  pm  away, 
is  out  of  the  field  of  view).  The  white  box  indicates  the  area 
imaged.  To  the  right  is  the  single  organic  matrix  factor  ob- 


■ 


low  high 

tained  for  this  image  together  with  its  score  image.  (B)  A  single 
organic  matrix  factor  is  also  obtained  for  the  area  imaged 
within  the  indent.  Image  size  in  both  cases  is  approximately 
200  x  35  pm.  For  clarity,  the  additional  background  factors 
resulting  from  luminescence  are  not  shown. 


area  than  at  the  edge  of  the  indent  or  in  the  control  areas. 
Table  1  shows  these  results,  with  the  bands  correspond¬ 
ing  to  new  factors  that  appeared  after  indentation  under¬ 
lined.  In  most  cases  a  lower-frequency  calcium  phosphate 
component  appeared;  in  one  case,  a  higher-frequency 
calcium  phosphate  component  appeared.  This  trend  is 
consistent  with  other  studies  from  this  laboratory  show¬ 
ing  an  increase  in  the  number  of  mineral  factors  present 
in  Raman  data  collected  from  areas  of  damage  [14,  43]. 

Discussion 

We  attribute  the  appearance  of  the  additional  collagen 
factor  to  damage.  This  factor  containing  shifted  amide 


bands  is  indicative  of  the  presence  of  collagen  that  has 
undergone  a  transformation  from  its  normal  triple-he¬ 
lical  state  to  a  disordered  state,  presumably  as  a  result  of 
the  loading  and  deformation  induced  by  the  indentation 
process.  The  appearance  of  this  factor  is  consistent  with 
observations  made  by  other  researchers  on  protein  films 
and  protein  solutions  [19].  The  higher-frequency  amide  I 
band  and  the  lower-frequency  amide  III  band  observed 
in  this  second  collagen  factor  are  known  spectral  sig¬ 
natures  of  collagen  in  which  crosslinks  have  been  re¬ 
duced  or  broken  [44];  similar  bands  have  also  been 
observed  in  studies  of  collagen  denaturation  [45-47]. 
The  breaking  of  crosslinks  in  mature  tissue  is  thought  to 
be  the  last  process  that  happens  before  fracture  [8]. 
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Fig.  6.  Brightfield  image  showing  the  edge  of  the  high-load 
indent  oil  specimen  #3.  The  white  box  indicates  the  area  im¬ 
aged.  To  the  right  are  the  two  organic  matrix  factors  obtained 
for  this  image  together  with  their  score  images.  The  factor  on 
the  top  is  the  Raman  spectrum  of  normal  collagen.  The  factor 
on  the  bottom  represents  a  Raman  spectrum  of  disorganized 
collagen,  showing  an  increase  in  the  low  frequency  component 


of  the  amide  III  band  and  high  frequency  component  of  the 
amide  1  band,  typical  of  collagen  that  has  lost  its  structure 
through  crosslink  rupture.  Note  that  the  score  images  are 
complementary.  Image  size  in  both  cases  is  approximately  200 
x  35  pm.  For  clarity,  the  additional  background  factors 
resulting  from  luminescence  are  not  shown. 


It  is  interesting  that  this  damaged  collagen  factor  is 
only  observed  at  the  edge  of  the  indents,  not  within  the 
indent  area  itself.  We  hypothesize  that  at  the  edge  of  the 
indents,  the  collagen  fibrils  experience  shear  forces 
caused  by  the  edge  of  the  indenter  pushing  into  the 
bone.  These  shear  forces  cause  the  rupture  of  crosslinks. 
Shear  behavior  in  the  area  of  indentations  is  well  char¬ 
acterized.  For  the  case  of  a  cylindrical  indenter,  a  simple 
analytical  description  of  the  elastic  and  plastic  stress 
field  below  the  punch  was  given  by  Sneddon  [48].  Unlike 
sharp  indenters,  friction  is  important  for  indenters  that 
are  flat  or  rounded.  With  friction,  the  pressure  needed  to 
penetrate  to  any  given  depth  increases,  and  shear  failure 
along  the  indenter/nraterial  interface  nucleates  cracks 
which  travel  outward  from  the  compressive  zone  directly 
under  the  indenter  [49],  These  cracks  are  clearly  visible 
in  the  SEM  image  of  the  high-load  indent  in  Figure  4.  In 


addition,  closer  observation  of  the  SEM  images  shows 
the  presence  of  a  shear  banding  pattern  at  the  indent 
edges.  Similar  shear  banding  patterns  have  been  ob¬ 
served  at  the  edges  of  spherical  indentations  [50],  and 
shear  banding  and  pile  up  at  the  edge  of  an  indentation 
have  also  been  generated  by  foreign  object-induced 
damage  in  metal  alloys  [51].  In  all  these  cases  the  ma¬ 
terial  was  exposed  to  an  indenting  process  that  allowed 
shearing  to  occur  under  the  indent.  In  the  case  of  the 
cylindrical  punch,  the  complete  region  below  the  indent 
is  subject  to  a  compressive  force  limiting  the  appearance 
of  local  shear  bands  to  just  outside  the  pile-up  zone  at 
the  edge  of  the  indent.  The  fact  that  the  damaged  col¬ 
lagen  factor  does  not  appear  in  the  data  acquired  from 
the  middle  of  the  indented  area  suggests  that  the  colla¬ 
gen  in  this  area  is  not  damaged  in  a  way  that  is  spec¬ 
troscopically  visible;  i.e.,  it  is  compressed  but  crosslinks 
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Raman  shift,  cm-1 


Raman  shift,  cnr1 

Fig.  7.  Example  of  multiple  mineral  factors  observed  in  in¬ 
dented  specimens.  (A)  Single  factor  obtained  from  analysis  of 
the  control  area  on  specimen  #3.  (B)  Factors  representing 
multiple  mineral  factors  obtained  from  analysis  of  the  high- 
load  indent  on  specimen  #3.  For  clarity,  the  additional  back¬ 
ground  factors  resulting  from  luminescenceare  not  shown. 

are  not  broken.  This  is  consistent  with  the  stress  be¬ 
havior  described  for  cylindrical  indentations. 

The  appearance  of  additional  mineral  factors  within 
the  indented  area  may  indicate  a  pressure-induced 
structural  transition  of  the  mineral.  The  addition  of  ei¬ 
ther  a  higher-  or  lower-frequency  component  is  consis¬ 
tent  with  high-pressure  spectroscopy  studies  on  calcium 
hydroxyapatite,  which  show  an  increase  in  frequency  of 
the  phosphate  V)  band  followed  by  a  decrease  with  in¬ 
creasing  pressure  [23,  25],  The  loads  applied  in  this  ex¬ 
periment  result  in  stresses  ranging  from  0.4  GPa,  yielding 
the  shallow  indents,  to  1.2  GPa,  yielding  the  deep  in¬ 
dents;  these  values  are  comparable  to  those  applied  in  the 
high-pressure  spectroscopy  studies  in  which  a  structural 
transition  was  observed.  We  therefore  believe  that  the 
additional  factors  represent  spatially  distinct  mineral 
species  formed  as  a  consequence  of  pressure-induced 
structural  transitions  within  the  bone  mineral  lattice  [23, 
25].  We  have  earlier  observed  the  appearance  of  addi¬ 
tional  mineral  factors  accompanying  permanent  defor¬ 
mation  in  bone;  for  example,  in  the  macroscopic  fracture 
of  mouse  femora  [43]  as  well  as  in  microcracks  within 
bovine  bone  [14].  However,  in  these  studies  we  could  not 
establish  whether  the  additional  mineral  factors  observed 
were  a  cause  or  an  effect  of  the  fracture  process.  The 
results  from  the  current  experiments  indicate  that  the 
additional  mineral  factors  are  in  fact  an  effect  of  me¬ 
chanical  loading.  The  fact  that  we  can  resolve  multiple 
factors  rather  than  seeing  one  mineral  factor  with  a 


shifted  phosphate  V!  band  indicates  that  not  all  the 
mineral  present  undergoes  the  transition,  possibly  be¬ 
cause  the  local  pressure  was  not  high  enough  to  drive  the 
transition  to  completion,  or  because  not  all  the  mineral 
species  present  were  able  to  undergo  the  transition. 

The  nature  of  the  structural  transition  occurring  in 
the  bone  mineral  is  not  known.  Xu  et  al.  [23]  postulate 
that  the  structural  transition  occurring  in  pure  hy¬ 
droxyapatite  is  the  result  of  a  rearrangement  of  hydroxyl 
groups.  However,  bone  mineral  is  thought  not  to  con¬ 
tain  hydroxyl  groups  in  any  appreciable  concentration 
[52-54].  In  addition,  the  additional  substituents  in  bone 
mineral  undoubtedly  affect  the  response  of  the  mineral 
to  applied  load.  It  is  known  that  carbonation  in  bone 
mineral  causes  vacancies  in  the  mineral  lattice,  creating 
an  opportunity  for  pressure-induced  movement  and 
structural  changes  [53].  Further  spectroscopic  and  X-ray 
diffraction  studies  are  required  to  better  elucidate  the 
nature  of  this  transformation. 

Conclusions 

Spectroscopic  changes  resulting  from  mechanical  load¬ 
ing  of  bone  were  observed  in  these  studies,  indicating  two 
different  damage  mechanisms  for  different  bone  tissue 
components.  The  organic  matrix  showed  clear  evidence 
of  crosslink  rupture,  consistent  with  what  is  known  of  the 
mechanical  behavior  of  other  collagenous  tissues,  par¬ 
ticularly  more  mature  tissues  which  contain  more 
crosslinks  [8,  9],  This  rupture  occurred  only  in  areas 
where  shear  forces  were  present  (indent  edges),  rather 
than  where  only  compression  occurred  (indent  floor). 
The  mineral  phase  of  the  tissue,  on  the  other  hand,  was 
found  to  be  altered  where  only  compression  had  oc¬ 
curred.  The  possibility  of  a  pressure-induced  structural 
transition  occurring  in  bone  mineral  is  particularly  in¬ 
triguing.  This  could  be  a  method  by  which  energy  is 
dissipated  with  minimal  damage  to  the  surrounding  tis¬ 
sue.  Similar  behavior  has  been  observed  in  ceramic  ma¬ 
terials  subjected  to  stress,  a  phenomenon  known  as 
“transformation  toughening’’  [55,  56].  These  processes 
have  also  recently  been  postulated  to  occur  in  bone  [57], 
We  propose  that  the  bone  mineral  structural  transitions 
that  we  have  observed  are  the  source  of  this  transfor¬ 
mation  toughening.  Aging,  as  well  as  bone  diseases  that 
alter  bone  mineral  crystallite  structure,  such  as  osteo¬ 
porosis,  could  hinder  this  transition  and  contribute  to 
increased  fragility.  If  this  is  the  case,  new  and  different 
treatment  options  may  be  realized  from  these  results. 
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Table  1.  Number  of  mineral  factors,  and  phosphate  Vi  band  positions  for  those  factors,  observed  for  all  specimens.  Factors  that 
appeared  after  indentation  (relative  to  the  control)  are  underlined. 


#1,  low  load 

#1,  high  load 

#2,  low  load 

#2,  high  load 

#3,  low  load 

#3,  high  load 

Control 

2  species 

955 

3  species 

954 

2  species 

959 

1  species 

959 

1  species 

959 

1  species 

959 

959 

959 

965 

962 

Edge  of  Indent 

Not  observed 
due  to  fluorescence 

1  species 

960 

2  species 

956 

960 

1  species 

960 

1  species 

960 

1  species 

960 

Within  Indent 

3  species 

954 

2  species 

955 

3  species 

954 

2  species 

953 

2  species 

957 

2  species 

957 

957 

959 

958 

959 

960 

959 

962 

962 
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ABSTRACT 

Osteogenesis  imperfecta  (OI)  is  genetic  defect  in  which  the  genes  that  code  for  the  at (I)  or  a2(I)  chains  of  type  I 
collagen  are  defective.  The  defects  often  result  in  substitution  of  a  bulky  amino  acid  for  glycine,  causing  formation  of 
collagen  that  can  not  form  the  normal  triple  helix.  Depending  on  the  details  of  the  defects,  the  outcomes  range  from 
controllable  to  lethal.  This  study  focuses  on  OI  type  IV,  a  more  common  and  moderately  severe  form  of  the  disease. 
People  with  the  disease  have  a  substantial  increase  in  the  risk  and  rate  of  fracture.  We  examine  the  spectroscopic 
consequences  of  these  defects,  using  a  mouse  model  (BRTL)  that  mimics  OI  type  IV.  We  compare  Raman  images  from 
tibial  cortical  tissue  of  wild-type  mice  and  BRTL  mice  with  single  copy  of  mutation  and  show  that  both  mineral  to 
matrix  ratios  and  collagen  inter-fibril  cross-links  are  different  in  wild-type  and  mutant  mice. 

Keywords:  Raman  spectroscopy,  Raman  Imaging,  Osteogenesis  Imperfecta,  bone  tissue,  collagen,  cross-links 


1.  INTRODUCTION 

Bone  is  composed  of  two  major  components:  inorganic  mineral  and  an  organic  matrix:  About  90%  of  the  organic  matrix 
is  type  I  collagen.  Collagen  provides  a  supporting  matrix  on  which  mineralization  occurs.  The  mineral  is  a  highly 
carbonated  form  of  apatite,  where  C032"  may  substitute  either  for  P043"  (B-type  carbonation)  or  OH'  (A-type 
carbonation).  A  variety  of  other  anions  and  cations  may  also  substitute  for  sites  of  the  crystal  lattice111. 

Osteogenesis  imperfecta  (OI),  also  known  as  brittle  bone  disease,  is  a  genetic  defect  affecting  bone  matrix  collagen1"1. 
As  a  result  of  mutations  in  one  of  the  two  genes  (COL1 A1  and  COL1 A2)  that  encode  the  strands  of  type  I  collagen,  risk 
of  fracture  significantly  increases.  Collagen  is  a  heterotrimer,  composed  of  two  copies  of  the  al  chain  and  one  copy  of 
the  a2  chain.  Each  chain  is  a  string  of  triplets  containing  glycine  followed  by  two  other  amino  acids,  frequently  proline 
and  hydroxyproline131.  The  mutations  causing  OI  most  commonly  result  from  a  substitution  of  glycine  by  a  larger  amino 
acid,  disrupting  the  triple  helix.  This  alteration  in  the  triple  helix  structure  may  alter  normal  collagen  cross  linking  and 
matrix  mineralization,  depending  on  the  location  of  the  defect  along  the  collagen  chain(s). 

Four  types  of  OI  have  been  attributed  to  collagen  mutations,  and  range  from  mild  (Type  I)  to  perinatal  lethal  (Type  II) [4]’ 
[51.  In  this  study  we  investigate  tissue  from  the  BRTL  mouse,  a  model  for  OI  type  IV,  a  moderately  severe  but  survivable 
form  of  the  disease.  BRTL  has  an  al  (I)  Gly  349  -  Cys  substitution  in  type  I  collagen,  mimicking  the  genetic  mutation 
found  in  an  OI  type  IV  patient  (Forlino  et  al,  1999)161191.  Just  as  humans  with  OI  get  fewer  fractures  as  they  age, 
biomechanical  tests  show  that  at  two  months  of  age,  BRTL  mice  have  lower  bone  strength  and  stiffness  than  wild-type 
mice,  but  at  six  months,  these  parameters  improve  to  WT  levels.  This  alteration  is  accomplished  without  corresponding 
improvements  in  the  geometry  of  the  bone,  implying  that  the  changes  are  due  to  improvements  in  the  material  properties 
of  the  bone  matrix  (Kozloff  et  al,  2004)[71’[8I'[I01. 
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Raman  microspectroscopy  has  matured  into  an  important  complement  to  traditional  methods  for  detection, 
quantification  and  imaging  of  local  variations  in  the  molecular  structures  of  bone  matrix  and  mineral.  While  similar 
information  is  available  from  infrared  spectroscopy,  Raman  spectroscopy  is  equally  applicable  to  specimens  as  thin  as  a 
few  microns  or  as  thick  as  intact  tissue  and  offers  superior  spatial  resolution.  As  a  microspectroscopy,  Raman  scattering 
provides  micron-level  spatial  resolution  and  can  be  used  to  generate  images  that  map  chemical  structure  or  physical 
property  variations. 


2.  APPARATUS  AND  METHODS 

2.1  Raman  microspectroscopic  imaging  system 

Our  epi-illumination  Raman  microspectroscopic  imaging  system  with  a  line-focused  laser  beam  is  shown  in  Figure  1. 
We  focus  an  external  cavity  785  nm  diode  laser  onto  the  specimen  through  a  10X/  0.5  NIR-coated  objective.  Raman 
scatter  is  collected  through  the  same  objective  and  focused  onto  the  entrance  slit  of  a  NIR  optimized  spectrograph  fitted 
with  a  cooled  back-illuminated  CCD.  By  using  a  line-focused  beam,  we  can  simultaneously  collect  a  line  of  254  spectra, 
and  then  by  translating  the  specimen  under  the  objective  in  small  increments,  an  image  is  built  up  line-by-line. 


Thermo-electrically  cooled 
CCD  detector 


Figure  1 .  Raman  microspectroscopic  imaging  system 


2.2  Data  analysis 

The  data  obtained  is  analyzed  using  both  univariate  analysis  and  multivariate  analysis  (i.e.,  band  target  entropy 
minimization,  or  BTEM)[111' [121.  In  BTEM,  singular  value  decomposition  is  used  to  calculate  the  eigenvectors  of  the 
spectroscopic  image  data  set.  Bands  in  non-noise  eigenvectors  that  would  normally  be  used  for  recovery  of  spectra  are 
examined  for  localized  spectral  features.  For  a  targeted  (identified)  band,  information  entropy  minimization  or  a  closely 
related  algorithm  is  used  to  recover  the  spectrum  containing  this  feature  from  the  non-noise  eigenvectors,  plus  the  next 
5-30  eigenvectors,  in  which  noise  predominates.  By  adding  or  subtracting  the  recovered  spectra,  known  as  “factors,”  in 
linear  combinations  and  regressing  them  back  into  the  original  data,  the  distribution,  or  “score,”  of  each  factor  in  the 
imaged  data  can  be  obtained.  The  univariate  data  analysis  is  performed  using  GRAMS32,  and  the  multivariate  data 
analysis  is  performed  using  MATLAB  (Version  6.1,  Mathworks,  Inc.). 


3.  EXPERIMENTS 

Tissue  from  previously  sacrificed  2  month  male  mice  with  single  copies  of  the  mutation  (BRTL,  n=5),  and  wild-type 
mice  (WT,  n=5)  were  employed  in  this  study.  Tibias  were  prepared  by  embedding  in  a  weakly  exothermic  polymer. 


86  Proc.  of  SPIE  Vol.  5321 


sectioning  at  the  midshaft  using  a  low-speed  diamond  saw,  and  polishing  to  near  the  tibio-fibular  junction.  Five 
transects,  each  containing  a  line  of  254  Raman  spectra,  were  taken  along  each  transverse  cross  section,  extending  from 
the  endosteal  to  the  periosteal  surface  of  the  posterior  aspect  of  the  tibia.  Following  baselining,  areas  under  the 
phosphate  Ui  (903-991  cm’1),  carbonate  Ui  (1046-1110  cm’1),  and  amide  I  (1616-1720  cm"1)  bands  were  measured. 
Phosphate  to  amide  I  band  area  ratio  was  calculated  as  a  measure  of  mineral  maturity  and  mineral  to  matrix  ratio. 
Additionally,  to  study  the  chemical  structure  variations,  a  254  x  150  pixel  image  of  an  approximately  330  x  200pm2  was 
taken  from  the  same  region  of  each  specimen.  After  recovery  of  the  bone  spectra  from  the  image  data,  curve-fitting  was 
used  to  resolve  the  underlying  bands  of  amide  I,  and  areas  under  the  amide  1  sub-bands  bands  were  obtained. 


4.  RESULTS  AND  DISCUSSION 

4.1  Abnormal  mineral  to  matrix  ratio  in  BRTL  mouse 

Two  month  old  BRTL  mice  show  approximately  15%  lower  mineral/matrix  ratios  than  those  of  the  wild-type  mice  at 
the  same  age.  As  shown  in  Figure  2  the  mineral/matrix  ratio  for  the  wild  type  mice  is  5.44  +  0.402,  while  in  the  BRTL  it 
is  4.64  ±  0.180.  The  lower  mineral  to  matrix  ratio  of  the  BRTL  mice  suggests  impaired  mineral  deposition  on  the 
defective  matrix,  caused  by  a  single  glycine  to  cysteine  point  substitution. 


Mineral  to  matrix  ratio 


Figure  2.  Mineral  to  matrix  ratio  of  2  month  old  WT  and  BRTL  mice 


4.2  Spectral  differences  and  separated  mineral  factors  found  in  the  BRTL  mouse 

In  the  resulting  data,  we  have  found  some  spectral  differences  between  BRTL  and  wild-type  mice.  The  Amide  I  band  of 
BRTL  mouse  slightly  shifts  to  higher  frequency,  while  the  Amide  III  band  shifts  to  lower  frequency.  The  spectral 
difference  on  the  matrix  is  what  we  would  expect  to  see  since  the  glycine  to  cysteine  substitution  disrupts  the  collagen 
triple  helix  structure  and  may  cause  collagen  to  unwind.  The  Raman  spectrum  of  bone  mineral  has  a  prominent 
phosphate  band  at  958-960  cm"1.  In  addition  to  the  spectral  differences  observed  on  the  amide  I  band,  band  target 
entropy  minimization  extracted  a  mineral  factor  at  952  cm'1,  which  likely  reflects  highly  disordered 
crystallites11311141  .While  A-type  carbonation  is  known  to  result  in  a  highly  deformed  lattice  bsmi  ],[i7]^  ^  wouj^  be  a 
simplification  to  attribute  this  band  to  crystallites  containing  A-type  carbonate.  The  mineral  factor  at  964  cm’1  is 
uncarbonated  phosphate,  implying  a  less  mature  tissue.  Figure  3  shows  the  three  mineral  factors  and  their  corresponding 
score  images  of  a  BRTL  mouse.  Background  factors  are  omitted  for  clarity.  The  images  give  the  spatial  distribution  of 
these  factors  in  the  imaged  area.  The  multiple  mineral  factors  found  in  the  BRTL  mouse  indicates  that  the 
mineralization  of  the  tissue  is  more  heterogeneous  than  that  of  the  wild-type  mouse. 
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Figure  3.  Factors  and  corresponding  score  images  of  a  BRTL  mouse,  (a)  the  major  bone  factor  with  phosphate 
at  958  cm'1,  carbonate  at  1070  cm'1,  amide  III  at  1245  cm'1,  CH2  wag  at  1450  cm'1,  and  amide  I  at  1665  cm'1. 

(b)  the  second  mineral  factor  with  phosphate  at  952  cm'1,  (c)  the  third  mineral  factor  with  phosphate  at  964  cm'1. 


4.3  Collagen  cross-links 

An  important  strength  of  both  Raman  and  infrared  spectroscopy  is  sensitivity  to  protein  secondary  structure.  Amide  I  is 
usually  an  envelope  that  can  be  resolved  into  two  or  more  bands  and  that  can  extend  over  much  of  the  1600-1700  cm"1 
region.  Similar  information  is  carried  in  amide  II  and  amide  III,  although  they  are  less  frequently  used  to  elucidate 
secondary  structure. 

Collagen  is  a  special  case  because  of  its  unusual  triple  helix  structure.  There  are  structural  differences  among  the  various 
types  of  collagen,  of  course,  but  the  general  patterns  are  the  same  throughout.  While  it  is  reasonable  to  assume  that 
collagen  amide  I  should  resemble  an  a-helix  spectrum,  the  actual  spectrum  is  more  complex.  In  uncrosslinked  collagen 
(e.g.  from  rat  or  calf  skin,  a  triple  helix  of  cq  chains)  the  amide  I  region  can  be  resolved  into  several  bands.  The  most 
prominent  in  the  FTIR  spectra  are  at  1630,  1660,  and  1670  cm"1,  although  these  bands  are  shifted  somewhat  depending 
on  the  type  of  collagen.  Based  on  studies  of  poly(gly-pro-pro),  Lazarev  and  co-workers  assign  these  bands  to  the 
carbonyl  stretches  of  hydrogen-bonded  glycine,  proline -2  and  proline-3  in  a  gly-pro-pro  motif.1201  However,  these 
assignments  must  be  considered  approximate  because  they  do  not  assume  the  presence  of  hydroxyproline,  let  alone 
other  amino  acids.  Payne  and  Veis  attribute  a  rat  skin  collagen  band  at  1660  cm"1  to  carbonyls  of  the  X  amino  acid 
hydrogen  bonded  to  glycine  amide  protons  of  an  adjacent  chain.  Additionally,  the  simplified  assignments  do  not  take 
cognizance  of  the  effects  of  interfibril  cross-linking  on  the  secondary  structure  of  collagen  and  the  consequent  changes 
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in  the  amide  I  envelope.  The  high  frequency  component  at  1685-1695  cm"1  is  associated  with  disruption  of  hydrogen 
bonding,  although  the  origin  of  the  effect  remains  unclear. 

Inter-  and  intra-fibril  cross-links  are  a  well-known  feature  of  bone  tissue  collagen.  Disruption  or  absence  of  these  cross¬ 
links  results  in  changes  in  the  relative  intensities  of  the  major  bands  of  amide  I.  The  disruptions  can  include  chemical 
perturbation  of  cross-links  and  their  chemical  cleavage  as  well  as  mechanical  rupture.  Similar  effects  are  observed  in 
immature  tissue,  where  the  crosslinks  have  not  yet  been  formed. 

It  has  been  shown  that  denaturing  collagen  to  gelatin  causes  a  high  frequency  shoulder  at  1685  cm"1  to  become  more 
prominent.  Mendelsohn  and  co-workers  recently  showed  that  the  1690  cm"1  band  intensity  in  bone  matrix  increases 
relative  to  the  intensity  of  the  1663  cm’1  band  when  dehydrodihydroxylysinonorleucine, 

dehydrohydroxylysinonorleucine  or  dehydrohistindohydroxymerodesmosine  cross-links  are  chemically  reduced [211.  Our 
own  work  with  fetal  murine  calvarial  tissue  shows  that  the  matrix  amide  I  band  in  newly  deposited  tissue  has  a 
prominent  shoulder  at  1685  -  1690  cm'1  that  becomes  smaller  as  the  tissue  ages  and  cross-links  are  formed.  Similarly  an 
increase  in  the  1690  cm’Vl660  cm"1  band  intensity  ratio  accompanies  mechanical  rupture  of  cross-links  by  shearing. 
Although  an  exact  assignment  is  not  possible,  the  correlation  of  a  matrix  collagen  band  at  1685-1690  cm’1  and  lack  of 
mature  intrafibril  cross-links  is  clear. 


1610  1630  1650  1670  1690  1710 

Raman 


Figure  4.  Curve-fitting  of  amide  I  band,  (a),  (b)  2  month  old  wild-type  mouse  (c),  (d)  2  month  old  BRTL  mouse 
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Here  we  have  performed  an  empirical  analysis  of  immature/mature  cross-link  ratios  as  the  ratio  of  the  areas  of  two 
amide  I  sub-bands  centered  at  -1685  cm"1  and  -1667  cm’1.  The  immature/mature  cross-link  ratio  of  BRTL  mice  is  32% 
higher  than  wild-type  mice.  As  shown  in  Figure  5,  the  1685/1667  band  area  ratio  for  BRTL  mice  is  0.98510.226,  while 


Figure  5.  1685  cm"Vl667  cm  1  band  area  ratio. 

in  the  wild-type  mice  is  0.74510.227.  In  Figure  4,  BRTL  mice  show  more  distinct  1685  cm"1  bands  than  those  in  wild- 
type  mice.  As  we  have  discussed,  an  increase  in  the  intensity  of  the  1685  cm’1  band  is  a  marker  for  diminished  interfibril 
cross-links.  This  data  suggests  that  the  glycine  to  cystein  mutation  disrupts  normal  collagen  interfibril  cross  linking. 

It  is  known  that  absence  of  collagen  intrafibril  cross-links  weakens  bone  tissue1221.  In  recent  studies  of  indentation  of 
bone  tissue,  we  have  shown  that  the  intrafibril  cross-links  are  easily  ruptured  by  shearing  but  are  not  measurably 
damaged  by  indentation  to  depths  of  100  pm[23l  Using  BRTL  mice,  in  which  substitution  of  cysteine  for  glycine-349 
disrupts  the  triple  helical  structure  we  have  shown  that  presence  or  absence  of  cross-links  may  be  an  important  chemical 
composition  contributor  to  the  mechanical  strength  of  bone  tissue. 


5.  CONCLUSIONS 

We  have  shown  that  the  cortical  bone  tissue  of  two  month  old  BRTL  mice  has  a  lower  mineral  to  matrix  ratio  than  that 
of  wild-type  mice  of  the  same  age.  Equally  important,  we  have  shown  that  the  mineral  in  BRTL  differs  from  that  in  the 
wild  type  and  that  the  differences  are  characteristic  of  less  mature  and  less  homogeneous  mineralization.  Similarly,  the 
matrix  collagen  of  BRTL  differs  from  that  of  the  wild  type  mouse.  Our  BRTL  collagen  amide  I  data  is  interpreted  as 
disruption  of  cross-links.  It  is  also  possible,  for  example,  that  BRTL  matrix  collagen  has  a  higher  ratio  of  immature  to 
mature  cross-links  than  wild  type  tissue.  A  more  quantitative  description  awaits  the  results  of  more  extensive  correlation 
of  bone  matrix  amide  I  envelope  data  and  independent  cross-link  analyses.  Such  measurements  are  underway  in  our 
laboratory.  Information  contained  in  the  bone  mineral  and  the  bone  matrix  vibrational  spectra  may  be  useful  predictors 
of  stiffness  and  strength  of  the  bone  tissue.  Raman  or  infrared  spectroscopy  of  bone  tissue  may  provide  important 
information  on  the  contribution  of  bone  matrix  to  the  mechanical  properties  of  the  tissue. 
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ABSTRACT:  The  goal  of  this  study  was  to  investigate  genetic  effects  on  mechanical  properties  of  the  mouse 
femur.  We  found  evidence  for  QTL  on  eight  chromosomes  that  affect  mechanical  traits.  Some  of  these  QTL 
may  have  primary  effects  on  body  weight  or  femoral  geometry,  and  others  seem  to  affect  bone  quality  directly. 

Introduction:  Previous  studies  have  shown  a  dependence  of  fragility-related  fracture  risk  on  genetic  background. 
Although  many  of  these  studies  investigated  the  effect  of  genetics  on  BMD,  basic  measures  of  bone  geometry  and 
mechanical  integrity  may  provide  a  more  comprehensive  characterization  of  the  genetic  effects  on  bone  fragility.  The 
purpose  of  this  study  was  to  identify  quantitative  trait  loci  (QTL)  that  affect  mechanical  and  material  properties  of 
cortical  bone  in  a  genetically  heterogeneous  mouse  population. 

Materials  and  Methods:  A  total  of  486  female  UM-HET3  mice  was  used  for  this  study.  UM-HET3  mice  are 
produced  as  the  offspring  of  (BALB/cJ  X  C57BL/6J)  F,  females  and  (C3H/HeJ  X  DBA/2J)  Fj  males.  Femurs  from 
18-month-old  mice  were  tested  to  failure  in  four-point  bending  to  assess  mechanical  properties  of  cortical  bone;  these 
properties  were  compared  with  genotype  data  from  185  biallelic  loci.  A  permutation-based  test  was  used  to  detect 
significant  associations  between  genetic  markers  and  mechanical  traits.  This  test  generates  p  values  that  account  for 
the  effect  of  testing  multiple  hypotheses.  Throughout  the  experiment,  p  <  0.05  was  considered  statistically 
significant.  Analysis  of  covariance  was  used  to  examine  possible  effects  of  body  weight  and  femoral  geometry. 
Results:  We  found  evidence  for  genes  on  maternal  chromosomes  11  and  13  and  paternal  chromosomes  2,  4,  7,  10, 
1 1 ,  and  1 7  that  affect  mechanical  and  material  properties  of  femoral  bone.  The  total  variance  explained  by  genetic 
effects  on  each  mechanical  trait  ranges  from  2.9%  to  15.4%.  Most  of  the  identified  polymorphisms  influence 
mechanical  traits  even  after  adjustment  for  body  weight.  Adjustment  for  femoral  geometry  reduces  the  effects  of 
some  of  the  QTL,  but  those  on  chromosomes  2  and  10  do  not  seem  to  be  influenced  by  femoral  geometry. 
Conclusions:Many  genes  and  chromosomes  are  involved  in  the  genetic  control  over  mechanical  integrity  of  cortical 
bone.  QTL  on  paternal  chromosomes  4  and  1 1  may  mediate  mechanical  properties,  at  least  in  part,  by  modulation 
of  femoral  geometry.  Other  QTL  identified  here  may  directly  affect  bone  tissue  quality. 
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INTRODUCTION 

The  National  Osteoporosis  Foundation  (NOF)  esti¬ 
mates  that  1.5  million  fragility-related  fractures  occur 
annually  in  the  United  States.  The  direct  costs  of  these 
fractures  were  an  estimated  $17  billion  in  2001,  and  this 
figure  is  rising.  In  addition  to  the  staggering  economic 
impact,  osteoporosis  often  causes  devastation  to  the  lives  of 
its  victims  and  their  families.  Only  15%  of  hip  fracture 
patients  can  walk  across  a  room  unaided  6  months  after  a 
fracture.  An  estimated  24%  of  hip  fracture  patients  >50 
years  of  age  die  within  1  year  of  their  fracture.* 11 
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The  risk  of  fracture  depends  on  factors  such  as  the 
amount  and  structure  of  bone  that  is  present  and  the  me¬ 
chanical  integrity  of  the  bone  tissue.  These  factors  may,  in 
turn,  be  affected  by  the  genetic  background,  hormonal  mi¬ 
lieu,  age,  disease  status,  and  physical  activity  of  the  indi¬ 
vidual.  At  this  time,  the  dependence  of  fracture  risk  on 
many  of  these  factors  is  not  well  understood.  A  number  of 
previous  studies  investigated  the  genetic  control  of  fracture 
risk;  however,  most  of  these  studies  focused  on  BMD, 
because  BMD  is  the  primary  factor  used  to  predict  fracture 
risk  clinically.  Twin  and  multigeneration  family  studies 
have  estimated  that  50-90%  of  the  normal  variability  in 
BMD  at  various  anatomical  locations  is  genetically 
determined.'2-61  Inbred  mouse  strains  have  also  been  used 
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to  investigate  the  genetics  of  bone  characteristics.  Beamer  et 
al.(7)  showed  differences  in  BMD  among  11  inbred  strains 
of  mice.  Other  studies  have  similarly  reported  differences  in 
BMD  between  various  inbred  strains. *8~U)  A  few  studies 
have  reported  specific  quantitative  trait  loci  (QTL)  that 
associate  with  BMD  in  mice.112  17)  Although  the  findings  of 
these  studies  are  important,  BMD  may  not  be  the  best 
outcome  measure  to  predict  fracture  risk.  We  believe  that  a 
more  complete  understanding  of  the  genetic  control  of  frac¬ 
ture  risk  can  be  attained  by  determining  the  effects  of 
genetic  background  on  specific  measures  of  bone  size,  struc¬ 
ture,  elasticity,  strength,  and  other  indices  of  mechanical 
integrity.  This  type  of  analysis  is  obviously  not  feasible  in 
humans,  but  inbred  mouse  strains  can  be  used  to  develop 
suitable  models,  because  detailed  geometric  and  mechanical 
analysis  can  be  performed  on  their  bones.  In  addition,  the 
genetic  backgrounds  of  these  animals  are  known,  and  dif¬ 
ferences  in  bone  properties  exist  between  strains. 

In  a  previous  study,  we  identified  14  QTL  that  regulate 
various  measures  of  femoral  geometry  in  a  genetically  het¬ 
erogeneous  mouse  population/18'  However,  femoral  geom¬ 
etry  cannot  completely  predict  fracture  risk,  because  femo¬ 
ral  strength  depends  on  both  the  geometric  and  material 
properties  of  the  bone.  A  few  previous  studies  have  begun 
to  investigate  the  genetic  determinants  of  cortical  bone 
mechanical  properties  in  mice.  Differences  in  various  me¬ 
chanical  properties  have  been  shown  between  inbred 
strains. (911,19'20)  Other  studies  have  identified  specific  QTL 
that  associate  with  femoral  strength121'  and  work  to  fail¬ 
ure*22'  in  mice.  The  purpose  of  this  study  was  to  identify 
QTL  that  affect  mechanical  and  material  properties  of  cor¬ 
tical  bone  in  a  genetically  heterogeneous  mouse  population 
to  enhance  our  understanding  of  the  genetic  control  of 
cortical  bone.  We  hypothesized  that  we  would  discover 
some  QTL  whose  effects  on  mechanical  properties  of  cor¬ 
tical  bone  would  be  mediated  by  effects  on  femoral  size  and 
shape.  We  also  predicted  that  other  QTL  would  influence 
mechanical  properties  independently  of  femoral  geometry, 
thus  representing  polymorphisms  that  directly  affect  the 
quality  of  cortical  bone  tissue  in  mice. 

MATERIALS  AND  METHODS 

Mice  and  husbandry 

The  animals  used  in  this  study  are  the  same  as  those  used 
in  a  previous  report  of  genetic  effects  on  geometric  traits  of 
cortical  bone.*18’  Briefly,  the  mouse  population  is  of  the 
UM-HET3  stock  and  is  derived  from  a  four-way  breeding 
among  four  inbred  strains:  BALB/cJ  (C),  C57BL/6J  (B6), 
C3H/HeJ  (C3),  and  DBA/2J  (D2).  The  experimental  ani¬ 
mals  are  the  female  progeny  of  (C  X  B6)  F,  females  and 
(C3  X  D2)  F[  males. 

The  Fj  breeding  animals  were  purchased  from  Jackson 
Laboratories  (Bar  Harbor,  ME,  USA)  and  were  used  to 
produce  a  total  of  525  UM-HET3  mice.  Animals  were 
housed  by  sex  in  a  single  suite  of  specific-pathogen-free 
(SPF)  rooms  and  were  exposed  to  identical  environmental 
conditions  (12: 12-h  light/dark  cycle,  23°C).  Mice  were 
given  ad  libitum  access  to  water  and  laboratory  mouse 
chow.  The  cages  were  covered  with  microisolator  tops  to 


minimize  the  spread  of  infectious  agents.  Sentinel  mice 
were  tested  every  3  months  to  verify  the  pathogen-free 
status  of  the  population.  All  such  tests  were  negative 
throughout  the  course  of  the  study.  UM-HET3  mice  were 
entered  into  the  study  in  a  staggered  fashion  at  a  rate  of 
25-35  mice/month.  Animals  were  killed  at  18  months  of 
age,  at  which  time  right  femurs  were  removed,  dissected 
free  of  soft  tissue,  and  frozen  in  a  PBS  solution.  This  work 
was  approved  by  the  Animal  Care  and  Use  Committee  at  the 
University  of  Michigan. 

Genotyping  assays 

Genomic  DNA  was  prepared  from  1-cm  sections  of  tail 
from  4-week-old  animals  using  a  standard  phenol-extraction 
method.*23'  Final  DNA  preparations  were  tested  for  concen¬ 
tration,  ability  to  sustain  PCR  amplification  under  standard 
conditions,  and  electrophoretic  size  distribution.  DNA  was 
genotyped  using  an  ALFexpress  automated  sequence  ana¬ 
lyzer  (Pharmacia,  Piscataway,  NJ,  USA);  the  details  of  this 
genotyping  method  have  been  described  previously.*24' 
Primer  pairs  were  purchased  from  MWG  Biotech  (High 
Point,  NC,  USA).  In  total,  185  markers  were  examined  from 
99  genetic  loci.  Of  the  99  loci,  86  markers  were  informative 
for  both  the  maternal-  and  paternal-derived  alleles,  and  13 
loci  were  only  informative  for  either  maternal  or  paternal 
alleles.  The  selection  of  genetic  loci  is  described  previously 
in  detail.* 18'  Chromosomal  localization  and  order  of  markers 
were  calculated  using  the  MapMaker  QTX  program  pack¬ 
age  (Whitehead  Institute,  MIT,  Boston,  MA,  USA). 

Mechanical  testing 

After  the  right  femurs  were  scanned  with  p,CT  for  a 
previous  geometric  analysis,  they  were  mechanically  tested 
in  four-point  bending.  A  servohydraulic  system  (858  Mini 
Bionix  II;  MTS,  Eden  Prairie,  MN,  USA)  with  a  custom 
testing  apparatus  (Fig.  1A)  was  used  to  test  the  femurs  to 
failure  at  a  constant  displacement  rate  of  0.5  mm/s.  The 
femurs  were  loaded  in  the  anterior-posterior  direction  so 
that  the  posterior  side  of  the  bone  was  in  tension  and  the 
anterior  side  was  in  compression  (Fig.  IB).  All  four  loading 
points  were  placed  in  contact  with  the  bone  by  adjusting  the 
heights  of  the  upper  two  points  independently.  A  50-lb  load 
cell  (Model  41;  Sensotec,  Columbus,  OH,  USA)  was  used  to 
measure  the  load  applied  to  the  bone,  and  the  mid- 
diaphyseal  displacement  was  measured  with  a  linear  vari¬ 
able  differential  transducer  (010  MHR;  Schaevitz  Engineer¬ 
ing,  Pennsauken,  NJ,  USA).  Load  and  displacement  data 
were  acquired  using  the  TestStar  IIs  system  (version  2.4; 
MTS)  at  a  sampling  frequency  of  2000  Hz.  MATLAB 
(version  6.5;  The  Mathworks,  Natick,  MA,  USA)  was  used 
with  a  custom  analysis  program  to  determine  stiffness,  yield 
load  and  displacement,  ultimate  load  and  displacement, 
postyield  displacement  and  energy  to  failure  from  the  load, 
and  displacement  data.  The  analysis  program  contained  an 
algorithm  that  determined  the  linear  region  of  the  load- 
displacement  curve;  stiffness  was  measured  as  the  slope  of 
this  linear  region.  The  yield  point  was  defined  as  the  point 
at  which  the  secant  stiffness  differed  by  10%  from  the  initial 
tangential  stiffness.  Ultimate  load  was  the  maximum  load 
attained  before  failure,  and  ultimate  displacement  was  the 
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FIG.  1.  (A)  Custom-testing  appara¬ 

tus  used  for  four-point  bending  tests 
of  femurs.  (B)  Femurs  were  loaded  in 
a  posterior  to  anterior  direction  so 
that  the  posterior  surface  was  under 
tension  and  the  anterior  surface  was 
under  compression.  (C)  Spatial  pa¬ 
rameters  used  to  estimate  elastic 
modulus  in  Eq.  1. 


corresponding  displacement.  Postyield  displacement  was 
calculated  as  the  displacement  at  failure  minus  the  displace¬ 
ment  at  the  yield  point.  Energy  to  failure  was  determined 
with  numerical  integration  as  the  area  under  the  load- 
displacement  curve  up  to  the  fail  point.  A  displacement  ratio 
was  calculated  as  the  ratio  of  ultimate  displacement  to  yield 
displacement  to  characterize  the  relative  magnitudes  of  elas¬ 
tic  and  plastic  deformation.  Additionally,  elastic  modulus 
(E)  and  ultimate  strength  (cr)  were  estimated  using  these 
equations: 


S{3ab2  -  4b3) 

12/ 

(1) 

MC 

(2) 

where  S  is  the  whole  bone  stiffness  determined  with  four- 
point  bending,  I  is  the  moment  of  inertia  with  respect  to  the 
medial-lateral  axis  of  the  bone,  a  and  b  relate  to  the  spacing 
of  the  contact  points  and  are  defined  as  shown  in  Fig.  1C,  M 
is  the  moment  at  maximum  load,  and  C  is  the  distance  from 
the  centroid  of  the  bone  to  the  periosteal  surface  in  the 
posterior  direction.  I  and  C  were  determined  for  each  femur 
from  p CT  data. 

Testing  was  performed  on  6  different  days.  Each  day,  a 
set  of  assays  was  preceded  by  tests  of  pieces  of  steel  music 
wire  (type  302/304  stainless  steel  spring  wire,  0.020  in 
diameter,  McMaster-Carr  product  8908K21;  Aurora,  OH, 
USA)  to  evaluate  the  consistency  of  the  assay  system.  This 
wire  was  chosen  because  it  has  a  stiffness  within  the  same 
range  as  the  murine  femurs.  Six  wire  specimens  were  eval¬ 
uated  at  the  beginning  of  each  day’s  work  and  every  2-3  h 
thereafter.  ANOVA  showed  no  significant  differences 


across  the  series  of  assays  (p  =  0.22).  The  means  and  SD 
for  the  6  testing  days  were  169.7  ±  2.6,  169.2  ±  3.0, 
168.6  ±  2.3,  167.9  ±  3.1,  167.4  ±  1.7, and  169.7  ±  3.3, 
respectively,  with  an  average  CV  of  1.62%  (1.56%,  1.75%, 
1.38%,  1.83%,  1.00%,  and  1.97%,  respectively). 

Statistical  methods 

A  single-point,  genome-wide  search  was  performed  for 
each  trait  to  detect  QTL  that  may  be  associated  with  the 
trait.  To  make  the  analysis  consistent  for  all  partially  and 
fully  informative  markers,  four-way  informative  markers 
were  split  into  two  sets  of  biallelic  markers  that  were 
informative  for  either  the  maternally  or  paternally  transmit¬ 
ted  alleles.  One-way  ANOVA  models,  with  a  trait  as  the 
dependent  variable  and  a  biallelic  marker  as  the  factor  with 
two  levels,  were  used  to  perform  genome-wide  searches  for 
all  185  biallelic  markers.  The  strength  of  linkage  associa¬ 
tions  between  genetic  markers  and  mechanical  traits  was 
evaluated  using  a  permutation-based  test  of  statistical  sig¬ 
nificance.  This  test  generates  an  “experiment-wise”  accep¬ 
tance  criterion  to  take  into  account  the  multiple  hypotheses 
that  were  tested  in  this  search  and  to  avoid  type  I  error 
inflation.'25'  A  null  distribution  for  permutation  analysis 
was  generated  based  on  1000  shuffles  of  original  phenotype 
data.  Throughout  the  experiment,  p  <  0.05  was  considered 
statistically  significant.  The  percent  of  variance  in  each 
mechanical  trait  that  can  be  explained  by  genetic  effects  was 
estimated  in  a  standard  way  from  corresponding  regression 
models. 

The  animals  used  in  this  study  were  of  various  sizes  and 
weights.  To  account  for  overall  animal  size,  covariance 
analyses  were  performed  with  a  trait  as  the  dependent 
variable,  a  biallelic  marker  as  the  factor  with  two  levels,  and 
body  weight  at  3  or  18  months  of  age  as  the  covariate. 
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Table  1.  Summary  of  Mechanical  Traits 


Trait 

Mean  ±  SD 

Stiffness  (N/mm) 

451.3  ±  94.0 

Yield  load  (N) 

33.3  ±  7.7 

Ultimate  load  (N) 

42.6  ±  10.4 

Yield  displacement  (mm) 

0.076  ±  0.015 

Ultimate  displacement  (mm) 

0.13  ±  0.04 

Postyield  displacement  (mm) 

0.056  ±  0.041 

Displacement  ratio  (unitless) 

1.79  ±  0.66 

Energy  to  failure  (N  •  mm) 

3.56  ±  1.96 

Predicted  modulus  (GPa) 

12.5  ±  3.0 

Predicted  ultimate  strength  (MPa) 

188.6  ±  47.3 

Femurs  were  tested  to  failure  in  four-point  bending,  and  10  properties 
were  determined  for  each  femur.  The  mean  value  of  each  trait  is  shown  with 
the  SD.  Displacement  ratio  is  the  ratio  between  ultimate  displacement  and 
yield  displacement.  Predicted  modulus  and  predicted  ultimate  strength  were 
calculated  according  to  Eqs.  1  and  2,  respectively. 


Permutation-based  statistical  tests  were  used  to  evaluate 
statistical  significance  as  described  above.  Because  whole 
bone  mechanical  properties  are  dependent  on  the  geometry 
of  the  bone,  a  second  set  of  covariance  analyses  was  per¬ 
formed  to  account  for  variations  in  geometric  properties  of 
the  femur.  Femoral  cross-sectional  area,  cortical  thickness, 
and  moment  of  inertia  with  respect  to  the  medial/lateral  axis 
of  the  bone  were  used  simultaneously  as  covariates.  These 
measures  of  femoral  geometry  were  determined  with  pCT 
as  described  previously.* 18) 

RESULTS 

Of  the  525  mice  used  in  this  study,  38  died  before  18 
months  of  age  and  were  not  included  in  this  analysis. 
Mechanical  data  could  not  be  obtained  on  1  additional 
mouse;  therefore,  the  final  population  consisted  of  486 
UM-HET3  mice.  On  average,  there  were  4.5  marker  loci  per 
autosome,  and  the  average  distance  between  marker  loci 
varied  from  15  to  20  cM. 

Ten  variables  were  determined  for  each  specimen  from 
the  four-point  bending  tests.  These  variables  are  shown  in 
Table  1  with  their  means  and  SD  among  the  UM-HET3 
mice  used  in  this  study.  Table  2  shows  Spearman  rank-order 
correlations  among  these  mechanical  traits  and  between  the 
traits  and  body  weight  measured  either  at  3  or  1 8  months  of 
age.  Not  surprisingly,  many  of  these  indices  show  signifi¬ 
cant  correlations  with  one  another,  although  the  many  cor¬ 
relations  below  r  =  0.5  suggest  that  most  of  these  measure¬ 
ments  provide  information  about  bone  properties  not  fully 
captured  by  any  other  single  variable.  Body  weight  mea¬ 
sures  show  weak  correlations  with  the  bone  measurements 
other  than  predicted  modulus  and  predicted  ultimate 
strength,  which  are  calculated  using  information  about  bone 
geometry  likely  to  be  modulated  by  overall  body  size  and 
weight. 

Table  3  shows  significant  linkage  associations  between 
mechanical  traits  and  genetic  markers  that  were  determined 
from  the  permutation  analysis.  Because  the  maternal  and 
paternal  alleles  were  separated  for  the  linkage  analysis,  the 
origin  of  each  marker  is  indicated  with  an  M  or  a  P  in  Table 


3.  The  chromosomal  position  of  each  marker  is  reported  in 
centimorgans  and  millions  of  basepairs  from  the  centro¬ 
mere.  The  estimated  genetic  effect  for  each  significant  as¬ 
sociation  is  given  with  the  SE  in  parentheses.  The  estimated 
genetic  effect  is  the  difference  between  mean  levels  of  the 
indicated  trait  in  the  two  groups  of  mice  differing  at  the 
indicated  marker  locus.  Units  of  estimated  effects  are  the 
same  as  the  variable  units  shown  in  Table  1.  The  column 
labeled  “Comparison”  indicates  the  background  allele  at  a 
particular  marker  that  results  in  the  larger  value  of  the  trait. 
The  maternally  inherited  alleles  were  either  C  or  B6  alleles, 
and  the  paternally  inherited  alleles  were  either  C3  or  D2 
alleles.  The  genome  survey  thus  indicated  the  presence  of 
QTL  on  eight  chromosomes,  including  maternal  chromo¬ 
somes  11  and  13  and  paternal  chromosomes  2,  4,  7,  10,  11, 
and  17. 

The  percent  variance  estimated  in  each  mechanical  trait 
that  can  be  attributed  to  genetic  effects  was  determined,  and 
the  results  are  shown  in  Table  4.  The  variance  explained  by 
individual  chromosomes  was  determined,  considering  only 
the  significant  marker/trait  associations  shown  in  Table  3. 
The  last  row  in  Table  4  shows  the  total  variance  that  can  be 
attributed  to  all  QTL  that  were  significantly  associated  with 
each  trait.  The  total  variance  explained  by  genetic  effects 
ranges  from  2.9%  to  15.4%. 

We  conducted  analyses  of  covariance  to  determine  if  the 
QTL  of  interest  listed  in  Table  3  might  directly  affect  body 
weight,  and  therefore,  indirectly  affect  mechanical  proper¬ 
ties.  Body  weights  at  3  and  18  months  of  age  were  used  as 
covariates,  and  the  results  are  shown  in  Table  5.  The  esti¬ 
mated  genetic  effects  on  the  mechanical  traits  are  shown 
before  and  after  adjustment  for  body  weight.  Estimated 
effects  are  reported  for  all  significant  linkage  associations 
shown  in  Table  3  as  well  as  for  a  few  associations  (chro¬ 
mosomes  2,  4,  5,  7,  and  11;  Table  5,  bold)  that  obtained 
significance  only  after  adjustment  for  one  of  the  weight 
measures.  Genome  scan-seeking  QTL  that  modulate  weight, 
per  se,  found  no  loci  on  any  of  these  chromosomes  that  had 
a  significant  effect  on  body  weight  at  3  or  18  months. 

A  similar  analysis  was  conducted  to  determine  if  the 
effects  presented  in  Table  3  were  potentially  mediated  by 
the  geometry  of  the  femoral  diaphysis.  The  results  are 
shown  in  Table  6.  Linkage  associations  involving  predicted 
elastic  modulus  and  predicted  ultimate  strength  were  not 
included  in  these  analyses,  because  femoral  geometry  has 
already  been  accounted  for  in  the  calculation  of  these  traits. 
The  regression  equation  included  cross-sectional  area 
(CSA),  cortical  thickness  (CT),  and  moment  of  inertia  with 
respect  to  the  medial-lateral  axis  of  the  bone  (/ML)  as 
covariates. 

DISCUSSION 

We  have  shown  evidence  for  QTL  on  maternal  chromo¬ 
somes  11  and  13  and  paternal  chromosomes  2,  4,  7,  10,  11, 
and  17  that  affect  whole  bone  mechanical  and  material 
properties  of  the  femur.  Together,  these  QTL  account  for 
2.9-15.4%  of  the  variance  estimated  in  the  mechanical 
traits.  Paternal  chromosomes  2,  4,  7,  and  11  contained 
multiple  markers  associated  with  the  mechanical  properties. 
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Table  2.  Correlations  Among  Mechanical  Factors,  Including  Body  Weight  at  3  and  18  Months 


Yield 

load 

Ultimate 

load 

Yield  Ultimate  Postyield 

displacement  displacement  displacement 

Displacement 

ratio 

Energy 

to 

failure 

Predicted 

modulus 

Predicted 

ultimate 

strength 

Weight 
at  3 

months 

Weight 
at  18 

months 

Stiffness 

0.72 

0.69 

-0.10 

-0.01 

0.04 

0.06 

0.29 

0.52 

0.43 

0.03 

0.07 

Yield  load 

0.67 

0.54 

0.05 

-0.13 

-0.22 

0.31 

0.46 

0.49 

0.10 

0.12 

Ultimate  load 

0.12 

0.61 

0.54 

0.48 

0.83 

0.46 

0.78 

-0.04 

0.02 

Yield  displacement 

0.07 

-0.26 

-0.42 

0.08 

0.05 

0.17 

0.12 

0.10 

Ultimate 

0.93 

0.84 

0.93 

0.08 

0.59 

-0.11 

-0.09 

displacement 

Postyield 

0.98 

0.87 

0.08 

0.52 

-0.16 

-0.13 

displacement 

Displacement  ratio 

0.79 

0.07 

0.45 

-0.15 

-0.15 

Energy  to  failure 

0.26 

0.73 

-0.10 

-0.06 

Predicted  modulus 

0.78 

-0.46 

-0.34 

Predicted  ultimate 

-0.39 

-0.31 

strength 

Weight  at  3  0.69 

months 


Values  tabulated  are  Spearman  rank  correlation  coefficients,  with  pairwise  exclusion  for  missing  data. 

Bold  indicates  p  <  0.001. 

N  ~  550  for  correlations  between  mechanical  traits.  N  ~  450  for  correlations  between  mechanical  traits  and  weight  at  3  months.  N  ~  205  for 
correlations  between  mechanical  traits  and  weight  at  18  months. 


Table  3.  Associations  with  Experimentwise  Significance  of  p  <  0.05 


Trait 

Marker 

Chromosomal 
position  (cM) 

Chromosomal 
position  (Mb) 

Estimated  genetic 
effect  (SE) 

Experiment-wise 
p  value 

Comparison 

Stiffness 

DllMit2  (P) 

2.4 

12.25 

39.4  (8.5) 

0.004 

C3  >  D2 

D11MU83  (P) 

16.0 

28.85 

43.1  (8.7) 

0.002 

C3  >  D2 

Yield  load 

D4MH55  (P) 

19.8 

44.54 

2.58  (0.71) 

0.048 

C3  >  D2 

D4MH84  (P) 

37.7 

74.51 

2.64(0.71) 

0.034 

C3  >  D2 

D4MH55  (P) 

49.6 

100.28 

2.77  (0.72) 

0.025 

C3  >  D2 

D11MH2  (P) 

2.4 

12.25 

3.35  (0.70) 

<0.0001 

C3  >  D2 

D11MU83  (P) 

16.0 

28.84 

3.15  (0.72) 

0.002 

C3  >  D2 

D11MU156  (P) 

34.0 

63.96 

2.92  (0.70) 

0.004 

C3  >  D2 

D11MU289  (P) 

55.0 

95.63 

2.71  (0.71) 

0.030 

C3  >D2 

Ultimate  load 

D4MH155  (P) 

49.6 

100.28 

3.59  (0.97) 

0.043 

C3  >  D2 

Ultimate  displacement 

D10Mit40  (P) 

29.0 

48.82 

0.015  (0.004) 

0.006 

C3  >D2 

Postyield  displacement 

D10MH40  (P) 

29.0 

48.82 

0.014  (0.004) 

0.036 

C3  >D2 

Displacement  ratio 

D2MU434  (P) 

41.5 

71.66 

0.23  (0.06) 

0.033 

C3  >D2 

D2MU285  (P) 

86.0 

153.66 

0.22  (0.06) 

0.049 

C3  >D2 

Predicted  modulus 

D11MU83  (M) 

16.0 

28.85 

1.3  (0.3) 

<0.001 

C>B6 

D11MH83  (P) 

16.0 

28.85 

1.3  (0.3) 

0.002 

C3  >D2 

D11MH156  (P) 

34.0 

63.96 

1.1  (0.3) 

0.016 

C3  >D2 

D11MH289  (P) 

55.0 

95.63 

1.2  (0.3) 

<0.001 

C3  >D2 

D13MH26  (M) 

38.0 

69.31 

1.1  (0.3) 

0.019 

C>B6 

D17MU185  (P) 

40.6 

67.59 

1.0  (0.3) 

0.017 

C3  >D2 

Predicted  ultimate  strength 

D7MU25  (P) 

16.0 

27.39 

17.3  (4.3) 

0.012 

D2  >  C3 

D7MH91  (P) 

28.1 

48.22 

19.9  (4.3) 

<0.001 

D2  >  C3 

D7MH253  (P) 

52.8 

105.19 

16.5  (4.3) 

0.027 

D2  >  C3 

DllMit289  (P) 

55.0 

95.63 

21.1  (4.3) 

<0.001 

C3  >  D2 

A  permutation-based  test  was  used  to  detect  significant  associations  between  genetic  markers  and  mechanical  traits.  This  test  generates  experiment-wise 
p  values,  which  account  for  the  effect  of  testing  multiple  hypotheses.  The  origin  of  the  marker  (maternal  or  paternal)  is  indicated  in  parentheses  in  the  Marker 
column.  The  position  of  the  marker  on  the  chromosome  is  given  in  centimorgans  (cM)  and  millions  of  base  pairs  (Mb)  from  the  centromere,  as  defined  in 
the  mouse  draft  sequence  build  MGSCv3.  The  estimated  genetic  effect  is  the  difference  between  mean  levels  of  the  indicated  trait  in  the  two  groups  of  mice 
differing  at  the  indicated  marker  locus.  SE  is  shown  in  parentheses,  and  units  for  each  variable  are  as  stated  in  Table  1.  Comparison  indicates  which 
background  allele  results  in  the  larger  value  of  a  particular  trait.  B6  and  C  are  maternal  alleles,  whereas  C3  and  D2  are  paternal  alleles. 


It  should  be  noted  that  this  does  not  necessarily  prove  the  show  associations  with  traits  because  of  linkage  to  a  single 
existence  of  multiple  genes  that  affect  mechanical  traits  on  effector  gene  locus.  This  initial  QTL  analysis  does  not  allow 
these  chromosomes.  Closely  linked  SSLP  markers  may  for  high-resolution  mapping  of  the  various  markers  to  their 
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Table  4.  Variance  Caused  by  Individual  and  Combined  Genetic  Effects 


Chromosome 

Stiffness 

Yield  load 

Ultimate 

load 

Ultimate 

displacement 

Postyield 

displacement 

Displacement 

ratio 

Predicted 

modulus 

Predicted 

ultimate 

strength 

2P  [D2Mit285,  434] 

— 

— 

— 

— 

— 

4.37 

— 

— 

4P  [D4Mit55,  84,  155] 

— 

3.72 

2.90 

— 

— 

— 

— 

— 

7P  [D7Mit25,  91,  253] 

— 

— 

— 

— 

— 

— 

— 

5.14 

10P  [D10Mit40] 

— 

— 

— 

3.52 

2.86 

— 

— 

— 

11M  [DllMit83] 

— 

— 

— 

— 

— 

— 

4.72 

— 

IIP  [DllMit2,  83,  156,  289] 

5.80 

7.27 

— 

— 

— 

— 

6.85 

4.96 

13M  [D13Mit26] 

— 

— 

— 

— 

— 

— 

3.12 

— 

17P  [D17Mitl85] 

— 

— 

— 

— 

— 

— 

3.07 

— 

Total  variance  explained 

5.80 

10.20 

2.90 

3.52 

2.86 

4.37 

15.44 

10.64 

The  amount  of  variance  observed  in  each  trait  that  can  be  explained  by  genetic  effects  is  shown  here,  expressed  as  a  percentage  of  total  variance  observed 
for  each  trait.  The  total  variance  explained  by  each  chromosome  was  determined  with  regression  analysis.  The  chromosome  number  is  listed  with  an  M 
or  a  P  to  indicate  maternal  or  paternal  origin.  The  relevant  markers  are  listed  in  brackets  beneath  the  chromosome  number.  Only  those  associations  for  which 
experiment-wise  p  <  0.05,  as  listed  in  Table  3,  are  included  in  this  analysis.  The  last  row  shows  the  total  variance  that  can  be  attributed  to  all  significant 
markers  associated  with  each  trait. 


Table  5.  Effect  Sizes  Before  and  After  Adjustment  for  Body  Weight 


Estimated  genetic  effect  (SE,  p[t]) 
Chromosomal  _ 


Trait 

Marker 

position  (cM) 

Unadjusted 

3 -month  weight 

18-month  weight 

Comparison 

Stiffness 

DllMit2  (P) 

2.4 

39.4  (8.5,  0.000) 

36.1  (9.4,  0.000) 

31.9(13.6,0.000) 

C3  >  D2 

DllMit83  (P) 

16.0 

43.1  (8.7,0.000) 

40.5  (9.8,  0.000) 

37.5  (13.5,0.006) 

C3  >  D2 

Yield  load 

D4MH55  (P) 

19.8 

2.58(0.71.0.000) 

2.85  (0.78,  0.000) 

2.79(1.13,0.014) 

C3  >  D2 

D4MH84  (P) 

37.7 

2.64(0.71.0.000) 

2.98  (0.79,  0.000) 

3.28  (1.12,0.004) 

C3  >  D2 

D4MH155  (P) 

49.6 

2.77  (0.72,  0.000) 

3.00  (0.80,  0.000) 

3.04(1.18,0.011) 

C3  >  D2 

D7Mit25  (P) 

16.0 

1.56  (0.71.0.029) 

0.83  (0.80,  0.301) 

4.16(1.09,0.000) 

D2  >  C3 

DllMit2  (P) 

2.4 

3.35  (0.70,  0.000) 

3.10(0.78,0.000) 

3.65  (1.10,0.001) 

C3  >  D2 

DllMit83  (P) 

16.0 

3.15  (0.72,  0.000) 

2.63  (0.82,  0.002) 

3.90(1.09,0.000) 

C3  >  D2 

DllMitl56  (P) 

34.0 

2.92  (0.70,  0.000) 

2.72  (0.78,  0.000) 

3.61  (1.10,0.001) 

C3  >  D2 

D11MH289  (P) 

55.0 

2.71  (0.71.0.000) 

2.29  (0.80,  0.044) 

3.78  (1.12,0.000) 

C3  >  D2 

Ultimate  load 

D4Mit84  (P) 

37.7 

3.45  (0.96,  0.000) 

3.66(1.09,0.000) 

4.82(1.30,0.000) 

C3  >  D2 

D4MH155  (P) 

49.6 

3.59  (0.97,  0.000) 

4.00(1.09,0.000) 

3.98  (1.36,0.004) 

C3  >  D2 

Yield  displacement 

D7Mit25  (P) 

16.0 

0.004(0.001,  0.006) 

0.003  (0.002,  0.054) 

0.008  (0.002,  0.000) 

D2  >  C3 

Ultimate  displacement 

D10Mit40  (P) 

29.0 

0.015  (0.004,  0.000) 

0.017  (0.004,  0.000) 

0.002  (0.005,  0.771) 

C3  >  D2 

Postfield  displacement 

D10Mit40  (P) 

29.0 

0.014  (0.004,  0.000) 

0.016  (0.004,  0.000) 

0.005  (0.005,0.310) 

C3  >  D2 

Displacement  ratio 

D2MH434  (P) 

41.5 

0.23  (0.06,  0.000) 

0.17(0.07,0.011) 

0.20  (0.08,0.012) 

C3  >  D2 

D2MH285  (P) 

86.0 

0.22  (0.06,  0.000) 

0.24  (0.07,  0.000) 

0.09  (0.08,  0.269) 

C3  >  D2 

Energy  to  failure 

D2Mit58  (M) 

51.4 

0.56  (0.18,  0.002) 

0.58  (0.21,0.005) 

0.85  (0.22,  0.000) 

B6  >  C 

Predicted  modulus 

DllMit83  (M) 

16.0 

1.3  (0.3,0.000) 

1.3  (0.3,0.000) 

1.6  (0.4,  0.000) 

C>B6 

DllMit83  (P) 

16.0 

1.3  (0.3,0.000) 

1.3  (0.3,0.000) 

1.4  (0.4,  0.000) 

C3  >  D2 

DllMitl56  (P) 

34.0 

1.1  (0.3,0.000) 

1.2  (0.3,0.000) 

1.6  (0.4,  0.000) 

C3  >  D2 

DllMit289  (P) 

55.0 

1.2  (0.3,0.000) 

1.3  (0.3,0.000) 

1.8  (0.4,  0.000) 

C3  >  D2 

D13MH26  (M) 

38.0 

1.1  (0.3,0.000) 

0.9  (0.3,0.001) 

1.2  (0.4,  0.007) 

C  >  B6 

Predicted  ultimate 

D17Mitl85  (P) 

40.6 

1.0  (0.3,0.000) 

0.7  (0.3,0.011) 

0.9  (0.4,  0.039) 

C3  >  D2 

strength 

D5MU292  (P) 

80.0 

11.4(4.6,  0.011) 

17.4  (4.6,  0.000) 

4.4  (6.5,  0.497) 

D2  >  C3 

D7Mit25  (P) 

16.0 

17.3  (4.3,  0.000) 

11.1  (4.5,  0.014) 

21.8  (6.0,  0.000) 

D2  >  C3 

D7Mit91  (P) 

28.1 

19.9  (4.3,  0.000) 

12.2  (4.5,  0.008) 

20.1  (6.1,0.001) 

D2  >  C3 

D7MH253  (P) 

52.8 

16.5  (4.3,  0.000) 

13.1  (4.5,0.004) 

12.7  (6.1,0.039) 

D2  >  C3 

DllMitl56  (P) 

34.0 

14.6(4.3,0.001) 

16.4  (4.4,  0.000) 

23.0  (6.0,  0.000) 

C3  >  D2 

DllMit289  (P) 

55.0 

21.1  (4.3,0.000) 

21.6  (4.4,  0.000) 

27.7  (6.0,  0.000) 

C3  >  D2 

Two  body  weight  measures,  obtained  at  3  and  18  months  of  age,  were  used  as  predictive  covariates  to  determine  if  differences  seen  in  mechanical  traits 
were  actually  caused  by  differences  in  body  weight.  The  estimated  genetic  effects  on  the  mechanical  traits  are  shown  before  and  after  adjustment  for  body 
weight.  The  estimated  genetic  effect  is  the  difference  between  mean  levels  of  the  indicated  trait  in  the  two  groups  of  mice  differing  at  the  indicated  marker 
locus.  SE  is  shown  in  parentheses,  along  with  probability  evaluated  by  the  /-statistic,  and  units  for  each  variable  are  as  stated  in  Table  1.  Estimated  effects 
are  reported  for  all  significant  associations  shown  in  Table  2,  as  well  as  a  few  associations  (shown  in  bold)  that  obtained  significance  only  after  adjustment 
for  one  of  the  weight  measures.  Comparison  indicates  which  background  allele  results  in  the  larger  value  of  a  particular  trait.  B6  and  C  are  maternal  alleles, 
whereas  C3  and  D2  are  paternal  alleles. 
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Table  6.  Effect  Sizes  Before  and  After  Adjustment  for  Geometry 


Trait 

Marker 

Positon  (cM) 

Estimated  genetic  effect  (SE) 

Unadjusted  Adjusted 

Comparison 

Stiffness 

DlNds2  (M) 

59.0 

19.6(8.7) 

27.2(7.1) 

B6  >  C 

DllMit2  (P) 

2.4 

39.4  (8.5) 

23.9  (7.0) 

C3  >  D2 

DllMit83  (P) 

16.0 

43.1  (8.7) 

25.3  (7.3) 

C3  >  D2 

Yield  load 

D4Mit55  (P) 

19.8 

2.58  (0.71) 

1.34  (0.55) 

C3  >  D2 

D4Mit84  (P) 

37.7 

2.64  (0.71) 

1.19(0.54) 

C3  >  D2 

D4Mitl55  (P) 

49.6 

2.77  (0.72) 

1.20(0.58) 

C3  >  D2 

DllMit2  (P) 

2.4 

3.35  (0.70) 

1.83  (0.55) 

C3  >  D2 

DllMit83  (P) 

16.0 

3.15  (0.72) 

1.26  (0.57) 

C3  >  D2 

DllMitl56  (P) 

34.0 

2.92  (0.70) 

0.58  (0.56) 

C3  >  D2 

DllMit289  (P) 

55.0 

2.71  (0.71) 

0.77  (0.57) 

C3  >  D2 

Ultimate  load 

D4Mitl55  (P) 

49.6 

3.59  (0.97) 

1.48  (0.75) 

C3  >  D2 

Ultimate  displacement 

D10Mit40  (P) 

29.0 

0.015  (0.004) 

0.015  (0.004) 

C3  >  D2 

Postyield  displacement 

D10Mit40  (P) 

29.0 

0.014  (0.004) 

0.014(0.004) 

C3  >  D2 

Displacement  ratio 

D2Mit434  (P) 

41.5 

0.23  (0.06) 

0.23  (0.06) 

C3  >  D2 

D2Mit285  (P) 

86.0 

0.22  (0.06) 

0.23  (0.06) 

C3  >  D2 

Three  measures  of  geometry,  cross-sectional  area  (CSA),  cortical  thickness  (CT),  and  moment  of  inertia  with  respect  to  the  medial  lateral  axis  of  the  bone 
(/ML),  were  used  simultaneously  as  predictive  covariates  to  determine  if  differences  seen  in  mechanical  traits  were  actually  caused  by  differences  in 
geometry.  Predicted  modulus  and  predicted  ultimate  strength  were  not  included  in  this  analysis  because  these  variables  have  already  been  adjusted  for  by 
geometry.  Genetic  effects  on  the  mechanical  traits  are  shown  before  and  after  adjustment  for  femoral  geometry.  The  estimated  genetic  effect  is  the  difference 
between  mean  levels  of  the  indicated  trait  in  the  two  groups  of  mice  differing  at  the  indicated  marker  locus.  SE  is  shown  in  parentheses,  and  units  for  each 
variable  are  as  stated  in  Table  1.  Genetic  effects  are  reported  for  all  significant  associations  shown  in  Table  3,  as  well  as  one  association  (shown  in  bold) 
that  obtained  significance  only  after  adjustment  for  geometry. 


chromosomal  positions.  Therefore,  multiple  markers  on  a 
chromosome  may  be  linked  to  a  single  gene  of  interest. 
Additional  genotyping  using  more  closely  spaced  markers  is 
now  underway  and  will  help  to  determine  the  chromosomes 
that  contain  more  than  one  locus  with  effects  on  the  me¬ 
chanical  properties  of  the  femur.  Although  we  have  used  a 
permutation-based  significance  test  to  control  type  I  error 
rate  and  produce  an  “experiment-wise”  p  value  for  each 
genome  scan,  we  have  not  adjusted  the  significance  criterion 
with  respect  to  the  number  of  traits  evaluated;  thus.  Table  3 
may  well  contain  a  small  proportion  of  false  positive  con¬ 
clusions.  In  addition,  we  note  that  a  study  of  500-600  mice 
has  only  limited  statistical  power.  Power  analysis  using  a 
synthetic  data  set  suggests  that  a  population  of  this  size  will 
detect  only  27%  of  the  QTL  with  heritability  of  0.10  and 
only  45%  of  QTL  with  heritability  0.15.  Thus,  there  is  a 
good  likelihood  that  our  compilation  (Table  3)  has  missed 
many  QTL  with  substantial  effects  on  the  bone  traits  we 
have  measured. 

The  mechanical  properties  studied  here  (stiffness,  yield 
load  and  displacement,  ultimate  load  and  displacement, 
postyield  displacement,  displacement  ratio,  and  energy  to 
failure)  depend  on  material  properties  of  the  bone  tissue 
(such  as  modulus  and  ultimate  strength  that  we  have  esti¬ 
mated  in  this  study)  and  the  size  and  cross-sectional  shape 
of  the  bone.  Our  previous  geometric  analysis  revealed  QTL 
that  strongly  associate  with  various  measures  of  femoral 
size  and  shape  on  maternal  chromosomes  1,  5,  6,  7,  11,  12, 
15,  and  17  and  paternal  chromosomes  1,  3,  4,  8,  9,  11,  14, 
and  15.' 1181  A  direct  comparison  of  these  results  with  the 
current  mechanical  trait  results  (not  including  material  prop¬ 
erties  like  modulus  and  ultimate  strength)  reveals  two  areas 
of  overlap:  paternal  chromosomes  4  and  11  have  QTL  that 


associate  with  both  femoral  geometry  and  mechanical  prop¬ 
erties.  It  is  possible  that  one  or  more  genes  on  these  chro¬ 
mosomes  actually  affect  the  geometry  of  the  femur,  which 
in  turn  affects  mechanical  properties.  Because  the  remaining 
QTL  that  associate  with  mechanical  properties  are  not 
strongly  associated  with  geometric  properties,  they  may 
instead  directly  affect  bone  tissue  quality  or  they  may  affect 
bones  indirectly  (e.g.,  through  the  alteration  of  a  hormone 
that  in  turn  influences  bone  quality). 

Material  properties  of  femoral  diaphyseal  cortical  bone  in 
the  UM-HET3  mice  were  estimated  by  calculating  elastic 
modulus  and  ultimate  strength  according  to  Eqs.  1  and  2, 
respectively.  These  estimated  material  properties  were  sig¬ 
nificantly  associated  with  markers  on  maternal  chromo¬ 
somes  11  and  13  and  paternal  chromosomes  7,  11,  and  17, 
indicating  that  genes  on  these  chromosomes  affect  the  ma¬ 
terial  properties  of  cortical  bone,  independent  of  the  femoral 
geometry.  However,  it  is  noteworthy  that  one  of  these 
markers,  DllMit83  on  the  maternal  chromosome  1 1,  is  also 
strongly  associated  with  various  femoral  geometry  traits. 
True  material  properties  should  be  independent  of  geome¬ 
try.  It  is  possible  that  the  elastic  modulus  and  ultimate 
strength  values  used  here  may  not  completely  account  for 
femoral  geometry,  because  they  are  estimates  based  on 
beam  theory.  Another  possibility  is  that  a  gene  near  the 
DllMit83  marker  affects  some  other  factor,  such  as  calcium 
metabolism  or  osteoblast  regulation,  that  could  indirectly 
affect  both  geometry  and  material  properties.  Future  plans 
include  producing  uniform  microbeams  from  femoral  cor¬ 
tical  tissue  of  the  same  animals  and  testing  the  beams  in 
micro-four-point  bending  to  determine  the  material  proper¬ 
ties  of  femoral  cortical  tissue  more  directly. 
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There  are  a  relatively  small  number  of  previous  studies 
that  have  identified  QTL  and  associate  it  with  femoral 
mechanical  properties  and  BMD  in  mice.  In  1999,  Beamer 
et  al.(13)  identified  four  QTL  on  chromosomes  1,  5,  13,  and 
15  that  were  strongly  linked  with  femoral  BMD  in  an  F2 
cross  between  C57BL/6J  and  CAST/EiJ  strains.  Using  a 
different  mouse  model  (F2  cross  between  C57BL/6J  and 
C3H/HeJ),  Beamer  et  al.<15)  performed  a  similar  study  in 
2001  and  found  QTL  on  chromosomes  1,  4,  6,  11,  13,  14, 
and  18  that  strongly  associate  with  femoral  BMD.  Klein  et 
al.n6)  measured  whole  body  BMD  in  mice  derived  from 
C57BL/6  and  DBA/2  strains  in  2001.  The  results  showed 
links  with  QTL  on  chromosomes  1,  2,  4,  and  ll.(16)  Com¬ 
paring  our  current  results  with  these  previous  BMD  studies, 
we  find  commonalities  on  chromosomes  2,  4,  11,  and  13. 
Inconsistencies  among  other  significant  QTL  may  reflect 
differences  in  outcome  measures  or  in  the  specific  genetic 
crosses  used;  these  inconsistencies  could  also  reflect  type  II 
errors  in  one  or  more  of  these  studies. 

Li  et  al.<21'22)  used  an  F2  cross  between  MRL/MpJ  and 
SJL/J  strains  to  identify  QTL-affecting  femoral  failure  load 
and  energy  to  failure.  The  results  identified  markers  on 
chromosomes  1.  2,  8,  9,  10,  and  17  that  associated  with 
failure  load  and  markers  on  chromosomes  2,  7,  8,  9,  and  X 
that  associated  with  energy  to  failure.  We  found  no  QTL 
that  strongly  associated  with  energy  to  failure  and  only  one 
on  chromosome  4  that  associated  with  ultimate  load;  this  is 
similar  to  failure  load  in  the  UM-HET3  mice.  The  lack  of 
consistency  between  our  results  and  the  findings  of  Li  et 
al  (2i,22)  ^  most  likely  attributable  the  different  mouse 
models  used  in  the  studies. 

We  considered  the  possibility  that  some  of  the  genetic 
effects  listed  in  Table  3  might  act  primarily  on  body  weight, 
with  only  secondary  effects  on  mechanical  properties.  To 
test  this  idea,  we  conducted  analyses  of  covariance  using 
body  weight  at  either  3  or  18  months  of  age  as  the  covari¬ 
ates.  The  results  shown  in  Table  5  show  that,  in  most  cases, 
the  genetic  effects  have  similar  magnitudes  before  and  after 
adjustment  for  body  weight.  However,  there  are  a  few 
exceptions.  First,  the  associations  between  the  paternal  al¬ 
leles  of  D10Mit40  and  ultimate  displacement  and  postyield 
displacement  were  reduced  after  adjustment  for  18-month 
body  weight.  It  is  possible  that  the  linkage  association  of 
this  QTL  with  displacement  measures  is  mediated  indirectly 
through  an  effect  on  body  weight  throughout  midlife.  It  is 
also  possible  that  some  third  factor,  perhaps  hormones, 
modulates  both  bone  properties  and  body  weight  and  is 
under  control  of  a  locus  on  paternal  chromosome  10.  Sim¬ 
ilarly,  the  association  between  displacement  ratio  and  the 
paternal  allele  at  D2Mit285  is  diminished  in  calculations 
involving  18-month  body  weight.  It  is  interesting  to  note 
that  the  relationship  between  D2Mit434  and  displacement 
ratio  is  not  similarly  affected  by  18-month  weight.  This 
difference  may  suggest  that  chromosome  2  harbors  two 
different  loci  with  effects  on  displacement  ratio — one  that  is 
associated  with  body  weight  dependence,  and  one  that  is 
not.  Further  studies  involving  a  higher  density  of  markers  on 
chromosome  2  are  required  to  confirm  this. 

Table  5  displays  six  linkage  associations  (bold)  in  which 
QTL  had  significant  increases  in  effect  size  only  after  ad¬ 


justment  for  body  weight.  Three  of  these  associations  in¬ 
volve  D4Mit84,  D7Mit25,  and  D1  lMitl56,  which  were 
already  noted  in  Table  3  because  of  their  associations  with 
other  mechanical  traits.  Two  other  loci,  however,  are  new. 
The  estimated  genetic  effect  of  the  gene  associated  with 
D2Mit58  on  energy  to  failure  increased  only  after  adjust¬ 
ment  for  body  weight  at  1 8  months.  It  is  possible  that  this 
maternal  QTL  modulates  the  process  by  which  bone  prop¬ 
erties  adjust  to  differences  in  body  weight  or  to  physiolog¬ 
ical  factors,  perhaps  hormonal,  that  modulate  both  body 
weight  and  bone  fragility.  The  other  new  QTL  modulates 
predicted  ultimate  strength  of  the  bone  and  is  associated 
with  the  paternal  alleles  at  D5Mit292.  Here,  the  effect  size 
increases  after  adjustment  for  weight  at  3  months.  However, 
the  effect  size  increase  is  much  less  dramatic  after  adjust¬ 
ment  for  weight  at  18  months. 

A  second  set  of  covariance  analyses  were  conducted  to 
determine  if  effects  noted  in  Table  3  were  potentially  me¬ 
diated  by  genetic  influences  on  the  geometry  of  the  femoral 
diaphysis.  These  results  are  presented  in  Table  6.  For  asso¬ 
ciations  involving  markers  on  paternal  chromosomes  4  and 
11,  the  estimated  genetic  effect  is  reduced  after  adjustment 
for  femoral  geometry.  These  results  are  consistent  with  the 
fact  that  paternal  chromosomes  4  and  1 1  contained  markers 
that  were  significantly  associated  with  both  femoral  geom¬ 
etry  and  mechanical  properties.  These  genetic  markers  may 
represent  QTL  whose  effects  on  mechanical  properties  are 
mediated,  at  least  in  part,  by  modulation  of  femoral  geom¬ 
etry.  However,  the  associations  shown  in  Table  6,  involving 
markers  on  paternal  chromosomes  2  and  10,  do  not  seem  to 
be  affected  by  femoral  geometry  and  presumably  reflect 
modulation  of  bone  quality.  Only  one  new  locus,  DlNds2, 
was  found  to  have  a  stronger  association  after  adjustment 
for  geometry.  In  this  case,  mice  inheriting  the  B6  allele  were 
found  to  show  greater  stiffness  in  calculations  adjusted  for 
bone  geometrical  traits,  suggesting  that  this  QTL  might 
modulate  compensation  of  bone  mechanical  properties  to 
bone  size  or  shape. 

In  conclusion,  we  have  shown  the  complex  nature  of  the 
genetic  control  of  cortical  bone  in  mice  that  is  caused  by  the 
multiple  QTL  that  were  found  to  associate  with  various 
measures  of  femoral  mechanical  and  material  properties. 
We  found  evidence  for  QTL  on  paternal  chromosomes  4 
and  1 1  whose  effects  on  mechanical  properties  may  be 
mediated,  at  least  in  part,  by  modulation  of  femoral  geom¬ 
etry.  Other  genetic  markers  identified  in  this  survey  may  be 
linked  to  QTL  that  directly  affect  bone  tissue  quality  and 
material  properties,  but  further  studies  are  needed  to  clarify 
the  role  of  these  QTL  in  the  modulation  of  femoral  mechan¬ 
ical  integrity.  Although  the  results  of  this  study  do  not  allow 
us  to  positively  identify  candidate  genes  that  affect  femoral 
mechanical  properties  in  mice,  we  now  have  an  appropriate 
starting  point  from  which  we  can  begin  the  search  to  find 
specific  genes.  In  addition,  further  analysis  of  this  four-way 
cross-population  should  help  to  delineate  interactions 
among  bone  traits,  underlying  hormonal  and  biochemical 
differences  among  mice,  and  the  genes  that  modulate  bones 
both  during  and  after  developmental  maturation.  Finally,  the 
results  of  this  study  may  have  important  implications  with 
respect  to  fragility  or  fracture  risk  in  humans.  The  data  in 
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this  study  relate  to  10-15  QTL,  each  of  which  corresponds 
to  one  or  two  human  syntenic  regions  with  200-500  known 
genes.  Testing  or  review  of  these  genes  is  beyond  the  scope 
of  this  study,  but  focused  searches  for  human  alleles  at  the 
same  (syntenic)  positions  of  the  mouse  genes  may  provide 
additional  insight  (beyond  BMD)  into  fracture  risk  predic¬ 
tion.  In  addition,  identification  of  the  true  effector  genes  and 
their  function  in  regulating  bone  extracellular  matrix  prop¬ 
erties  in  the  mouse  may  support  the  development  of  new 
therapeutic  strategies  for  treating  or  preventing  fragility  in 
humans. 
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ABSTRACT:  The  Brtl  mouse  model  for  type  IV  osteogenesis  imperfecta  improves  its  whole  bone  strength  and 
stiffness  between  2  and  6  months  of  age.  This  adaptation  is  accomplished  without  a  corresponding  improve¬ 
ment  in  geometric  resistance  to  bending,  suggesting  an  improvement  in  matrix  material  properties. 

Introduction:  The  Brittle  IV  (Brtl)  mouse  was  developed  as  a  knock-in  model  for  osteogenesis  imperfecta  (OI)  type 
IV.  A  Gly349Cys  substitution  was  introduced  into  one  collal  allele,  resulting  in  a  phenotype  representative  of  the 
disease.  In  this  study,  we  investigate  the  effect  of  the  Brtl  mutation  on  whole  bone  architecture,  strength,  and 
composition  across  a  range  of  age  groups. 

Materials  and  Methods:  One-,  2-,  6-,  and  12-month-old  Brtl  and  wildtype  (WT)  mice  were  analyzed.  Femurs  were 
assessed  at  the  central  diaphysis  for  cortical  geometric  parameters  using  p,CT  and  were  subsequently  mechanically 
tested  to  failure  by  four-point  bending.  Matrix  material  properties  were  predicted  using  p.C'T  data  to  normalize  data 
from  mechanical  tests.  Raman  spectroscopy  and  DXA  were  used  to  assess  matrix  composition. 

Results:  Our  findings  show  a  postpubertal  adaptation  in  which  Brtl  femoral  strength  and  stiffness  increase  through 
a  mechanism  independent  of  changes  in  whole  bone  geometry.  These  findings  suggest  an  improvement  in  the 
material  properties  of  the  bone  matrix  itself,  rather  than  improvements  in  whole  bone  geometry,  as  seen  in  previous 
mouse  models  of  OI.  Raman  spectroscopic  results  suggest  these  findings  may  be  caused  by  changes  in  mineral/matrix 
balance  rather  than  improvements  in  mineral  crystallinity. 

Conclusions:  Our  findings  parallel  the  currently  unexplained  clinical  observation  of  decreased  fractures  in  human  OI 
patients  after  puberty.  The  Brtl  mouse  remains  an  important  tool  for  investigating  therapeutic  interventions  for  OI. 
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INTRODUCTION 

Osteogenesis  imperfecta  (OI)  is  an  autosomal  dominant 
genetic  disorder  of  the  extracellular  matrix.’ 1,2)  The 
hallmark  feature  of  OI  is  bone  fragility  and  susceptibility  to 
fracture  from  minimal  trauma.  Patients  with  OI  are  custom¬ 
arily  grouped  in  types  I  to  IV,  according  to  the  Sillence 
classification,  based  on  clinical  and  radiographic  features. (1) 
About  90%  of  OI  cases  are  caused  by  mutations  in  type  I 
collagen,  the  major  structural  component  of  the  extracellu¬ 
lar  matrix  of  bone  and  skin.  The  severity  of  the  disorder  is 
related  to  the  location  and  type  of  mutation  in  each  collagen 
alpha  chain.  Recently,  a  small  percentage  of  OI  patients 
(about  10%)  have  been  reclassified  as  types  V-VII  based  on 
skeletal  histology  and  radiographic  features;  these  cases  do 
not  result  from  collagen  mutations. t3~6)  Type  IV  OI  is  the 
moderately  severe  form  of  the  disorder,  with  phenotypes 
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ranging  between  the  milder  phenotype  of  OI  type  I  and  the 
more  severe  complications  of  OI  type  III.’ 11  Infants  are  often 
bom  with  fractures,  but  fracture  rate  can  be  quite  variable 
throughout  childhood.’ 11  Interestingly,  patients  with  type  IV 
OI  show  a  marked  decrease  in  fracture  rate  near  the  age  of 
puberty ,(1,7)  suggesting  the  existence  of  an  unidentified  ad¬ 
aptation  mechanism  which  increases  whole  bone  structural 
integrity. 

Although  rearrangements  and  exon  skipping  mutations 
have  been  described,  OI  type  IV  typically  (85%)  results 
from  a  glycine  substitution  in  one  of  the  collagen  I  chains. (1) 
Understanding  how  a  single  point  mutation  in  the  collagen 
gene  can  alter  the  susceptibility  of  an  entire  bone  to  fracture 
requires  observations  over  multiple  hierarchical  levels  of 
the  tissue.  Furthermore,  the  development  of  effective  ther¬ 
apeutic  agents  for  OI  may  require  manipulation  at  one  or 
more  of  these  hierarchical  levels.  Thus,  a  critical  examina¬ 
tion  of  an  appropriate  and  reproducible  model  for  OI  is 
required. 


’Orthopaedic  Research  Laboratories,  Department  of  Orthopaedic  Surgery,  University  of  Michigan,  Ann  Arbor,  Michigan,  USA; 
department  of  Chemistry,  University  of  Michigan,  Ann  Arbor,  Michigan,  USA;  'Bone  and  Extracellular  Matrix  Branch,  National 
Institute  of  Child  Health  and  Human  Development,  Bethesda,  Maryland,  USA. 
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FIG.  1.  (A)  ju-CT  image  of  central  diaphyseal 

region  of  analysis.  (B)  Femurs  were  subjected  to 
four-point  bending  with  the  posterior  surface 
under  tension. 


Several  murine  models  of  OI  have  been  reported,*8  1,1 
however,  none  of  these  mice  result  from  dominant  negative 
transmission  of  a  typical  glycine  substitution  mutation. 
Only  naturally  occurring  and  knock-in  mutations  will  have 
both  the  physiological  tissue  distribution  of  the  mutant 
collagen  transcripts  and  the  physiological  levels  of  mutant 
alpha  chain,  because  their  expression  is  under  the  control  of 
endogenous  promoters.  For  these  reasons,  the  Brittle  IV,  or 
Brtl  mouse,  was  developed  as  the  first  knock-in  murine 
model  for  OI.' 121  Brtl  is  heterozygous  for  a  Gly349— >Cys 
(G349C)  substitution  in  one  collal  allele,  reproducing  the 
genetic  mutation  found  in  a  type  IV  child.* 131  Brtl  mice 
show  a  moderately  severe  skeletal  phenotype,  with  long 
bone  deformities  and  bone  fragility.  Furthermore,  the  Brtl 
mouse  reproduces  the  molecular,  biochemical,  and  radio- 
graphic  features  of  OI  type  IV,(  121  making  it  an  excellent 
model  for  studying  this  disease. 

The  purpose  of  this  investigation  was  to  evaluate  the 
effects  of  a  single  glycine  substitution  in  the  alpha  1(1) 
collagen  chain  on  the  whole  bone  properties  of  the  Brtl 
mouse  model  for  OI.  This  involves  examination  of  the 
inherent  material  properties  of  the  bone  matrix,  the  relative 
levels  of  mineral  and  extracellular  matrix  constituencies, 
and  the  structural  organization  of  the  femur.  In  addition, 
these  studies  were  performed  in  mice  over  a  range  of  ages, 
from  skeletally  immature  to  mature  adult,  to  observe  any 
changes  in  parameters  over  time. 

MATERIALS  AND  METHODS 

Specimens 

Seventy  nine  male  Brtl  and  WT  mice  were  killed  at  1 
(Brtl,  n  =  10;  WT,  n  =  14),  2  (Brtl,  n  =  10;  WT,  n  =  12), 
6  (Brtl,  n  =  11;  WT,  n  =  12),  and  12  (Brtl,  n  =  5;  WT,  n  = 
5)  months  of  age,  corresponding  to  developmental  status  of 
childhood,  puberty,  mature  young  adult,  and  adult.  Brtl 
mice  have  a  mixed  background  of  Svl29/CD-1/C57BL/6S 
and  are  bred  by  crossing  heterozygous  Brtl  with  WT,  each 
generated  through  separate  Brtl  X  WT  matings.*12’  After 
death,  femurs  were  harvested,  stripped  of  soft  tissue,  and 
stored  frozen  in  saline  until  ready  for  testing.  In  one-half  of 


the  specimens,  left  femurs  were  analyzed  for  geometric 
properties,  and  right  femurs  were  mechanically  tested.  In 
the  remainder  of  the  population,  bone  geometry  and  me¬ 
chanical  properties  were  both  measured  in  left  femurs.  In  a 
subpopulation  of  1 1  mice,  both  femurs  were  evaluated  to 
assess  left-right  differences  in  structural  properties.  A  sep¬ 
arate  group  of  43  Brtl  and  WT  mice  were  used  to  determine 
BMD  of  the  left  femur.  All  experiments  were  performed  in 
accordance  with  approval  from  the  NICHD  ACUC  commit¬ 
tee. 

Evaluation  of  femoral  geometry  by  pCT 

A  custom  i±CT  scanner  developed  and  validated  in  our 
laboratory*  14~16’  was  used  to  scan  femurs  at  25-mm  voxel 
resolution  to  assess  geometric  properties  of  the  central  di- 
aphysis  (Fig.  1A).  The  3D  data  set  generated  by  pCT  is 
organized  as  a  series  of  slices  25  jum  thick,  oriented  along 
the  long  axis  of  the  bone.  A  threshold  was  applied  to  each 
slice,*15’  and  transverse  cross-sections  were  analyzed  over  2 
(1  month)  or  3  mm  (2,  6,  and  12  months)  segments  starting 
at  the  distal  insertion  of  the  third  trochanter  and  extending 
toward  the  distal  end  of  the  bone.  Each  slice  was  analyzed 
for  its  cross-sectional  area,  cortical  thickness,  and  bending 
moment  of  inertia  in  the  anterior-posterior  direction  about 
the  medial-lateral  axis.  Data  within  each  femur  was  aver¬ 
aged  over  the  2-  or  3 -mm  segment  to  provide  representative 
values  over  the  entire  3D  mid-diaphysis. 

Whole  bone  mechanical  properties 

Femurs  were  loaded  to  failure  in  four-point  bending  at  0.5 
mm/s  in  the  anterior-posterior  direction  using  a  servohy- 
draulic  testing  machine  (810  Material  Test  System;  MTS 
Systems  Corporation,  Eden  Prairie,  MN,  USA;  Fig.  IB). 
Regions  loaded  in  four-point  bending  corresponded  to  those 
measured  by  juCT.  Load  head  displacement  was  monitored 
using  an  external  linear  variable  differential  transducer 
(LVDT;  Lucas  Schavitts,  Hampton,  VA,  USA),  and  load 
data  were  collected  with  a  50-lb  load  cell  (Sensotec,  Co¬ 
lumbus,  OH,  USA).  Data  were  sampled  at  200  Hz  on  a 
Macintosh  Ilci  computer  through  an  A/D  system  (Lab view; 
National  Instruments,  Austin,  TX,  USA).  Load- 
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displacement  curves  were  analyzed  for  whole  bone  yield 
load,  stiffness,  yield  energy,  maximum  load,  failure  energy, 
and  post-yield  displacement. 

Predicted  matrix  properties 

Results  obtained  from  whole  bone  mechanical  tests  are 
influenced  by  the  material  properties  of  the  bone  matrix,  as 
well  as  by  the  structural  organization  of  the  whole  bone. 
Therefore,  to  separate  the  relative  contribution  of  matrix 
material  properties  from  whole  bone  structural  organization, 
material  properties  of  the  extracellular  matrix  were  pre¬ 
dicted  using  beam  theory' 171  based  on  the  geometric  prop¬ 
erties  of  the  femur  measured  by  pCT.  Predicted  ultimate 
stress  is  assumed  to  be  greatest  at  the  posterior  surface, 
which  is  the  area  under  greatest  tension  during  bending 
experiments,  and  is  calculated  as 

Fac 


whereas  the  equation  for  predicted  elastic  modulus  is 

a2 

E  =  5  y^(3 L  ~  4 a)  (2) 

where  F  equals  the  measured  maximum  load,  c  is  the 
distance  from  the  centroid  of  the  bone  to  the  posterior 
periosteal  surface,  I  is  the  moment  of  inertia  of  bending  in 
the  anterior-posterior  direction,  S  is  the  initial  stiffness  of 
the  bone,  and  L  and  a  are  geometric  constants  based  on  the 
four-point  bending  setup  (3.57  and  0.99  mm). 

These  equations  predict  inherent  material  properties  of 
the  extracellular  matrix,  assuming  a  bending  mode  of  fail¬ 
ure. 

Raman  spectroscopy 

Vibrational  spectroscopy  is  ideal  for  probing  the  material 
properties  of  biological  specimens  nondestructive^.'18’ 
Both  Raman  and  infrared  spectroscopy  yield  simultaneous 
information  about  both  the  inorganic  and  organic  compo¬ 
nents  of  bone,  yet  leave  the  specimen  intact  for  further 
testing.  Because  it  is  a  light  scattering  measurement,  Raman 
spectroscopy  is  especially  advantageous  for  the  thick  spec¬ 
imens  used  in  this  study.  An  epi-illumination  excitation/ 
collection  optical  configuration  is  used,  and  specimens  do 
not  have  to  be  transparent.  As  a  preliminary  investigation 
into  the  collagen-mineral  ultrastructure  of  the  Brtl  bone 
matrix,  12  specimens  from  2-  and  6-month-old  animals  ( n  = 
3  per  group)  were  prepared  for  Raman  spectroscopic  imag¬ 
ing  from  contralateral  femurs  by  using  a  custom  microm¬ 
illing  device  to  create  barbell- shaped  specimens. <19-21) 
Gauge  sections  were  ~200  X  300  pm  in  cross-section  and 
1000  pm  long. 

The  Raman  imaging  system  has  been  previously 
described. 122-25 1  Raman  scatter  was  excited  using  a  high- 
aspect  ratio  rectangular  laser  beam  at  785  nm  (Kaiser  Op¬ 
tical  Systems).  The  total  laser  power  at  the  specimen  surface 
was  — 150  mW.  The  beam  was  focused  onto  the  center  of 
the  specimen’s  gauge  length  through  a  10X/0.5  NA  objec¬ 


FIG.  2.  Raman  spectrum  of  a  6-month-old  WT  mouse  femur.  Impor¬ 
tant  bands  discussed  in  this  paper  are  indicated  together  with  their 
vibrational  assignments.  The  amide  II,  CH2  wag,  and  amide  I  bands 
result  from  vibrations  of  the  collagenous  matrix,  whereas  the  phosphate 
i>i  and  carbonate  bands  arise  from  the  bone  mineral. 


tive  (Zeiss),  and  the  resulting  Raman  scatter  was  collected 
in  an  epi-illumination  mode  through  the  same  objective.  In 
this  manner,  Raman  spectra  were  excited  at  every  point 
along  the  ~200-p,m-long  laser  rectangle.  The  Raman  scat¬ 
ter  was  focused  into  an  axial-transmissive  spectrograph 
(Kaiser  Optical  Systems)  and  detected  using  a  thermoelec- 
trically  cooled  CCD  detector  (Andor  Technologies).  To 
increase  the  number  of  spectra  available  for  analysis,  the 
specimen  was  translated  in  l- pm  steps  using  a  computer- 
controlled  X-Y  translation  stage  (New  England  Affiliated 
Technologies),  and  a  new  line  of  spectra  was  acquired  at 
each  new  specimen  position.  This  was  repeated  10  times  for 
each  specimen,  creating  a  small  Raman  image  —  1 84  X  11 
pm  in  size  and  containing  a  total  of  1260  Raman  spectra. 
The  total  acquisition  time  for  each  image  was  ~50  minutes. 
Specimens  were  kept  frozen  before  and  after  these  experi¬ 
ments. 

Two  different  approaches  were  used  for  the  analysis  of 
the  Raman  spectroscopic  data.  A  multivariate  approach 
known  as  factor  analysis'26,27’  was  used  to  determine 
whether  different  types  of  mineral  and/or  collagenous  ma¬ 
trix  were  present  within  the  areas  imaged.  This  methodol¬ 
ogy  is  used  to  extract  Raman  spectra  from  large  data  sets, 
such  as  the  hyperspectral  images  presented  here.  The  details 
of  application  to  bone  tissue  spectra  have  been  previously 
reported. (28,29)  Differences  in  mineral  or  organic  matrix 
composition  affect  the  Raman  spectrum  in  subtle  but  ob¬ 
servable  ways.  If  the  composition  of  the  mineral  and/or 
matrix  varies  across  the  area  imaged,  factor  analysis  is  able 
to  separate  the  spectra  of  the  different  mineral  and/or  matrix 
species.  Thus,  multiple  mineral  or  matrix  species,  each 
represented  by  its  own  Raman  spectrum,  may  be  observed 
within  a  single  image  area.  Alternatively,  if  there  is  no 
variation  in  mineral  and/or  matrix  composition  within  the 
area  imaged,  only  a  single  mineral  and/or  matrix  species 
will  be  observed.  Univariate  analysis,  in  the  form  of  band 
height  ratios,  was  also  performed  on  the  Raman  spectra,  a 
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FIG.  3.  Brtl  femoral  areal  BMD  is  significantly  less  than  WT  at  each 
age  (*p  <  0.05).  Additionally,  the  gain  in  Brtl  aBMD  between  2  and  6 
months  of  age  is  significantly  greater  than  that  of  WT  ( *p  <  0.05). 


sample  of  which  is  shown  in  Fig.  2.  Carbonate-to-phosphate 
ratios  were  computed  from  the  Uj  band  heights  of  the 
relative  carbonate  (~1070  cm-1)  and  phosphate  (~959 
cm-1)  peaks.  Full-width  half-maximum  (FWHM)  values  of 
the  u,  phosphate  band  were  calculated  to  provide  a  quali¬ 
tative  indication  of  the  quality,  or  crystallinity,  of  the  min¬ 
eral  phase/30-31 1  Representative  mineral-to-matrix  ratios 
were  calculated  by  ratioing  the  height  of  the  phosphate  peak 
to  the  height  of  the  amide  I  peak  (~1675  cm”1). 

Evaluation  of  aBMD  by  Piximus 

Male  mice  aged  1  (Brtl,  n  =  5;  WT,  n  =  5),  2  (Brtl,  n  = 
7;  WT,  n  =  8),  and  6  (Brtl.  n  —  9;  WT,  n  =  8)  months  were 
killed  by  lethal  injection.  Areal  BMD  (aBMD)  scans  of  Brtl 
mouse  femurs  were  compared  with  those  of  age-matched, 
littermate  controls  using  a  GE  Lunar  PIXImus2  (Madison, 
WI,  USA)  on  the  whole  mouse  specimen  and  the  internal 
calibration  standards  provided  by  GE. 

Statistical  analysis 


Table  1.  Geometric  Parameters  Derived  From  /aCT 


Age 

Genotype 

Cross-sectional  area 
(mm2) 

Cortical  thickness 
(mm) 

n 

1  Month 

Bill 

0.58  ±  0.11* 

0.16  ±  0.02* 

8 

WT 

0.73  ±  0.09 

0.19  ±  0.02 

13 

2  Month 

Brtl 

0.75  ±  0.07* 

0.20  ±  0.02* 

10 

WT 

1.05  ±  0.11 

0.24  ±  0.03 

12 

6  Month 

Brtl 

1.03  ±  0.17* 

0.28  ±  0.04 

11 

WT 

1.26  ±  0.23 

0.28  ±  0.04 

12 

12  Month 

Brtl 

1.18  ±  0.20 

0.26  ±  0.05 

5 

WT 

1.22  ±  0.13 

0.26  ±  0.04 

5 

Values  are  means  ±  SD. 
*p  <  0.05;  Brtl  vs.  WT. 


FIG.  4.  Bending  moments  in  the  anterior-posterior  direction  show 
that  the  structure  of  the  Brtl  femur  at  1,2,  and  6  months  is  less  resistant 
to  fracture  than  WT  (*/?  <  0.05). 


6  months  of  age,  Brtl  femurs  have  significantly  lower 
aBMD  than  WT.  Both  genotypes  show  significant  increases 
in  aBMD  at  each  age  point;  however.  Brtl  significantly 
improves  its  rate  of  mineralization  from  72.8%  and  76.2% 
of  WT  values  at  1  and  2  months  of  age,  respectively,  to 
88.9%  of  WT  aBMD  at  6  months  of  age  ( p  <  0.05). 


pG T  and  four-point  bending  data  were  statistically  analyzed 
using  multivariate  general  linear  models  for  each  age  group, 
with  significance  attributed  to  p  <  0.05.  A  nonparametric 
repeated  measures  model  (Wilcoxon  signed-rank  test)  was 
used  to  test  for  within-mice  differences  in  the  subset  of  /rCT 
data.  Raman  spectra  were  analyzed  using  factor  analysis  to 
identify  numbers  of  mineral  and  organic  species  in  the  matrix. 
Univariate  analysis  was  performed  for  phosphate  to  carbonate 
ratios  and  measures  of  full  width  half-maximum  of  the  ul 
phosphate  band.  aBMD  values  of  Brtl  and  WT  age-matched 
littermates  were  analyzed  using  a  Mann- Whitney  U  with  exact 
significances  test  for  genotypic  differences,  and  differences  in 
BMD  gain  between  2  and  6  months  were  assessed  using  a 
general  linear  model  and  testing  for  interactions  between  age 
and  genotype.  Data  are  presented  as  mean  ±  SD. 

RESULTS 

aBMD  of  Brtl  femurs 

Figure  3  shows  a  comparison  of  left  femoral  bone  den¬ 
sities  of  male  Brtl  femurs  and  WT  counterparts.  At  1,  2,  and 


Intraspecimen  variation  in  biomechanical  studies 

No  significant  differences  were  found  between  left  and 
right  femurs  in  the  subset  of  mice  analyzed  by  pCT.  There¬ 
fore,  where  necessary,  geometric  data  from  left  femurs  were 
used  to  predict  matrix  material  properties  from  right  femoral 
four-point  bending  data  in  the  same  mouse.  Two  1-month- 
old  Brtl  specimens  were  excluded  from  /u,CT  analysis  be¬ 
cause  of  artifact  during  scanning.  One  1 -month-old  Brtl 
specimen  was  fractured  during  the  preload  segment  of  load¬ 
ing  and  was  excluded  from  the  results. 

Structural  analysis 

Results  from  pGT  show  that  Brtl  mouse  femurs  have 
significantly  smaller  cross-sectional  areas  at  1,  2,  and  6 
months  than  do  their  WT  counterparts  (Table  1).  Cortical 
thickness  values  at  1  and  2  months  were  also  significantly 
reduced  (Table  1),  but  at  the  postpubertal  age  of  6  months, 
the  Brtl  mouse  femur  achieved  parity  with  controls.  Results 
from  bending  moments  of  inertia  in  the  anterior-posterior 
direction  (Fig.  4)  show  significant  differences  in  Brtl  com- 
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1  Month 


2  Months 


6  Months  12  Months 


FIG.  5.  Average  cross-sectional  views  of  WT 
femur  diaphysis  (heavy  outline)  superimposed 
on  top  of  Brtl  femur  diaphysis  (gray)  (2-mm- 
diameter  circle  shown  in  light  gray  for  refer¬ 
ence).  From  left  to  right,  1,  2,  6,  and  12  months. 


Table  2.  Mechanical  Properties  Derived  From  Four-Point  Bending  Tests 


Age 

Genotype 

Post-yield 

displacement 

(mm) 

Yield  energy  (N-mm) 

Failure  energy 
(N-mm) 

Yield  load  (N) 

1  Month 

Brtl 

0.1946  ±  0.077* 

0.27  ±  0.19 

1.90  ±  0.71* 

6.11 

±  2.37* 

(p  =  0.087) 

WT 

0.3829  ±  0.138 

0.51  ±  0.37 

5.96  ±  2.06 

9.26 

±  2.62 

2  Month 

Brtl 

0.1759  ±  0.071* 

0.32  ±  0.19 

3.03  ±  1.49* 

8.84 

±  1.85 

(P  = 

0.093) 

WT 

0.3016  ±  0.128 

0.34  ±  0.20 

7.37  ±  3.14 

11.24 

±  3.93 

6  Month 

Brtl 

0.1244  ±  0.056 

0.55  ±  0.23* 

3.96  ±  2.05 

17.12 

±  3.85 

WT 

0.1214  ±  0.88 

1.15  ±  0.92 

4.97  ±  3.33 

22.67 

±  10.52 

12  Month 

Brtl 

0.250  ±  0.014 

0.63  ±  0.46 

1.09  ±  0.71 

17.81 

±  8.93 

WT 

0.0572  ±  0.037 

0.61  ±  0.18 

2.16  ±  1.12 

19.07 

±  5.46 

Values  are  means  ±  SD. 
*p  <  0.05;  Brtl  vs.  WT. 


pared  with  controls  at  1,  2,  and  6  months.  These  results 
suggest  that,  based  on  structural  geometry  alone,  indepen¬ 
dent  of  matrix-level  material  properties,  the  Brtl  femur 
should  be  less  resistant  to  an  applied  load  than  WT  through 
6  months  of  age.  Average  cross-sectional  views  shown  in 
Fig.  5  summarize  these  results. 

Mechanical  testing 

Femurs  from  young  Brtl  mice  show  a  more  brittle  phe¬ 
notype  compared  with  WT,  as  shown  by  reduced  post-yield 
displacement  and  reduced  failure  energy  at  1  and  2  months 
(Table  2).  Brtl  mice  show  reduced  yield  (Table  2)  and 
maximum  loads  (Fig.  6)  at  1  and  2  months  of  age  but  equal 
WT  values  at  6  and  12  months.  Brtl  femurs  are  significantly 
less  stiff  at  2  months  but  comparable  with  WT  at  6  and  12 
months  (Fig.  7).  These  findings  show  a  postpubertal  alter¬ 
ation  whereby  mechanical  integrity  of  the  femur  is  substan¬ 
tially  improved  to  the  level  of  the  WT.  These  findings  are 
contrary  to  geometric  results  from  pCT,  which  suggested 
that  based  on  structural  geometry  alone,  Brtl  mice  at  6 
months  should  be  less  stiff  and  less  strong  than  WT.  There¬ 
fore,  we  hypothesize  that  alterations  occur  in  the  Brtl  matrix 
after  puberty  that  dramatically  improve  the  matrix  material 
properties  measured  at  6  months  of  age. 

Predicted  matrix  properties 

When  four-point  bending  results  are  normalized  by  bend¬ 
ing  moments  of  inertia  to  predict  matrix  material  properties, 
predicted  elastic  modulus  of  the  6-month-old  Brtl  femur 
was  significantly  increased  compared  with  WT  (Table  3). 


FIG.  6.  Maximum  load  is  reduced  in  Brtl  femurs  compared  with  WT 
at  1  and  2  months,  while  at  6  and  12  months,  no  significant  differences 
were  found.  (*  p  <  0.05) 


Additionally,  6-month-old  Brtl  femurs  showed  a  strong 
trend  toward  increased  predicted  ultimate  strength  over  WT 

{p  =  0.068). 

Raman  spectroscopy 

Predicted  matrix  properties  suggest  a  dynamic  period  of 
alteration  leading  to  increased  material  properties  at  the 
matrix  level.  To  evaluate  the  potential  mechanisms  causing 
this  change,  specimens  from  2-  and  6-month-old  animals 
were  evaluated  by  Raman  spectroscopy  (Table  4).  One 
2-month-old  Brtl  specimen  was  dropped  from  the  analysis 
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FIG.  7.  Stiffness  of  the  Brtl  femur  achieves  equivalence  to  WT  at  6 
months  of  age.(*  p  <  0.05) 


Table  3.  Predicted  Ultimate  Strength  and 
Predicted  Elastic  Modulus 


Age 

Genotype 

Ultimate  strength 
(MPa) 

Elastic  modulus 
(MPa) 

1  Month 

Brtl 

43.03  ±  10.28 

612.65  ±  192.42 

WT 

48.58  ±  7.21 

624.98  ±  189.92 

2  Month 

Brtl 

58.43  ±  12.12 

862.44  ±  267.28 

WT 

57.85  ±  12.66 

645.24  ±  159.65 

6  Month 

Brtl 

74.20  ±  22.27 

( p  =  0.068) 

1180.99  ±  322.81* 

WT 

58.33  ±  17.20 

686.01  ±  133.9 

12  Month 

Brtl 

29.83  ±  12.89 

611.35  ±  147.15 

WT 

44.04  ±16.15 

660.35  ±  268.39 

Values  are  means  ±  SD. 
*p  <  0.05;  Brtl  vs.  WT. 


because  of  a  poor  signal-to-noise  ratio,  making  baselining  of 
the  spectra  difficult.  One  6-month-old  Brtl  specimen  was 
dropped  from  the  analysis  because  of  a  failed  Q-test.  Anal¬ 
ysis  of  both  the  mineral  subregion  and  organic  subregion  of 
the  Raman  spectra  revealed  no  trends  in  the  number  of 
mineral  or  matrix  factors  with  respect  to  disease  status. 
Six-month-old  Brtl  specimens  seemed  to  have  higher 
carbonate-to-phosphate  ratios  than  WT  specimens,  and  this 
trend  was  not  evident  at  2  months.  Measures  of  phosphate 
band  width  show  no  difference  in  mineral  crystallinity  with 
respect  to  age  or  disease  state.  However,  results  suggest  that 
the  Brtl  mouse  has  a  significantly  reduced  mineral-to-matrix 
ratio  compared  with  WT  at  6  months  of  age. 

DISCUSSION 

The  Brtl  mouse  was  generated  as  a  model  for  01  type  IV 
by  introducing  a  classic  glycine  substitution  in  one  collal 
allele. <12)  Individuals  with  type  IV  OI  have  moderately 
severe  skeletal  fragility,  with  decreased  fracture  susceptibil¬ 
ity  after  puberty.  The  Brtl  mouse  results  presented  here 
model  the  clinical  course  of  the  disorder  in  children.  The 
femurs  of  prepubertal  1 -month-old  and  pubertal  2-month- 
old  Brtl  mice  were  structurally  less  able  to  withstand  a  load 
than  age-matched  WT  controls.  This  finding  was  verified 


through  whole  bone  mechanical  tests.  At  the  postpubertal 
age  of  6  months,  however,  the  Brtl  mouse  femurs  had  whole 
bone  strength  and  stiffness  comparable  with  WT,  despite 
having  a  bone  geometrical  organization  less  resistant  to 
applied  loads.  Predicted  material  properties  of  elastic  mod¬ 
ulus  and  ultimate  strength  are  increased,  accounting  for 
increases  in  whole  bone  strength  through  an  increase  in 
matrix-level  material  properties  without  corresponding 
structural  adaptations.  These  findings  are  bolstered  by  re¬ 
sults  from  the  aBMD  measurements.  Brtl  aBMD  was  sig¬ 
nificantly  lower  than  WT  at  all  ages;  however,  Brtl  aBMD 
increased  between  2  and  6  months  at  a  rate  greater  than 
associated  with  the  pubertal  changes  in  WT  mice.  Raman 
results  suggest  a  change  in  the  balance  of  ultrastructural 
components,  as  mineral  matrix  ratios  adjust  in  Brtl  from 
being  nonsignificantly  larger  than  WT  at  2  months  to  being 
significantly  lower  at  6  months.  The  2-month-old  Raman 
data  do  not  contradict  the  findings  of  reduced  Brtl  aBMD  at 
all  ages,  because  Raman  measures  the  ratio,  not  total 
amount,  of  mineral  and  collagen  constituencies  at  the  ultra- 
structural  level.  If  Brtl  produces  less  collagen  than  WT  at 
the  ultrastructural  level,  it  could  still  show  a  higher  mineral: 
matrix  ratio  by  Raman,  but  a  lower  aBMD  because  of  the 
difference  in  the  two  measurement  techniques. 

Together,  these  findings  suggest  a  period  of  adaptation 
between  2  and  6  months  of  age  in  which  bone  stiffness  and 
strength  attain  WT  values  through  a  strengthening  mecha¬ 
nism  at  the  level  of  the  extracellular  matrix  rather  than 
through  alterations  in  whole  bone  structural  layout. 

Significant  further  gains  in  bone  strength  are  not  demon¬ 
strated  beyond  6  months  of  age.  This  finding  is  consistent 
with  other  studies  showing  minimal  gains  in  bone  strength 
and  stiffness  beyond  18  weeks  of  age  in  the  C57BL/6J 
mouse* 321  and  constant  peak  torque  between  6  and  12 
months  in  C57BL/6J  X  C3H/He  mice.(33)  Furthermore,  WT 
ductility  is  decreased  at  6  and  12  months  compared  with 
earlier  ages,  rendering  differences  in  failure  energy  and 
post-yield  displacement  no  longer  significant.  This  suggests 
a  natural  aging  process  in  the  WT  mouse  leading  to  age- 
related  decreases  in  ductility.  Similar  decrements  in  femoral 
ductility  have  been  shown  in  SAMR1  control  mice  and  the 
SAMP6  mouse  model  for  senile  osteoporosis  between  4  and 
12  months  of  age.*341 

Other  models  of  OI  have  demonstrated  adaptations  in 
response  to  a  matrix-altering  mutation.  The  Movl3  mouse 
carries  a  pro  virus  preventing  transcription  of  the  a  1(1)  col¬ 
lagen  gene.'35’  Mice  heterozygous  for  the  Movl3  mutation 
show  a  50%  decrease  in  normal  type  I  collagen  and  serve  as 
a  model  for  OI  type  I,*81  in  which  patients  produce  a  reduced 
quantity  of  structurally  normal  type  I  collagen.  At  8  weeks 
of  age,  Movl3  femurs  are  more  brittle  than  WT  because  of 
microstructural  differences  between  the  two  animals. (36) 
Furthermore,  microbeams  machined  from  the  Movl3  femur 
show  decreased  matrix  material  properties  and  a  reduced 
fatigue  life  compared  with  controls.*211  To  compensate  for 
these  mechanical  impairments,  between  the  ages  of  8  and  15 
weeks,  Movl3  mice  deposit  new  bone  matrix  on  the  peri¬ 
osteal  surfaces  of  the  femur,  significantly  increasing  their 
cortical  area  and  bending  moment  of  inertia.*371  This  struc¬ 
tural  adaptation  significantly  improves  the  whole  bone 
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Table  4.  Results  From  Raman  Spectroscopy  Suggest  an  Alteration  in  Mineral-to-Matrix  Ratios  at  6  Months  of  Age 


Age 

Genotype 

Mineral  factors 

C03:P04  ratio 

FWHM  vl 

P04  band 

Mineral/matrix  ratio 

2  Month 

Brtl 

1. 1, 1 

0.144  ±  0.011 

16.25  ±  1.17 

12.50  ±  4.79 

WT 

1. 1, 1 

0.162  ±  0.068 

16.94  ±  0.46 

8.45  ±  1.65 

6  Month 

Brtl 

1, 2, 1 

0.209  ±  0.027 

17.15  ±  2.16 

5.75  ±  0.11* 

WT 

2, 1, 1 

0.193  ±  0.056 

17.13  ±  0.65 

11.17  ±  2.99 

*  p  <  0.05;  Brt  vs.  WT. 


bending  strength  of  the  Movl3  mouse  to  a  level  greater  than 
control  specimens. 

The  oim/oim  mouse  has  a  naturally  occurring  mutation 
that  prevents  pro-a2(I)  chains  from  incorporating  into  type 
I  collagen.  The  oim/oim  extracellular  matrix  does  not  con¬ 
tain  any  normal  type  I  collagen  but  is  instead  composed  of 
al(I)  homotrimers.'9’  This  model  for  01  type  III  show 
reduced  indices  of  whole  bone  strength  compared  with  WT 
at  all  ages/33'38’  Whereas  geometric-based  indices  of  bone 
strength  suggest  an  improved  bone  geometry  in  older  oim/ 
oim  mice,  this  alteration  fails  to  improve  bone  strength  to 
the  level  of  the  WT  femur.<33) 

Our  current  findings  show  improvements  in  whole  bone 
strength  through  improved  matrix  level  material  properties 
rather  than  geometrical  adaptation.  Raman  spectroscopy 
was  used  as  a  means  to  investigate  specific  mineral-matrix 
interactions  to  help  determine  what  alterations  lead  to  the 
improvement  in  whole  bone  strength.  Our  results  yielded  no 
difference  in  number  of  mineral  factors  or  mineral  crystal¬ 
linity  with  respect  to  disease  state.  Rather,  a  higher  C03: 
P04  ratio  and  reduced  mineral-to-matrix  ratios  were  ob¬ 
served  in  6-month-old  Brtl  specimens  compared  with  WT. 
At  this  hierarchical  level,  it  is  unclear  how  these  fundamen¬ 
tal  changes  might  affect  the  mechanical  properties  of  the 
tissue.  Bohic  et  al.'39’  have  shown  similar  increases  in 
carbonation  of  tooth  enamel  from  a  patient  with  dentino¬ 
genesis  imperfecta  compared  with  an  age-matched  control. 
These  findings  differ  from  investigations  in  the  oim/oim 
mouse.  Camacho  et  al/40’  used  Fourier-transform  infrared 
spectroscopy,  an  imaging  modality  that  provides  similar 
information  to  Raman  spectroscopy,  to  show  differences  in 
phosphate  species'40’  and  crystallinity138’  compared  with 
controls,  oim/oim  mice  show  trends  opposite  Brtl  in  C03: 
P04  and  minerahmatrix  ratios  compared  with  the  milder 
phenotype  of  oim/+  at  6  months  and  compared  with  WT  at 
1  year/40’  Differences  in  these  findings  may  be  because  of 
differences  in  the  collagen  defect.  Matrix  containing  only 
a  1(1)  homotrimeric  molecules  in  the  oim/oim  mice  may 
induce  a  different  mineralization  template  than  formed  by 
Brtl  matrix,  which  contains  a  mixture  of  normal  collagen 
and  mutant  collagen  with  a  glycine  to  cysteine  substitution. 
Neither  Movl3  nor  oim/oim  mice  show  a  compensatory 
mechanism  at  the  level  of  the  extracellular  matrix  capable  of 
significantly  improving  whole  bone  mechanical  strength  to 
the  level  of  controls,  independent  of  changes  in  whole  bone 
geometry.  The  femur  of  the  Brtl  mouse  demonstrates  this 
type  of  alteration. 

Vetter  et  al/41’  have  demonstrated  small  but  significant 
decreases  in  apatite  crystal  size  in  children  and  adolescents 


with  type  IV  01  compared  with  controls.  In  the  current 
study,  no  changes  in  mineral  crystallinity  were  found  be¬ 
tween  Brtl  and  WT  mice;  however,  we  cannot  definitively 
rule  out  this  possibility  because  no  studies  have  yet  cali¬ 
brated  Raman  band  width  to  mineral  crystallinity  as  mea¬ 
sured  by  X-ray  diffraction. 

Another  limitation  to  this  study  is  the  small  sample  size 
used  in  the  Raman  analysis.  Additionally,  milling  of  the 
microspecimens  may  have  introduced  artifact  because  of 
mechanical  deformation  of  the  bone/22'25’  Further  analysis 
of  the  Brtl  extracellular  matrix  is  warranted.  Changes  in 
tissue  microarchitecture,  such  as  percent  porosity  and  rela¬ 
tive  balance  of  woven  versus  lamellar  bone,  have  been 
demonstrated  to  be  important  mediators  of  mechanical  be¬ 
havior  in  the  Movl3  mouse/36’  Alterations  in  porosity  or 
remodeling  spaces,  as  seen  in  human  cases  of  OI,  may  lead 
to  the  inclusion  or  deduction  of  stress  risers  within  the  bone 
matrix,  which  would  be  unaccounted  for  in  our  current 
calculations  for  predicted  strength  and  modulus.  Further¬ 
more,  woven  and  lamellar  bone  distribution  may  change  as 
Brtl  ages,  altering  the  mechanical  profile  of  the  bone.  In 
addition  to  microarchitecture,  further  investigation  of  the 
mineral  and  organic  phase  of  the  Brtl  bone  through  methods 
such  as  small-angle  X-ray  scattering'42’  and  transmission 
electron  microscopy'43’  may  lend  insight  into  mineral  crys¬ 
tal  size,  shape,  orientation,  and  arrangement  within  the 
matrix.  These  methods  may  be  needed  to  lend  further  un¬ 
derstanding  into  the  results  generated  through  Raman  spec¬ 
troscopy. 

Results  from  this  study  and  others'37’  suggest  two  general 
avenues  of  therapeutics  for  the  disease:  (1)  improving  geo¬ 
metric  resistance  to  fracture,  as  in  the  Movl3  mouse,  and 
(2)  improving  the  material  properties  of  the  extracellular 
matrix  itself,  as  shown  here  in  the  Brtl  mouse.  Hormones 
such  as  parathyroid  hormone  (PTH)  or  other  anabolic  drugs 
may  be  useful  as  a  means  of  improving  geometric  resistance 
to  fracture  by  adding  bone  to  the  cortex  as  seen  in  both 
animal'4445’  and  human'46’  histomorphometric  studies. 
Bisphosphonates  have  been  hypothesized  to  change  the 
material  properties  of  bone,  possibly  through  increases  in 
secondary  matrix  mineralization.'47-49’  Indeed,  alendronate, 
a  third-generation  bisphosphonate,  has  been  demonstrated 
to  reduce  the  number  of  fractures  and  increase  femoral 
metaphyseal  density  in  the  growing  oim/oim  mouse  model 
for  OI/50’  The  Brtl  mouse  is  an  appropriate  model  in  which 
to  address  these  two  avenues  of  therapeutics  and  their 
potential  implications  on  OI  bone  mechanics  and  geometry. 

The  findings  presented  here  parallel  clinically  observed 
decreases  in  OI  patient  fracture  rate  after  puberty  without 
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administration  of  drug  therapies. (7)  The  mechanism  of  this 
improvement  in  patient  bone  integrity  is  not  well  under¬ 
stood.  The  postpubertal  decrease  in  fracture  rate  occurs  with 
both  glycine  substitutions  and  exon  skipping  defects  of 
collagen.  Our  results  are  the  first  to  suggest  a  possible 
matrix-level  adaptation  leading  to  improved  whole  bone 
strength  without  corresponding  increases  in  structural  ge¬ 
ometry  in  a  model  of  01.  It  is  apparent  that  nature  has 
developed  mechanisms  to  adapt  whole  bone  integrity  to  the 
fragility  induced  by  type  I  collagen  mutations.  Further  in¬ 
vestigations  are  underway  to  determine  the  molecular  mech¬ 
anisms  triggering  this  adaptation  in  Brtl  matrix  properties 
and  the  nature  of  the  changes  in  postpubertal  matrix  com¬ 
position.  Exploiting  these  naturally  occurring  adaptation  strat¬ 
egies  may  provide  novel  therapeutic  interventions  for  01. 

ACKNOWLEDGMENTS 

This  study  was  supported  by  National  Institutes  of  Health 
AR46024  (SAG),  National  Institutes  of  Health  AR47969 
(MDM)  and  NICHD,  National  Institutes  of  Health  Intramu¬ 
ral  Funding  (JCM).  KMK  is  supported  by  a  Whitaker  Foun¬ 
dation  Graduate  Student  Fellowship.  We  thank  EP  Franken- 
burg,  M  Spurchise,  R  Rae  Jr,  J  Hall,  S  Volkmann,  D  Kayner 
and  C  Roehm  for  their  valuable  contributions  to  this  study. 

REFERENCES 

1.  Byers  PH  2002  Osteogenesis  imperfecta.  In:  Royce  PM,  Stein- 
mann  B  (eds.)  Connective  Tissue  and  Its  Heritable  Disorders,  2nd 
ed.  Wiley-Liss,  Inc.,  New  York,  NY,  USA,  pp.  385-430. 

2.  Kuivaniemi  H,  Tromp  G,  Prockop  DJ  1991  Mutations  in  collagen 
genes:  Causes  of  rare  and  some  common  diseases  in  humans. 
FASEB  J  5:2052-2060. 

3.  Glorieux  FH,  Rauch  F,  Plotkin  H,  Ward  L,  Travers  R,  Roughley  P, 
Lalic  L,  Glorieux  DF,  Fassier  F,  Bishop  NJ  2000  Type  V  osteo¬ 
genesis  imperfecta:  A  new  form  of  brittle  bone  disease.  J  Bone 
Miner  Res  15:1650-1658. 

4.  Labuda  M,  Morissette  J,  Ward  LM,  Rauch  F,  Lalic  L,  Roughley  PJ, 
Glorieux  FH  2002  Osteogenesis  imperfecta  type  VII  maps  to  the 
short  arm  of  chromosome  3.  Bone  31:19-25. 

5.  Ward  LM,  Rauch  F,  Travers  R,  Chabot  G,  Azouz  EM,  Lalic  L, 
Roughley  PJ,  Glorieux  PH  2002  Osteogenesis  imperfecta  type  VII: 
An  autosomal  recessive  form  of  brittle  bone  disease.  Bone  31:12- 
18. 

6.  Glorieux  FH,  Ward  LM,  Rauch  F,  Lalic  L,  Roughley  PJ,  Travers 
R  2002  Osteogenesis  imperfecta  type  VI:  A  form  of  brittle  bone 
disease  with  a  mineralization  defect.  J  Bone  Miner  Res  17:30-38. 

7.  Marini  JC  1988  Osteogenesis  imperfecta:  Comprehensive  manage¬ 
ment.  Adv  Pediatr  35:391-426. 

8.  Bonadio  J,  Saunders  TL,  Tsai  E,  Goldstein  SA,  Morris-Wiman  J, 
Brinkley  L,  Dolan  DF,  Altschuler  RA,  Hawkins  JE  Jr,  Bateman  JF 
1990  Transgenic  mouse  model  of  the  mild  dominant  form  of 
osteogenesis  imperfecta.  Proc  Natl  Acad  Sci  USA  87:7145-7149. 

9.  Chipman  SD,  Sweet  HO,  McBride  DJ  Jr,  Davisson  MT,  Marks  SC 
Jr,  Shuldiner  AR,  Wenstrup  RJ,  Rowe  DW,  Shapiro  JR  1993 
Defective  pro  alpha  2(1)  collagen  synthesis  in  a  recessive  mutation 
in  mice:  A  model  of  human  osteogenesis  imperfecta.  Proc  Natl 
Acad  Sci  USA  90:1701-1705. 

10.  Saban  J,  Zussman  MA,  Havey  R,  Patwardhan  AG,  Schneider  GB, 
King  D  1996  Heterozygous  oim  mice  exhibit  a  mild  form  of 
osteogenesis  imperfecta.  Bone  19:575-579. 

11.  Khillan  JS,  Olsen  AS,  Kontusaari  S,  Sokolov  B,  Prockop  DJ  1991 
Transgenic  mice  that  express  a  mini-gene  version  of  the  human 
gene  for  type  I  procollagen  (COL1A1)  develop  a  phenotype  re¬ 
sembling  a  lethal  form  of  osteogenesis  imperfecta.  J  Biol  Chem 
266:23373-23379. 

12.  Forlino  A,  Porter  FD,  Lee  EJ,  Westphal  H,  Marini  JC  1999  Use  of 
the  Cre/lox  recombination  system  to  develop  a  non-lethal  knock-in 
murine  model  for  osteogenesis  imperfecta  with  an  a  1(1)  G349C 


substitution.  Variability  in  phenotype  in  BrtllV  mice.  J  Biol  Chem 
274:37923-37931. 

13.  Sarafova  AP,  Choi  H,  Forlino  A,  Gajko  A,  Cabral  WA,  Tosi  L, 
Reing  CM,  Marini  JC  1998  Three  novel  type  I  collagen  mutations 
in  osteogenesis  imperfecta  type  IV  probands  are  associated  with 
discrepancies  between  electrophoretic  migration  of  osteoblast  and 
fibroblast  collagen.  Hum  Mutat  11:395-403. 

14.  Feldkamp  LA,  Goldstein  SA,  Parfitt  AM,  Jesion  G,  Kleerekoper  M 
1989  The  direct  examination  of  three-dimensional  bone  architec¬ 
ture  in  vitro  by  computed  tomography.  J  Bone  Miner  Res  4:3-11. 

15.  Kuhn  JL,  Goldstein  SA,  Feldkamp  LA,  Goulet  RW,  Jesion  G  1990 
Evaluation  of  a  microcomputed  tomography  system  to  study  tra¬ 
becular  bone  structure.  J  Orthop  Res  8:833-842. 

16.  Goulet  RW  1993  The  quantification  of  the  structure  and  mechan¬ 
ical  properties  of  trabecular  bone.  Ph.D.  thesis,  University  of 
Michigan,  Ann  Arbor,  MI. 

17.  Turner  CH,  Burr  DB  1993  Basic  biomechanical  measurements  of 
bone:  A  tutorial.  Bone  14:595-608. 

18.  Carden  A,  Morris  MD  2000  Application  of  vibrational  spectros¬ 
copy  to  the  study  of  mineralized  tissues.  J  Biomed  Opt  5:259-268. 

19.  Choi  K,  Goldstein  SA  1992  A  comparison  of  the  fatigue  behavior 
of  human  trabecular  and  cortical  bone  tissue.  J  Biomech  25:1371- 
1381. 

20.  Kuhn  JL,  Goldstein  SA,  Choi  K,  London  M,  Feldkamp  LA,  Mat¬ 
thews  LS  1989  Comparison  of  the  trabecular  and  cortical  tissue 
moduli  from  human  iliac  crests.  J  Orthop  Res  7:876-884. 

21.  Jepsen  KJ,  Schaffler  MB,  Kuhn  JL,  Goulet  RW,  Bonadio  J,  Gold¬ 
stein  SA  1997  Type  I  collagen  mutation  alters  the  strength  and 
fatigue  behavior  of  Movl3  cortical  tissue.  J  Biomech  30: 1 141 — 
1 147. 

22.  Carden  A,  Rajachar  RM,  Morris  MD,  Kohn  DH  2003  Ultrastruc- 
tural  changes  accompanying  the  mechanical  deformation  of  bone 
tissue:  A  Raman  imaging  study.  Calcif  Tissue  Int  72:166-175. 

23.  Pezzuti  JA,  Morris  MD,  Bonadio  JF,  Goldstein  SA  1998  Hyper- 
spectral  Raman  imaging  of  bone  growth  and  regrowth  chemistry. 
Proc  SPIE  3261:270-276. 

24.  Timlin  JA,  Carden  A,  Morris  MD,  Bonadio  JF,  Hoffler  CE,  Ko- 
zloff  KM,  Goldstein  SA  1999  Spatial  distribution  of  phosphate 
species  in  mature  and  newly  generated  mammalian  bone  by  hy- 
perspectral  Raman  imaging.  J  Biomed  Opt  4:28-34. 

25.  Timlin  JA,  Carden  A,  Moms  MD,  Rajachar  RM,  Kohn  DH  2000 
Raman  spectroscopic  imaging  markers  for  fatigue-related  micro¬ 
damage  in  bovine  bone.  Anal  Chem  72:2229-2236. 

26.  Malinowski  ER  1991  Factor  Analysis  in  Chemistry,  2nd  ed.  John 
Wiley  &  Sons,  Inc.,  New  York,  NY,  USA. 

27.  Reyment  R,  Joreskog  KG  1996  Applied  Factor  Analysis  in  the 
Natural  Sciences.  Cambridge  University  Press,  New  York,  NY, 
USA. 

28.  Morris  MD,  Crane  NJ,  Gomez  LE,  Ignelzi  MA  Jr  2004  Compat¬ 
ibility  of  staining  protocols  for  bone  tissue  with  Raman  imaging. 
Calcif  Tissue  Int  (in  press). 

29.  Tamowski  CP,  Ignelzi  MA  Jr,  Morris  MD  2002  Mineralization  of 
developing  mouse  calvaria  as  revealed  by  Raman  microspectros¬ 
copy.  J  Bone  Miner  Res  17:1118-1126. 

30.  Freeman  JJ,  Wopenka  B,  Silva  MJ,  Pasteris  JD  2001  Raman 
spectroscopic  detection  of  changes  in  bioapatite  in  mouse  femora 
as  a  function  of  age  and  in  vitro  fluoride  treatment.  Calcif  Tissue 
Int  68:156-162. 

31.  de  Mul  FF,  Hottenhuis  MH,  Bouter  P,  Greve  J,  Arends  J,  ten  Bosch 
JJ  1986  Micro-Raman  line  broadening  in  synthetic  carbonated 
hydroxyapatite.  J  Dent  Res  65:437-440. 

32.  Ferguson  VL,  Ayers  RA,  Bateman  TA,  Simske  SJ  2000  Develop¬ 
ment  of  endogenous  osteoporosis  in  male  C57BL/6J  mice.  J  Bone 
Miner  Res  15:S461. 

33.  McBride  DJ  Jr,  Shapiro  JR,  Dunn  MG  1998  Bone  geometry  and 
strength  measurements  in  aging  mice  with  the  oim  mutation.  Calcif 
Tissue  Int  62:172-176. 

34.  Silva  MJ,  Brodt  MD,  Ettner  SL  2002  Long  bones  from  the  senes¬ 
cence  accelerated  mouse  SAMP6  have  increased  size  but  reduced 
whole-bone  strength  and  resistance  to  fracture.  J  Bone  Miner  Res 
17:1597-1603. 

35.  Schnieke  A,  Dziadek  M,  Bateman  J,  Mascara  T,  Harbers  K, 
Gelinas  R,  Jaenisch  R  1987  Introduction  of  the  human  pro  alpha 
1(1)  collagen  gene  into  pro  alpha  l(I)-deficient  Mov-13  mouse  cells 
leads  to  formation  of  functional  mouse-human  hybrid  type  I  col¬ 
lagen.  Proc  Natl  Acad  Sci  USA  84:764-768. 


KOZLOFF  ET  AL. 


622 


36.  Jepsen  KJ,  Goldstein  SA,  Kuhn  JL,  Schaffler  MB,  Bonadio  J  1996 
Type-I  collagen  mutation  compromises  the  post-yield  behavior  of 
Movl3  long  bone.  J  Orthop  Res  14:493-499. 

37.  Bonadio  J,  Jepsen  KJ,  Mansoura  MK,  Jaenisch  R,  Kuhn  JL, 
Goldstein  SA  1993  A  murine  skeletal  adaptation  that  significantly 
increases  cortical  bone  mechanical  properties.  Implications  for 
human  skeletal  fragility.  J  Clin  Invest  92:1697-1705. 

38.  Camacho  NP,  Hou  L,  Toledano  TR,  Ilg  WA,  Brayton  CF,  Raggio 
CL,  Root  L,  Boskey  AL  1999  The  material  basis  for  reduced 
mechanical  properties  in  oim  mice  bones.  J  Bone  Miner  Res 
14:264-272. 

39.  Bohic  S,  Heymann  D,  Pouezat  JA,  Gauthier  O,  Daculsi  G  1998 
Transmission  FT-IR  microspectroscopy  of  mineral  phases  in  cal¬ 
cified  tissues.  Comptes  Rendus  de  l’Academie  des  Sciences  Serie 
III,  Sciences  de  la  Vie  321:865-876. 

40.  Camacho  NP,  Landis  WJ,  Boskey  AL  1996  Mineral  changes  in  a 
mouse  model  of  osteogenesis  imperfecta  detected  by  Fourier  trans¬ 
form  infrared  microscopy.  Connect  Tissue  Res  35:259-265. 

41.  Vetter  U,  Eanes  ED,  Kopp  JB,  Termine  JD,  Robey  PG  1991 
Changes  in  apatite  crystal  size  in  bones  of  patients  with  osteogen¬ 
esis  imperfecta.  Calcif  Tissue  Int  49:248-250. 

42.  Fratzl  P,  Paris  O,  Klaushofer  K,  Landis  WJ  1996  Bone  mineral¬ 
ization  in  an  osteogenesis  imperfecta  mouse  model  studied  by 
small-angle  x-ray  scattering.  J  Clin  Invest  97:396-402. 

43.  McKee  MD,  Nanci  A  1995  Postembedding  colloidal-gold  immu- 
nocytochemistry  of  noncollagenous  extracellular  matrix  proteins  in 
mineralized  tissues.  Microsc  Res  Tech  31:44-62. 

44.  Hirano  T,  Burr  DB,  Turner  CH,  Sato  M,  Cain  RL,  Hock  JM  1999 
Anabolic  effects  of  human  biosynthetic  parathyroid  hormone  frag¬ 
ment  (1-34),  LY333334,  on  remodeling  and  mechanical  properties 
of  cortical  bone  in  rabbits.  J  Bone  Miner  Res  14:536-545. 

45.  Wronski  TJ,  Yen  CF  1994  Anabolic  effects  of  parathyroid  hor¬ 
mone  on  cortical  bone  in  ovariectomized  rats.  Bone  15:51-58. 

46.  Dempster  DW,  Cosman  F,  Kurland  ES,  Zhou  H,  Nieves  J, 
Woelfert  L,  Shane  E,  Plavetic  K,  Muller  R,  Bilezikian  J,  Lindsay 
R  2001  Effects  of  daily  treatment  with  parathyroid  hormone  on 


bone  microarchitecture  and  turnover  in  patients  with  osteoporosis: 
A  paired  biopsy  study.  J  Bone  Miner  Res  16:1846-1853. 

47.  Chavassieux  PM,  Arlot  ME,  Reda  C,  Wei  L,  Yates  AJ,  Meunier  PJ 
1997  Histomorphometric  assessment  of  the  long-term  effects  of 
alendronate  on  bone  quality  and  remodeling  in  patients  with  os¬ 
teoporosis.  J  Clin  Invest  100:1475-1480. 

48.  Meunier  PJ,  Boivin  G  1997  Bone  mineral  density  reflects  bone 
mass  but  also  the  degree  of  mineralization  of  bone:  Therapeutic 
implications.  Bone  21:373-377. 

49.  Bourrin  S,  Ammann  P,  Bonjour  JP,  Rizzoli  R  2002  Recovery  of 
proximal  tibia  bone  mineral  density  and  strength,  but  not  cancel¬ 
lous  bone  architecture,  after  long-term  bisphosphonate  or  selective 
estrogen  receptor  modulator  therapy  in  aged  rats.  Bone  30:195- 
200. 

50.  Camacho  NP,  Raggio  CL,  Doty  SB,  Root  L,  Zraick  V,  Ilg  WA, 
Toledano  TR,  Boskey  AL  2001  A  controlled  study  of  the  effects  of 
alendronate  in  a  growing  mouse  model  of  osteogenesis  imperfecta. 
Calcif  Tissue  Int  69:94-101. 


Address  reprint  requests  to: 
Steven  A  Goldstein,  PhD 
Orthopaedic  Research  Laboratories 
Department  of  Orthopaedic  Surgery 
University  of  Michigan 
Room  G-161 
400  North  Ingalls 
Ann  Arbor,  MI  48109-0486,  USA 
E-mail:  Steve gld@  umich.edu 

Received  in  original  form  May  21,  2003;  in  revised  form  October 
28,  2003;  accepted  November  14,  2003. 


(guiui)  VSO  (ZUIUI)  VSO 


Femur  Mid-diaphyseal  Cross-sectional  Area  by  Diet  and  Exercise 


Cross  Sectional  Area  by  Diet  and  Exercise  (Female  C57  Only) 


1 

Low  |  Nor 

Y 

4* 

3S 

mal  |  High 

0 

Cross  Sectional  Area  by  Diet  and  Exercise  (Male  C57  Only) 


1.2 


Cortical  Thickness  (mm)  Cortical  Thickness  (mm) 


Femur  Mid-diaphyseal  Cortical  Thickness  by  Diet  and  Exercise 


Cortical  Thickness  by  Diet  and  Exercise  (Female  C57  Only) 


Cortical  Thickness  by  Diet  and  Exercise  (Male  C57  Only) 


Cortical  Thickness  by  Diet  and  Exercise  (Female  C3  Only) 


Cortical  Thickness  by  Diet  and  Exercise  (Male  C3  Only) 


0.4 


0.4 


0.35 

0.3 

0.25 


0.2 


0.15 

0.1 

0.05 

0 


Normal 


Yes 


I 


Normal 

No 


High 


lyy  (mm4)  lyy  (mm4) 


Femur  Mid-diaphyseal  Bending  Moment  of  Inertia  by  Diet  and  Exercise 


Bending  Moment  of  Inertia  by  Diet  and  Exercise  (Female  C57  Only) 


Bending  Moment  of  Inertia  by  Diet  and  Exercise  (Male  C57  Only) 


Bending  Moment  of  Inertia  by  Diet  and  Exercise  (Female  C3  Only) 


Bending  Moment  of  Inertia  by  Diet  and  Exercise  (Male  C3  Only) 


Shaping  Factor  Shaping  Factor 


Femur  Mid-diaphyseal  Shaping  Factor  by  Diet  and  Exercise 


Shaping  Factor  by  Diet  and  Exercise  (Female  C57  Only) 


Shaping  Factor  by  Diet  and  Exercise  (Male  C57  Only) 


Shaping  Factor  by  Diet  and  Exercise  (Female  C3  Only) 


Shaping  Factor  by  Diet  and  Exercise  (Male  C3  Only) 


0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 


i 


Normal 

Yes 


High 


Low 


BMD  (mg/cc)  BMD  (mg/cc) 


Femur  Mid-diaphyseal  Bone  Mineral  Density  by  Diet  and  Exercise 

BMD  by  Diet  and  Exercise  (Female  C57  only)  BMD  by  Diet  and  Exercise  (Male  C57  only) 


600 


500 


400 


BMD  by  Diet  and  Exercise  (Female  C3  Only) 


600 


500 


400 


BMD  by  Diet  and  Exercise  (Male  C3  Only) 


TMD  (mg/cc)  TMD  (mg/cc) 


Femur  Mid-diaphyseal  Tissue  Mineral  Density  by  Diet  and  Exercise 


Tissue  Mineral  Density  by  Diet  and  Exercise  (Female  C57  Only) 


Tissue  Mineral  Density  by  Diet  and  Exercise  (Male  C57  Only) 


1200 


1000 


800 


600 


400 


200 


0 


Tissue  Mineral  Density  by  Diet  and  Exercise  (Female  C3  Only) 


1200 


Tissue  Mineral  Density  by  Diet  and  Exercise  (Male  C3  Only) 


Yield  Load  (N)  Yield  Load  (N) 


Femur  Yield  Load  by  Diet  and  Exercise 


Yield  Load  by  Exercise  and  Diet  (Female  C57  Only) 


Yield  Load  by  Exercise  and  Diet  (Male  C57  Only) 


Ultimate  Load  (N)  Ultimate  Load  (N) 


Femur  Ultimate  Load  by  Diet  and  Exercise 


Ultimate  Load  by  Exercise  and  Diet  (Female  C57  Only)  Ultimate  Load  by  Exercise  and  Diet  (Male  C57  Only) 


Ultimate  Load  by  Exercise  and  Diet  (Female  C3  Only) 


Stiffness  (N/mm)  Stiffness  (N/mm) 


Femur  Stiffness  by  Diet  and  Exercise 


Stiffness  by  Exercise  and  Diet  (Female  C57  Only) 

600  -i - 


500 


400 


300 


Stiffness  by  Exercise  and  Diet  (Female  C3  Only) 

600  -i - 


500 


Stiffness  by  Exercise  and  Diet  (Male  C57  Only) 

600  -i - 


500 


400 


E 

E 


Stiffness  by  Exercise  and  Diet  (Male  C3  Only) 

600  -i - 


Displacement  Ratio  (mm/mm)  Displacement  Ratio  (mm/mm) 


Femur  Displacement  Ratio  by  Diet  and  Exercise 


Displacement  Ratio  by  Exercise  and  Diet  (Female  C57  Only) 


Displacement  Ratio  by  Exercise  and  Diet  (Male  C57  Only) 


Femur  Phosphate/Amide  1  Peak  by  Diet  and  Exercise 


Mineral  to  Matrix  Ratio  by  Exercise  and  Diet  (Female  C57  Only) 


Mineral  to  Matrix  Ratio  by  Exercise  and  Diet  (Female  C3  Only) 


Femur  Phosphate/Carbonate  Peak  by  Diet  and  Exercise 


Phosphate  to  Carbonate  Peak  by  Exercise  and  Diet  (Female  C57  Only) 


Phosphate  to  Carbonate  Peak  by  Exercise  and  Diet  (Female  C3  Only) 


(EvUiw/EvUiui)  A±  A3  (E»uiui/E»uiui)  A1  Aa 


Proximal  Vertebral  Bone  Volume  Fraction  by  Diet  and  Exercise 

Bone  Volume  Fraction  by  Diet  and  Excercise  (Female  C57  Only)  Bone  Volume  Fraction  by  Diet  and  Excercise  (Male  C57  Only) 

0.6  t -  0.6  t - 


Bone  Volume  Fraction  by  Diet  and  Excercise  (Female  C3  Only)  Bone  Volume  Fraction  by  Diet  and  Excercise  (Male  C3  Only) 


TB_N  (1/mm)  TB_N  (1/mm) 


Proximal  Vertebral  Trabecular  Number  by  Diet  and  Exercise 


Trabecular  Number  by  Diet  and  Exercise  (Female  C57  Only) 


Trabecular  Number  by  Diet  and  Exercise  (Male  C57  Only) 


T rabecular  Number  by  Diet  and  Exercise  (Female  C3  Only)  T rabecular  Number  by  Diet  and  Exercise  (Male  C3  Only) 


Normal 

Yes 


High 


Normal 

No 


Proximal  Vertebral  Trabecular  Thickness  by  Diet  and  Exercise 

T rabecular  Thickness  by  Diet  and  Exercise  (Female  C57  Only)  T rabecular  Thickness  by  Diet  and  Exercise  (Male  C57  Only) 


0.12 


0.08 


E 

E 


0.06 


0 


Low 


Normal  High 

Yes 


Low  Normal 

No 


Trabecular  Thickness  by  Diet  and  Exercise  (Female  C3  Only) 


T rabecular  Thickness  by  Diet  and  Exercise  (Male  C3  Only) 


0.12  -  0.12 


Proximal  Vertebral  Degree  of  Anisotropy  by  Diet  and  Exercise 

Degree  of  Anisotropy  by  Diet  and  Exercise  (Female  C57  Only)  Degree  of  Anisotropy  by  Diet  and  Exercise  (Male  C57  Only) 


2.5 


t 

Low  Normal 

Yes 

II 

High  Low  Normal 

No 

1 1  ii  i 

High  Low  Normal  High  Low  Normal  High 

Yes  No 

Degree  of  Anisotropy  by  Diet  and  Exercise  (Female  C3  Only)  Degree  of  Anisotropy  by  Diet  and  Exercise  (Male  C3  Only) 


2.5 


Low  Normal 

Yes 


High 


Low 


I 


Normal 

No 


2 


Normal 

Yes 


I 


Normal 

No 


BMD  (mg/cc)  BMD  (mg/cc) 


Proximal  Vertebral  Bone  Mineral  Density  by 


Bone  Mineral  Density  by  Diet  and  Exercise  (Female  C57  Only) 


Diet  and  Exercise 


Bone  Mineral  Density  by  Diet  and  Exercise  (Male  C57  Only) 


Bone  Mineral  Density  by  Diet  and  Exercise  (Female  C3  Only)  Bone  Mineral  Density  by  Diet  and  Exercise  (Male  C3  Only) 


Proximal  Vertebral  Tissue  Mineral  Density  by  Diet  and  Exercise 


Tissue  Mineral  Density  by  Diet  and  Exercise  (Female  C57  Only) 


Tissue  Mineral  Density  by  Diet  and  Exercise  (Male  C57  Only) 


Stiffness  (N/mm)  Stiffness  (N/mm) 


Vertebral  Stiffness  by  Diet  and  Exercise 

Stiffness  by  Exercise  and  Diet  (Female  C57  Only)  Stiffness  by  Exercise  and  Diet  (Male  C57  Only) 


Yield  Load  (N)  Yield  Load  (N) 


Vertebral  Yield  Load  by  Diet  and  Exercise 

Yield  Load  by  Exercise  and  Diet  (Female  C57  Only)  Yield  Load  by  Exercise  and  Diet  (Male  C57  Only) 


Ultimate  Load  (N)  Ultimate  Load  (N) 


Vertebral  Ultimate  Load  by  Diet  and  Exercise 

Ultimate  Load  by  Exercise  and  Diet  (Female  C57  Only)  Ultimate  Load  by  Exercise  and  Diet  (Male  C57  Only) 


Dratio  (mm/mm)  Dratio  (mm/mm) 


Vertebral  Displacement  Ratio  by  Diet  and  Exercise 


Displacement  Ratio  by  Exercise  and  Diet  (Female  C57  Only) 


Displacement  Ratio  by  Exercise  and  Diet  (Male  C57  Only) 
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Displacement  Ratio  by  Exercise  and  Diet  (Female  C3  Only) 


Displacement  Ratio  by  Exercise  and  Diet  (Male  C3  Only) 


ii 

1  i 

Low  |  Normal  |  High 

Yes 


Low  |  Normal  |  High 

No 


I  1 
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mi 

At 

Low  |  Normal 

Yes 

High 

Low  |  Normal  |  High 

No 

Tibial  Mid-diaphyseal  Cortical  Parameters  by  Sex  and  Exercise 


Cross  Sectional  Area  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Cortical  Thickness  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


0.25 


Female 


Male 


Bending  Moment  of  Inertia  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Shaping  Factor  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


(£uiui/6ui)  QiAia 


Tibial  Mid-diaphyseal  Cortical  Parameters  by  Sex  and  Exercise 


Bone  Mineral  Density  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Volumetric  Bone  Mineral  Density  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 
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1540 
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1500 
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Tibial  Proximal  Metaphyseal  Trabecular  Parameters  by  Sex  and  Exercise 


Tibia  Metaphyseal  Bone  Volume  Fraction  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Tibia  Metaphyseal  Trabecular  Number  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Tibia  Metaphyseal  Trabecular  Thickness  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Tibia  Metaphyseal  Trabecular  Spacing  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


0.16 


Tibial  Proximal  Metaphyseal  Trabecular  Parameters  by  Sex  and  Exercise 


Tibia  Metaphyseal  Bone  Mineral  Density  by  Sex  and  Exercise  Tibia  Metaphyseal  Tissue  Mineral  Density  by  Sex  and  Exercise 

(Normal  Diet,  C57  Strain)  (Normal  Diet,  C57  Strain) 


Twist  (deg)  Stiffness  (N-mm/deg) 


Tibial  Torsion  Parameters  by  Sex  and  Exercise 


Stiffness  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Failure  Torque  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Failure  Twist  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


Energy  by  Sex  and  Exercise 
(Normal  Diet,  C57  Strain) 


2T  d3 


M  C57  Low  Yes 


128T  d3 


M  C57  High  Yes 


% 


184T 

^3  Low  Yes 


198T  dl 
M  C3  Low  No 


226T  ml 

F  C3  Normal  Yes 


261T  dl 

F  C3  Normal  No 


296T  d3 
M  C3  High  Yes 


31 7T  d3 
M  C57  High  No 


366T  d3 
F  C57  Low  No 


41 5T  d3 

M  C57  Normal  Yes 


443T  dl 

M  C57  Normal  No 


478T  d3 
F  C57  High  No 


